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OBJECTIVES: 


Adhesively  bonded  joints  have  potential  for  wide  usage  in  a  variety  of  structural  applications. 
In  many  present  day  applications,  conventional  fasteners  such  as  bolts,  rivets,  welds  etc.,  are  found 
to  be  unsuitable,  particularly  when  the  components  are  made  of  polymeric  or  composite  material. 
Penetration  methods  (i.e.  drilling  holes,  etc.)  give  rise  to  high  stress  concentrations  and,  in  the  case 
of  composites,  sever  the  fiber  reinforcement,  causing  a  reduction  in  joint  strength.  These  problems 
can  be  overcome  in  a  simple  and  inexpensive  way  by  replacing  conventional  fasteners  with  adhesive 
bonds.  Unfortunately,  the  promise  of  adhesive  bonding  has  existed  for  over  a  decade,  with  slow 
progress  being  made  towards  their  use  in  engineering  designs.  While  adhesive  bonds  are  being  used, 
they  are  largely  confined  to  low-stress  applications.  They  have  yet  to  find  their  place  in  highly 
stressed  primary  structural  applications  in  aerospace,  automobile,  composites  or  similar  industries 
wherein  lie  their  huge  potential. 

Why  has  industry  not  fully  exploited  the  potential  of  adhesive  bonding?  One  reason  is  that  their 
integrity  over  the  lifetime  of  a  structure  -  say  5  to  20  years  of  service  -  cannot  yet  be  predicted  in  a 
satisfactory  manner.  This  problem  can  be  traced  to  the  fact  that  adhesive  bonds  have  properties 
that  are  time  dependent.  Adhesives  are  known  to  creep  in  service  and  hostile  environmental  pa¬ 
rameters,  such  as  temperature  and  moisture,  tend  to  accelerate  the  degradation  process,  making  it 
rather  difficult  to  predict  the  structural  integrity  in  a  design  lifetime. 


OBJECTIVES: 


1 


Moisture  induced  failure  of  polymer-to-metal  interfaces  is  a  severe  problem  in  applications 
involving  coatings,  sealants  and  adhesives.  Such  effects  are  particularly  pronounced  when  the 
adherend  and/or  the  adhesive  are  subject  to  moisture-induced  degradation,  and  when  the  displace¬ 
ment  of  the  bond-line  by  water  is  energetically  favorable.  Thermodynamics  of  surface  energies  and 
work  of  adhesion  has  proved  quite  successful  for  predicting  the  ultimate  stability  or  lack  of  stability 
of  a  joint.  Unfortunately,  these  approaches  have  not  been  able  to  incorporate  kinetics  to  provide 
predictions  for  the  time  to  failure.  A  significant  challenge  is  to  develop  a  rational  basis  to  combine 
thermodynamics,  mechanics  and  transport  kinetics  into  a  tool  for  an  understanding  of  bond  du¬ 
rability.  Since  the  problem  of  durability  is  strongly  interdisciplinary  in  nature,  a  scheme  is  needed 
in  which  ideas  from  all  the  key  disciplines  involved  arc  combined  in  an  integrated  computational 
tool.  Such  a  scheme  would  provide  a  predictive  capability  by  which  the  results  of  small  scale  tests 
could  be  used  to  determine  the  lifetime  of  a  given  joint.  A  need  therefore  exists  for  improved  ex¬ 
perimental  characterization  techniques  and  for  advanced  computational  procedures.  The  successful 
conclusion  of  such  a  research  program  would  give  the  designer  and  manufacturers:  a  better  char¬ 
acterization  of  their  products,  the  possibility  of  enlarging  markets  by  developing  better  user  confi¬ 
dence  in  structural  adhesives,  and  lowered  costs  through  reduction  in  unrealistic  safety  factors. 

The  following  generic  research  objectives  are  proposed: 

( 1 )  Develop  or  improve  phenomenological  descriptions  of  the  mechanical  response  of  adhesives, 
and  of  moisture  diffusion  in  adhesives. 

(2)  Develop  adequate  characterization  techniques  for  the  mechanical  behavior  of  the  adhesive,  for 
moisture  diffusion  along  the  bond  line  and  for  environmentally-induced  failure. 

(3)  Incorporate  the  resulting  predictive  capabilities  into  a  well  documented  computational  tool. 

This  dissertation  addresses,  in  varying  degrees  of  depth,  the  three  suggested  research  avenues. 
As  we  enter  more  detailed  theoretical  discussions,  an  effort  will  be  made  not  to  lose  sight  of  the  final 
goal,  viz:  providing  an  engineering  solution  to  the  durability  problem.  This  means  that  although 
an  effort  will  be  made  to  understand  mechanisms  at  the  molecular  level,  a  lack  of  complete  under¬ 
standing  of  a  fundamental  question  will  not  deter  us  from  describing  the  material  response  in  a 
mathematical  manner.  Following,  is  a  brief  outline  of  the  proposed  study  of  the  durability  problem: 
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We  first  observe  that  moisture  plays  a  central  role  in  the  degradation  process  of  adhesive 
joints.  Obviously,  there  are  two  main  rate  limiting  parameters:  the  transport  rate  of  moisture  and 
the  rate  of  environmental  degradation.  For  this  reason,  a  durability  study  should  naturally  proceed 
along  two  different  directions:  (1)  the  study  of  moisture  diffusion  and  (2)  the  determination  of  a 
failure  criterion.  While  a  general  theoretical  approach  can  be  developed  for  the  former,  no  general 
method  exists  for  the  latter  because  failure  mechanisms  arc  too  specific  to  a  given  adherend- 
interphase-adhcsive  system.  Therefore,  a  general  theoretical  approach  will  be  proposed  for  the 
diffusion  problem,  while  an  empirical  approach  will  be  favored  for  the  failure  study. 

'ITie  underlying  idea  for  modelling  diffusion  will  be  based  on  the  free  volume  concept,  i.e.,  that 
the  motion  of  diffusing  molecules  depends  on  the  space  available  to  them  in  the  form  of  vacancies 
between  polymer  chains.  Since  this  free  volume  is  also  known  to  control  the  rate  of  creep  or  re¬ 
laxation,  a  unified  mechanical-diffusional  theory  can  be  used.  While  our  reasoning  initially  will 
occur  at  the  molecular  level,  the  proposed  model  will  be  entirely  phenomenological,  mainly  to  re¬ 
duce  material  characterization  to  an  engineering  level.  The  problem  of  the  transport  rate  will  be 
addressed  by  adapting  and  enhancing  existing  diffusion  and  viscoelasticity  theories  based  on  free 
volume  concepts. 

The  second  stage  in  the  study  of  durability,  involves  the  loss  of  mechanical  strength.  It  has 
been  established  that  there  are  two  main  failure  modes  in  adhesive  joints:  cohesive  (inside  the  ad¬ 
hesive)  and  interfacial  (along  the  boundary  between  the  adhesive  and  the  adherend).  Whereas  cri¬ 
teria  for  cohesive  failure  are  already  available,  much  work  remains  to  be  done  on  the  problem  of 
interfacial  failure.  The  difficulty  arises  from  the  fact  that  interfacial  failures  usually  occur  in  moist 
environment  as  a  result  of  a  chemical  or  physical  attack  on  the  interphase.  Since  these  degradation 
phenomena  are  difficult  to  represent  analytically,  empirical  methods  are  often  used.  One  such 
method,  based  on  the  critical  concentration  concept  will  be  reviewed  and  tested  experimentally. 

The  limit  of  the  general  methodology  proposed  above,  is  found  in  the  case  of  joints  exposed 
to  unusually  harsh  conditions.  This  will  be  illustrated  with  the  case  of  rubber-to-stecl  bonds 
undergoing  cathodic  debonding  in  seawater.  An  alternate  methodology  will  be  proposed  and  a 
description  of  a  durability  specimen  will  be  given. 
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1.  A  MODEL  FOR  DIFFUSION  IN 
POLYMERS.  GOVERNING  EQUATIONS: 


Introduction: 


Moisture  may  enter  a  bond  by  diffusion  through  the  adhesive,  through  the  adherends  or  by 
moving  along  an  interphase  region.  If  one  assumes  that  water  molecules  have  the  same  degradation 
effect  on  the  adhesive  and  interphase  region,  regardless  of  their  migration  path,  then  it  becomes  of 
ultimate  concern  to  be  able  to  understand  and  predict  the  kinetics  of  moisture  ingression  via  the 
different  paths.  This  chapter  investigates  an  important  mode  of  moisture  intrusion  in  a  joint:  bulk 
diffusion  through  the  adhesive. 

Earlier  work  has  demonstrated  that  the  diffusion  rate  in  a  polymer  can  be  increased  by  a 
number  of  factors  such  as  temperature,  stress  and  damage.  In  particular,  "non-Fickian'  kinetics 
often  associated  with  a  sharp  diffusion  front,  has  been  attributed  to  a  strong  dependence  of  the 
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diffusion  coefficient  on  concentration  and/or  moisture  induced  damage.  Externally  applied  stresses, 
residual  stresses,  and  swelling  stresses  are  known  to  affect  diffusion  rate.  These  various  features 
(with  the  exception  of  damage)  are  reviewed  in  light  of  their  usefulness  to  the  specific  problem  of 
moisture  diffusion  in  adhesive  joints.  Although  damage  can  significantly  enhance  moisture 
ingression,  it  is  reasonable  to  exclude  it  from  this  preliminary  investigation  because  joints  in  service 
in  moist  environments  are  normally  loaded  to  a  small  fraction  of  their  ultimate  breaking  strength 
(dry),  thereby  minimizing  damage  arising  solely  from  the  applied  load. 

In  Chapter  1,  the  theoretical  background  leading  to  the  diffusion  governing  equations,  will 
be  reviewed  in  detail.  Although  the  key  concepts  underlying  the  proposed  theoretical  treatment  have 
been  studied  by  several  other  workers,  some  novelty  is  introduced  in  the  mathematical  formulation 
of  the  constitutive  relations.  One  of  the  major  objectives  of  this  study  is  to  offer  a  treatment  uni¬ 
fying  both  the  diffusion  behavior  and  mechanical  behavior  into  a  single  phenomenological  model. 
It  will  be  shown  that  the  unified  approach  not  only  facilitates  the  description  of  the  coupling  be¬ 
tween  diffusion  and  stress  relaxation,  but  also  leads  to  a  common  material  characterization.  Many 
workers  have  recognized  the  importance,  in  the  case  of  polymers,  of  the  interdependence  between 
the  diffusion  process  and  stress  relaxation.  Unfortunately,  the  solutions  proposed  to  date  have  either 
dealt  with  specific  boundary  value  problems,  or  have  failed  to  include  in  one  single  treatment,  all 
the  known  forms  of  coupling.  It  follows  that  the  solutions  offered  were  often  severely  restricted  in 
their  conditions  of  applicability.  By  contrast,  this  treatment  is  meant  to  be  as  general  as  possible, 
and  to  be  implemented  in  a  finite  element  code.  Unlike  closed  form  solutions,  the  Finite  Element 
Method  can  deal  with  almost  any  kind  of  boundary  value  problem  of  interest  to  designers.  In  the 
field  of  adhesion  engineering,  our  goal  is  to  provide  designers  with  a  general  computational  method 
to  predict  moisture  intrusion  into  adhesive  bond  lines  or  into  composite  matrices,  when  they  are 
subjected  to  a  combination  of  external  loading  and  internal  swelling. 

Chapter  1  is  introductory  and  thus  should  be  appraised  only  in  the  context  of  the  entire  work. 
In  particular,  Chapter  1  is  not  sufficient  to  fully  understand  the  effect  of  stress  and  stress  relaxation 
on  diffusion.  The  constitutive  relations  for  the  mechanical  behavior  and  swelling  behavior  are  dis¬ 
cussed  in  Chapter  3. 
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General  Form  of  the  Diffusion  Governing  Equation: 


Let  us  consider  a  binary  mixture  consisting  of  a  penetrant  of  small  size  (e.g.:  water)  perme¬ 
ating  a  polymer.  The  polymer  is  isotropic,  amorphous  and  can  be  in  either  the  glassy,  leathery  or 
rubbery  state.  It  is  natural  to  assume  that  the  small  penetrant  molecules  are  much  more  mobile 
than  the  polymeric  macromoleculcs.  It  follows  that  mass  transport  in  this  system  can  be  described 
in  terms  of  one  flux  only,  namely  that  of  the  penetrant,  because  the  polymer  can  be  considered  as 
a  fixed  reference.  In  general,  the  mass  flux  of  penetrant  can  be  separated  in  four  components,  each 
corresponding  to  a  class  of  thermodynamic  driving  force1  .  Each  driving  force  gives  rise  to  small 
perturbations  in  the  random  walk  of  the  diffusing  species.  Assuming  the  driving  torces  are  uncou¬ 
pled,  and  the  departure  from  randomness  small,  the  resulting  macroscopic  fluxes  are  additive: 

j  =  jc  +  jr+.r  +  j*  (i) 

where:  Jc  =  Fickian  mass  flux  (entropy  driven) 

JT  =  thermally  driven  mass  flux 

J"  =  stress  induced  mass  flux 

J*  =  forced  mass  flux  (driven  by  an  external  force) 

(Bold  characters  will  represent  vectorial  quantities  throughout  this  dissertation) 

The  Fickian  flux  arises  as  a  result  of  concentration  gradients  and  is  given  by: 

JC  =  -  DVC  (2) 

where:  D  =  diffusion  coefficient 

C  =  penetrant  concentration 
V  =  gradient  operator 
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Expression  (2)  is  known  as  Fick's  First  Law  as  long  as  D  is  a  constant  or  a  function  of  C  only. 
It  will  be  shown  later  that  the  diffusion  coefficient  is  a  function  of  temperature  and  stress  (or  strain) 
as  well: 

D  =  D(C,  T,  a)  (3) 

where:  T  =  temperature 

a  =  stress. 

The  thermal  flux  is  driven  by  temperature  gradients  and  is  known  as  the  Soret  Effect1: 

t  DtC 

J =~RT  VLnT  (4) 

where:  DT  =  thermal  diffusion  coefficient 

Thermal  diffusion  becomes  a  significant  component  of  the  overall  flux  when  an  adhesive 
undergoes  frequent  temperature  fluctuations.  Sudden  variations  in  external  temperature  give  rise 
to  temperature  gradients  internally  which  may  persist  for  a  long  time,  especially  when  the  structure 
has  a  low  thermal  conductivity. 

The  stress  induced  flux  is  produced  by  a  non-uniform  stress  field2  and  is  given  by: 

J *  =  -Jj£rVWp  (5) 

where:  R  =  gas  constant 

Wp  —  potential  energy  function 

There  is  no  general  form  for  WP.  It  is  a  function  of  the  local  entropy  change  of  the  penetrant 
when  the  polymer  is  subject  to  a  stress  (strain)  field.  Since  the  stress  distribution  within  a  joint  is 
rarely  uniform,  proper  attention  should  be  given  to  this  effect  a  priori.  Furthermore,  if  the  adhesive 
tends  to  swell  by  a  large  amount  in  the  presence  of  the  penetrant,  J'  could  conceivably  become 
quite  large  because  differential  swelling  across  the  adhesive  is  generally  constrained  by  stiff 
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adherends,  leading  to  large  pressure  gradients  internally,  under  either  sorption  or  desorption  proc¬ 
esses.  A  possible  form  for  Wp  will  be  suggested  in  this  chapter. 


The  forced  diffusion  flux  is  of  importance  in  cases  when  the  penetrant  molecules  are  subject 
to  an  external  force  field  (e.g.,  electrical  or  inertial).  In  an  adhesive  joint  such  a  field  may  exist  lo¬ 
cally  in  the  region  known  as  the  interphase.  Close  to  a  hydrophilic  oxide  layer,  the  intermolecular 
field  experienced  by  a  water  molecule  can  be  more  intense  than  in  the  bulk  polymer.  This  disparity 
is  expected  to  cause  local  perturbations  in  the  diffusive  flux.  It  should  be  noted  that,  in  this  article, 
J*  refers  to  a  component  of  the  flux  arising  from  any  external  force  field,  with  the  exception  of  stress 
fields.  The  ufcct  of  stress  fields  is  already  incorporated  in  the  J°  component  of  the  flux. 

Let  us  now  consider  the  total  flux  J.  Since  conservation  of  mass  must  be  satisfied,  we  must 

have: 


~  =  -  V.J 
dt 


(6) 


where  t  is  time,  leading  to  the  final  governing  equation: 

^-  =  V.|d(C,T,<7)[vC  +  -^t  VH'/,]  +  Jr+J*j  (7) 

In  the  remaining  discussion,  the  problem  will  be  reduced  to  that  of  the  diffusion  in  a  medium 
of  uniform  temperature  and  in  the  absence  of  an  external  force  field  (Jr  =  J*  =  0).  Expression  (7) 
shows  that  the  effect  of  stress  is  twofold:  (1)  stress  field  gradients  give  rise  to  a  thermodynamic  force 
V  WP  and  (2)  stress  also  affects  the  mobility  of  the  penetrant,  quantified  here  by  the  diffusion  coef¬ 
ficient.  The  effect  of  stress  on  mobility  does  not  disappear  in  a  uniform  stress  field,  although  the 
stress  induced  driving  force  does  vanish. 

At  this  point,  it  is  useful  to  recall  that  polymeric  adhesives  are  viscoelastic  in  nature.  Thus, 
our  governing  equation  is  valid  instantaneously  only.  Stress  and  stress  gradients  need  to  be  updated 
constantly  in  order  to  properly  account  for  the  rheology  of  the  material.  The  problem  of  the  cou¬ 
pling  between  the  diffusion  boundary  value  problem  and  the  viscoelasticity  boundary  value  problem 
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will  be  addressed  in  more  detail  in  Chapter  3.  For  the  time  being,  we  just  need  to  k.  o  in  mind  that 
the  governing  equation  is  implicitly  time  dependent  on  the  right  hand  side  of  the  equation. 

A  thorough  literature  review2-7  led  us  to  conclude  that  expression  (7)  contains  many  of  the 
known  features  necessary  for  modelling  non-Fickian  transport  in  the  absence  of  damage.  These 
features  are: 

1.  the  history  dependence  of  the  diffusion  coefficient,  in  this  case,  through  the  coupling  with 
the  viscoelastic  response, 

2.  the  effect  of  a  non-uniform  stress  distribution,  and 

3.  the  concentration  dependence  of  the  diffusion  coefficient.  (Although  listed  as  a  source 
of  anomalous  behavior,  this  concentration  dependence  does  not  in  fact  violate  Fick's 
Second  Law  in  its  most  general  form). 

Note  that  the  proposed  governing  equation  is  by  no  means  the  only  possible  theoretical  model. 
It  just  contains  provision  for  important  features  often  found  in  separate  theories,  but  not  in  a  uni¬ 
fied  approach. 

The  stress  dependence  of  the  solubility  (accompanied  here  by  an  implicit  time  dependence) 
has  also  been  identified  as  a  source  of  non-Fickian  behavior7-9.  This  difficulty  can  be  addressed  in 
the  framework  of  a  numerical  analysis,  by  normalizing  the  concentration  with  respect  to  the  satu¬ 
ration  level  at  the  current  time  step,  and  by  providing  an  adequate  model  for  the  stress  dependence 
of  solubility.  Such  a  model  will  be  discussed  later  in  this  chapter. 
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A  Model  for  the  Diffusion  Coefficient  Based  on  the 
Cohen-  Turnbull  Theory: 


The  purpose  of  this  section  is  to  establish  a  general  theory  for  predicting  the  diffusion  coeffi¬ 
cient  of  small  molecules  in  polymeric  materials.  Particular  attention  will  be  devoted  to  the  tem¬ 
perature  dependence,  stress  dependence  and  penetrant  concentration  dependence  of  D  . 


The  Cohen-Turnbull  Model: 


The  Cohen-Turnbull  model1011  was  originally  developed  to  describe  self-diffusion  in  an  ideal 
liquid  made  of  hard  spheres,  but  it  has  been  extended  to  concentrated  solutions  and  undiluted 
polymers.  The  model  is  based  on  the  idea  that  molecular  transport  occurs  by  the  movement  of 
molecules  into  voids,  with  a  size  greater  than  some  critical  value.  Voids  are  formed  by  the  statistical 
redistribution  of  the  free  volume.  Free  volume  is  defined  as: 

Yf=  V-  V0  (8) 


where:  V  =  specific  volume 

V„  =  occupied  volume 

Furthermore,  it  is  assumed  that: 

1.  Free  volume  is  continually  redistributed  with  time,  and  that  no  local  free  energy  is  required  for 
redistribution.  As  a  result,  distribution  is  random. 

2.  Molecular  transport  occurs  by  the  movement  of  diffusing  molecules  into  voids  of  at  least  their 
molecular  size,  which  are  formed  as  a  result  of  the  random  redistribution  of  free  volume. 


I.  A  MODEL  FOR  DIFFUSION  IN  POLYMERS.  GOVERNING  EQUATIONS: 


10 


In  the  case  of  the  self  diffusion  (or  viscous  flow)  of  a  polymer,  the  molecular  jump  represen¬ 
tation  must  be  replaced  by  the  concept  of  segmental  rearrangements  which  can  occur  whenever  the 
surrounding  free  volume  reaches  some  critical  value  Vc.  In  other  words,  in  order  for  a  portion  of 
the  polymeric  chain  to  take  a  new  conformation,  the  "cage"  formed  by  nearest  neighbors  must  reach 
a  critical  size. 

The  probability  that  n  increments  of  free  volume  of  average  size  Vf  would  accumulate  into  a 

V 

void  of  size  Ve  is  proportional  to  exp(  —  9  ),  where  0  is  an  overlap  factor  arising  from  the  fact 

V, 

that  the  same  void  is  available  to  more  than  one  molecule.  In  the  case  of  viscous  flow,  critical 
volume  Vc  must  be  large  enough  to  allow  a  local  rearrangement  of  the  chain  conformation,  whereas 
in  the  case  of  the  diffusion  of  a  small  penetrant,  Vc  must  be  large  enough  to  allow  a  local  molecular 
jump  of  the  penetrant.  Thus  Vc  is  a  characteristic  of  the  polymer  only  in  the  case  of  self  diffusion 
(or  viscous  processes)  and  Ve  is  a  characteristic  of  a  given  polymer-penetrant  pair  in  the  case  of  the 
diffusion  of  a  small  penetrant  molecule.  From  the  form  of  the  probability  distribution,  it  follows 
that  the  diffusion  coefficient  contains  an  activation  volume  term: 


D  =  D0  exp(  —  6  -y- ) 


Equation  (9)  implies  that  free  volume  is  the  only  parameter  involved  in  describing  transport 
phenomena.  This  theory  is  not  sufficient  however,  at  temperatures  too  low  relative  to  the  glass 
transition  temperature  T(,  when  motions  of  the  chains  within  the  free  volume  are  too  slow.  Nor 
does  it  apply  at  very  high  temperatures  where  activation  energy  barriers  must  be  overcome13 .  Note 
that  in  the  framework  of  the  Cohen-Tumbull  theory  the  free  volume  is  not  the  true  free  volume  in 
the  geometric  sense,  but  only  that  portion  of  the  total  free  volume  which  can  be  redistributed  with 
no  local  change  in  free  energy.  It  follows  that  the  occupied  volume  V„  defmed  here  is  not  the  hard 
shell  volume,  but  a  larger  quantity  which  can  be  changed  by  a  stress  field  (i.e.  the  occupied  volume 
is  compressible). 
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The  Macedo-Litovitz  and  Vrentas-Duda  Theories: 


Macedo  and  Litovitz12  introduced  an  expression  for  the  self-diffusion  coefficient  very  similar 
to  that  of  Cohen  and  Turnbull,  except  that  it  contains  an  activation  energy  factor: 


D  =  D0  exp 


(10) 


It  was  postulated  that  in  addition  to  the  minimum  hole-size  requirement,  an  energy  barrier  must 
be  surmounted.  When  inverted,  expression  (10)  describes  a  viscosity.  The  obtained  form  is  virtually 
identical  to  a  known  extension  of  the  WLF  equation  including  an  Arrhenius-like  temperature 
dependence13.  In  the  glassy  state  (non-equilibrium),  the  temperature  dependence  predicted  by  the 
elementary  WLF  equation  is  known  to  be  incorrect  and  an  Arrhenius  temperature  dependence  has 
been  widely  reported.  This  finding  is  consistent  with  expression  (10)  and  will  be  the  object  of  fur¬ 
ther  discussion  later. 

Vrentas  et  a/.14  developed  a  constitutive  behavior  for  the  self  diffusion  in  polymers  based  on 
equation  (10).  Their  treatment  of  the  free  volume  is  more  complex  than  the  one  used  in  this  work. 
In  addition  to  the  free  volume  treatment  of  Cohen  and  Turnbull,  the  Flory-Huggins  polymer  sol¬ 
ution  theory  and  some  aspects  of  the  entanglement  theory  of  Bueche  are  employed.  An  average 
free  volume  Vf  as  well  as  an  average  critical  volume  Vc  are  defined  for  the  binary  mixture  sorbent 
(subscript  1)  -  polymer  (subscript  2): 


Vf-  &ua>](&2i  +  T-  Tgi)  -I-  knw2(k22  +  T—  Tg2) 


Vc  =  “>\V\  +(02iV2 


where: 


co,  =  mass  fraction  of  component  i; 

Tfl  =  glass  transition  temperature  of  component  i; 

ku,kj\  —  free  volume  parameters  for  the  solvent; 

kuk22  =  free  volume  parameters  for  the  polymer. 


(ID 

(12) 
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t,  =  ratio  of  the  critical  molar  volume  of  the  penetrant  jumping 

unit  to  the  critical  molar  volume  of  the  polymer  jumping  unit 
V]  =  specific  critical  hole  free  volume  of  component  i. 

In  the  limit,  when  the  mass  fraction  of  penetrant  is  small,  l^and  Vc  become  mainly  polymer 
dependent  and  the  above  approach  becomes  equivalent  to  the  or.c  shown  herein.  Our  simpler  ap¬ 
proach,  in  which  the  free  volume  of  the  mixture  is  dominated  by  that  of  the  polymer  is  very  rea¬ 
sonable  because  polymers  with  a  large  moisture  solubility  are  not  suitable  for  the  structural  adhesive 
applications  we  are  concerned  with. 


Extension  of  the  Doolittle  Theory  to  the  Problem  of  Diffusion  of  Small 
Molecules: 


Doolittle  proposed  an  empirical  equation  relating  polymer  viscosity  to  the  free  volume,  which 
was  found  to  apply  to  polymers  in  their  rubbery  range13.  Fluidity  F,  the  inverse  of  viscosity  is  given 

by: 


F=  A~] 


exp(  -  B- 


V- Vr 


(13) 


where  A  and  B  are  empirical  constants.  The  parallelism  between  equations  (9)  and  (13)  is  striking. 
This  was  explained  earlier  by  the  fact  that  the  elementary  process  for  viscous  flow  and  the  elemen¬ 
tary  process  for  sorbent  diffusion  arc  identical.  In  order  to  prove  that  equations  (9)  and  (13)  are 
identical  mathematically,  one  can  use  1  instcin's  equation  for  the  diffusion  coefficient  and  Stokes' 
equation  for  the  friction  constant.  According  to  Finstein's  equation  for  the  diffusion  coefficient15: 
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D  = 


kT 

fn 


(14) 


where: 


D  =  diffusion  coefficient 
k  =  Boltzman  constant 

T  =  temperature 
f  =  friction  constant 


Equation  (14)  holds  for  any  thermally  driven  process  governed  by  the  three-dimensional  random 
walk  of  diffusing  species.  Stokes'  equation  states  that  the  friction  constant  is  proportional  to 
viscosity15  : 


fn  =  *1 


(15) 


As  reviewed  by  Ferry13,  the  work  of  Rouse  and  Bueche  demonstrated  that  equation  (15)  holds 
for  polymers  and  that  the  proportionality  coefficient  X  is  a  function  of  density,  molecular  weight, 
and  molecular  geometry13  .  In  this  case,  f„  is  a  measure  of  segmental  frictional  resistance.  Com¬ 
bining  (14)  and  (15)  yields: 

D  =  X~]kTF  (16) 

And  wc  find  a  diffusion  coefficient  of  the  form: 


D  =  ixT 


-6 


V'Vf\ 

Vf  ) 


(17) 


where  fj.  i:  a  material  constant.  Let  rj  and  rj,  be  the  viscosities  of  the  polymer, /and  f„  the  fractional 
free  volumes  at  a  given  temperature  T  and  a  reference  temperature  T„ ,  respectively.  (The  fractional 
free  volume  is  defined  as  the  ratio  VJV).  Doolittle,  in  his  work  on  the  viscous  flow  of  polymers 
showed  that13: 
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where  B  is  the  empirical  constant  mentioned  earlier.  Assuming  that  /  changes  linearly  with  tem¬ 


perature, 

f=fo  +  «oO'-T0)  (19) 

and  combining  with  expression  (18)  leads  to  the  well  known  WLF  equation  for  the  shift  factor  a: 

jr{T-T0) 

In  a=-  - -  (20) 

—  +  r  _  t 

The  WLF  equation  concerns  itself  with  the  viscoelastic  response.  It  was  noted  earlier  that  in 
general,  the  critical  hole  size  Vc  for  the  local  rearrangement  of  a  polymer  segment  is  not  the  same 
as  the  cntical  hole  size  for  penetrant  jump.  It  can  be  seen  in  expression  (9)  that  a  change  in  Vc  can 
be  expressed  mathematically  by  varying  0,  with  Ve  remaining  constant.  Noting  that  0  and  B  arc 
inversely  related,  it  is  now  possible  to  conveniently  extend  the  above  results  to  the  diffusion  of  a 
penetrant  of  small  size,  simply  by  replacing  B  by  BD,  a  numerical  parameter  inversely  related  to  the 
minimum  hole  size  for  the  jump  of  a  small  penetrant  molecule. 

By  combining  expressions  (14)  (15;  and  (18)  a  diffusion  counterpart  of  the  Doolittle  equation 
is  found: 

'>■4  (2,) 
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Derivation  of  the  Constitutive  Equation  for  the  Diffusion  of  a  Penetrant  of 
Small  Size: 


Several  investigations  in  the  field  of  rheology  have  suggested  that  free  volume  is  a  good  uni¬ 
fying  parameter  to  describe  changes  in  the  time  scale  of  material  response.  In  the  WLF  equation, 
the  effect  of  temperature  was  incorporated  by  stating  that  the  fractional  free  volume  was  a  linear 
function  of  temperature  in  the  rubbery  range.  Likewise,  the  effect  of  stress  and  solvent  concen¬ 
tration  can  be  introduced  by  stating  as  a  first  approximation  that  the  fractional  free  volume  is  a 
linear  function  of  the  dilatational  stress  component  and  the  solvent  concentration.  Knauss  and 
Emri16  used  this  concept  to  develop  a  nonlinear  viscoelastic  constitutive  behavior  centered  on  free 
volume  effects.  They  postulated  that  the  change  in  fractional  free  volume  due  to  each  variable  was 
additive: 


f—fo  +  ciAT+^  +  yC  (22) 

where:  a  -  coefficient  (volumetric)  of  thermal  expansion  of  the  free  volume 

y  =  coefficient  (volumetric)  of  swelling 

=  volume  dilatation  of  the  free  volume  due 
to  external  loads 

However,  as  pointed  out  by  Knauss  and  Emri  it  is  more  appropriate  to  state  that  the  free  volume 
depends  on  the  temperature  history,  strain  history  and  swelling  history  if  we  are  to  extend  this 
concept  to  the  glassy  range.  Thus  / should  be  written  in  the  form  of  a  sum  of  convolution  integrals: 

f=fo  +  «(0  *dT  +  4(0  +  y(0  *dC  (23) 

At  first,  we  will  limit  ourselves  to  the  simpler  case  where  a,  M  and  y  are  time  independent  and 
equation  (23)  reduces  to  (22).  These  assumptions  are  known  to  be  correct  above  Tt  where  the 
system  is  always  close  to  thermodynamic  equilibrium.  The  linearity  in  concentration  is  unlikely  to 
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be  obeyed  when  strong  diluents  are  used  and  high  concentrations  are  approached.  For  this  reason, 
it  will  be  assumed  that  the  change  in  free  volume  due  to  the  solvent  is  better  described  by: 


A/=  yCN  (24) 

And  (22)  then  becomes: 

f=fo  +  ct&T+/kk  +  yCN  (25) 

The  exponent  N  must  be  considered  as  an  adjustable  parameter  relating  to  the  occurrence  of  non¬ 
zero  volume  change  on  mixing.  Substituting  (25)  in  (21)  yields  an  expression  for  the  diffusion  co¬ 
efficient: 


D  =  T  exp^ 


fo 


a(T—  Ta)  +  <4*  +  yCN  1 

fo  +  a(7’“  Ta)  +  4/c  +  yCN  J 


(26) 


Although  there  are  eight  parameters  in  this  theory,  each  has  a  precise  physical  significance.  Sub¬ 
script  o,  on  any  variable  refers  to  a  reference  temperature  7'0.  In  order  to  use  BD,f,  and  a  from  the 
WLF  theory,  it  is  convenient  to  use  T,  =  Tg.  The  domain  of  validity  of  equation  (26)  is  the  same 
as  that  of  the  WLF  equation,  that  is,  roughly:  {Tg  Tg  +  100). 

As  pointed  out  oreviously,  an  activation  energy  must  be  introduced  if  one  wishes  to  extend 
(26)  above  Tt  +  100.  This  temperature  range  however,  is  not  considered  in  the  current  discussion. 


tension  to  the  Glassy  State,  Physical  Aging: 


Expression  (26)  was  originally  derived  for  small  penetrant  molecules  and  for  temperatures 
corresponding  to  the  rubbery  range  of  the  material.  Local  motions  of  the  main  chains  or  of  the  side 
groups,  as  well  as  brownian  motions  of  the  penetrant  molecules,  are  known  to  persist  in  the  glassy 
state  as  well.  A  list  of  possible  local  chain  motions  is  given  in  Reference  13.  Since  such  local  seg¬ 
mental  rearrangements  can  drive  the  diffusion  of  small  penetrant  molecules,  the  phenomenological 
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model  derived  earlier  still  applies  in  principle  below  the  glass  transition  temperature,  except  for  the 
temperature  dependence.  The  reason  the  temperature  term  is  not  valid  any  more  is  that  a  de¬ 
scribes  the  dilatation  of  the  free  volume  above  Tt  only,  (a  is  equal  to  the  difference  between  the 
volume  thermal  expansion  coefficient  of  the  specific  volume  and  that  of  the  occupied  volume). 
Below  Tv  the  fractional  free  volume  is  much  less  temperature  dependent  (although  it  is  still  time 
dependent  due  to  aging).  It  has  been  shown  by  Matsuoka  et  al.  that  the  abrupt  change  in  the 
temperature  dependence  at  the  Glass  Transition  can  be  explained  by  using  Adam-Gibbs'  formula 
for  the  relaxation  behavior17.  Adam-Gibbs'  formula  expresses  the  shift  fr.ctor  in  terms  of  the  con¬ 
figurational  entropy  of  the  polymer,  rather  than  in  terms  of  the  free  volume.  The  configurational 
entropy  is  used  as  a  measure  of  the  size  of  the  cooperatively  rearranging  region  and  seems  to  be  a 
better  description  of  the  state  of  the  polymer  below  Tt  than  the  free  volume.  This  entropy  based 
approach  leads  to  an  Arrhenius-like  equation  if  the  entropy  is  frozen  at  some  value,  as  in  the  rapidly 
quenched  glassy  state,  or  leads  to  a  Vogel-Fulcher-like  equation  (same  as  WLF)  above  Tt.  In  order 
to  remain  consistent  with  the  above  results  as  well  as  with  expression  (10),  the  following  form  is 
proposed  for  the  diffusion  coefficient  below  Tg: 


D  =  Dg  exp(  — 


AH 

RT 


) 


i kk  +  yCN  1 
fg  +  Jck  +  ycN  J 


(27) 


where:  AH  -  activation  energy 

R  =  the  gas  constant 


Expression  (27)  states  that  the  free  volume  continues  to  govern  the  effect  of  stress  and  swelling 
below  the  glass  transition  temperature,  implying  that  these  forces  act  primarily  on  very  local 
volumetric  properties  such  as  the  critical  void  size.  By  contrast  temperature  appears  to  affect  chain 
motions  over  larger  volumes  encompassing  cooperative  segmental  motions.  The  Arrhenius  tem¬ 
perature  dependence  shown  in  expression  (27)  has  been  observed  experimentally  on  numerous 
penetrant-polymer  systems  below  the  Glass  Transition  Temperature1.  Often,  the  Arrhenius  form  is 
also  used  above  Tt,  with  an  activation  energy  differing  from  that  below  Tt.  This  is  simply  due  to 
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the  fact  that  the  WLF  equation  can  be  mathematically  fitted  as  an  Arrhenius  equation  by  using  an 
apparent  activation  energy. 

Of  course, /is  not  the  same  at,  and  below  the  glass  transition  temperature.  Thus  expression 
(27)  by  itself  gives  an  incomplete  description  of  the  material  behavior.  This  apparent  inconsistency 
is  dealt  with  by  introducing  physical  aging  into  the  model.  We  will  see  that  the  use  of  a  shift  factor 
(acceleration  factor)  is  phenomenologically  equivalent  to  having  a  reduced  tfkk- 

In  the  rubbery  state,  equilibrium  is  reached  very  rapidly  in  response  to  variations  in  temper¬ 
ature,  stress,  and  penetrant  concentration.  This  is  due  to  the  fact  that  the  free  volume  response  is 
quasi-instantaneous  under  these  conditions.  By  contrast,  a  material  in  the  glassy  state  is  not  in 
thermodynamic  equilibrium,  and  the  response  of  the  free  volume  to  changes  in  external  conditions 
is  delayed.  For  this  reason,  a  time  dependent  diffusion  coefficient  containing  hereditary  integrals 
is  more  promising  for  the  glassy  range.  Expression  (27)  is  still  valid  instantaneously  however,  and 
can  be  used  in  an  iterative  numerical  scheme.  Glassy  polymers  have  a  volume  enthalpy  and  entropy 
which  arc  larger  that  they  would  be  in  the  equilibrium  state.  This  metastable  (or  supercooled)  state 
causes  the  free  volume  to  slowly  collapse  with  time  until  equilibrium  is  reached.  The  phenomenon 
is  known  as  physical  aging  "  .  Successful  diffusion  modelling  in  the  glassy  range  must  incorporate 
this  important  effect. 

Struik18  proposed  a  simple  method  to  incorporate  aging  in  existing  constitutive  equations. 
The  time  dependent  reduction  of  the  free  intermolecular  space  causes  relaxation  processes  to  take 
place  over  a  longer  time.  Thus,  an  acceleration  factor  was  defined  in  order  to  relate  actual  time  to 
effective  time: 


dX  =  a(t)dt 


(28) 


where: 


X  =  effective  time 

t  =  actual  time 

a{t)  =  acceleration  factor 


The  acceleration  factor  may  also  be  defined  in  terms  of  molecular  mobility  M(t): 
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where: 


t,  =  aging  time  at  the  start  of  service  life  or  testing 

t  =  service  life  or  test  time 


Struik  proposed  the  following  power  law  for  the  acceleration  factor: 

“W  “  ( 1J+T  )*  <3°> 

where  fi  is  a  constant  and:  Q  <,  n  <,  1.  Expression  (30)  was  derived  theoretically  and  confirmed 
by  long  term  experiments18  .  We  conclude  that  expression  (23)  can  be  readily  modified  to  incor¬ 
porate  aging,  by  replacing  actual  time  t  by  effective  time  X  in  the  integral  symbols.  Moreover,  we 
have  established  from  the  above  discussion  that  X  is  related  to  t  by: 


At  this  stage,  it  must  be  noted  that  there  exists  a  secondary  transition  temperature  T„  in  the 
glassy  range,  below  which  physical  aging  ceases,  due  to  extremely  low  molecular  mobility.  It  follows 
that  the  above  treatment  no  longer  applies  below  T„.  Adamson19  has  proposed  both  a  conceptual 
model  and  an  experimental  method  to  estimate  the  free  volume  in  the  entire  glassy  region,  including 
the  truly  glassy  region  below  Tw.  Further,  his  data  confirm  that  the  Free  Volume  Theory  of  mo¬ 
lecular  mobility  -with  some  minor  revisions-  still  explains  observed  sorption  behavior  in  the  glassy 
domain.  It  is  true  that  the  physics  of  free  volume  is  basically  incorrect  in  a  non-equilibrium  state 
such  as  the  glassy  state.  It  all  remains  that  the  basic  concepts  of  the  theory  are  still  extremely  helpful 
in  deriving  correct  predictive  models  for  the  behavior  under  stress.  As  demonstrated  by  Knauss16  , 
in  the  field  of  Viscoelasticity,  these  models  can  be  valuable  tools  for  life  prediction,  in  spite  of  their 
inherent  imperfection. 
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A  Potential  Field  Function  for  the  Stress-Induced  Diffusion 
Flux: 


It  was  shown  earlier  that  stress  (strain)  can  affect  the  diffusion  rate  in  two  ways:  (1)  by  al¬ 
tering  the  mobility  of  the  penetrant  and  (2)  by  producing  an  extra  driving  force.  By  relating  the 
diffusion  coefficient  to  strain,  we  addressed  the  mobility  effect.  We  now  must  establish  a  theoretical 
background  to  estimate  the  magnitude  of  the  stress-induced  driving  force.  A  primary  cause  for  this 
driving  force  in  a  joint  is  uneven  swelling  in  the  direction  of  the  moisture  flux.  Uneven  swelling 
gives  rise  to  pressure  gradients  along  constrained  bond  lines.  In  contrast  to  fluids  where  pressure 
energy  is  dissipated,  viscoelastic  solids  (eg:  adhesives)  will  store  some  elastic  energy,  causing  the 
chemical  potential  of  the  penetrant  to  be  changed.  The  potential  energy  function  introduced  in 
expression  (5)  is  a  measure  of  this  free  energy  change. 

Various  forms  have  been  suggested  for  potential  field  function  Wr  by  workers  interested  in  the 
diffusion  of  interstitials  in  stressed  metals20121,  in  the  stress-induced  diffusion  of  macromolecules2122  , 
and  in  the  stress  assisted  diffusion  of  small  penetrants  in  polymers9.  In  each  case,  the  method 
consisted  of  finding  the  entropic  potential  of  the  diffusing  species  under  the  stress  conditions  of  in¬ 
terest.  For  the  problem  of  transport  of  small  penetrant  molecules  in  polymers,  we  will  take  the  free 
energy  of  mixing  of  a  sorbent  (subscript  1)  with  a  polymer  (subscript  2),  for  our  entropic  potential 
function.  From  the  Flory-Huggins  theory23,  we  know  that  the  chemical  potential  of  the  sorbent  is 
given  by: 

Mi  —  M?  +  RT  [xi  v|  +  Ln  Vj  +  v2]  (32) 

where:  v,  =  volume  fraction  of  penetrant 

Vj  =  volume  fraction  of  polymer 
X,  =  Flory-Huggins  interaction  parameter 
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Let  Vx  and  V2  be  the  volume  of  penetrant  and  polymer  respectively,  in  the  binary  mixture. 
The  fractional  free  volume  may  then  be  defined  as: 

Vf 

f~  ~7r~  (33) 

Using  the  same  notation  as  in  the  section  dealing  with  diffusivity,  and  assuming  that  the  penetrant 
tends  to  fill  the  available  free  volume,  the  volume  of  penetrant  may  be  approximated  by: 


V i  =  V2f 


(34) 


The  volume  fraction  of  polymer  under  stress  may  be  approximated  by: 

y  _ _ 1 _ 

1  +fo  +  *4*  +  yCN 


(35) 


Likewise,  the  volume  fraction  of  penetrant  under  stress  (before  the  system  has  recovered  an 
equilibrium  state)  may  approximated  by: 


f0  +  yc 


N 


1  +fo  +  ®Lt  +  yC; 


N 


(36) 


Since  the  formation  of  extra  "holes'  in  the  presence  of  mechanical  deformation  does  not  affect 
the  sum  of  all  pair-wise  interaction  energies  between  the  polymer  and  the  penetrant  molecules,  the 
enthalpy  term  in  expression  (32)  is  unchanged  under  stress.  Thus  the  chemical  potential  fi\  of  the 
penetrant  in  a  strained  polymer  is  given  by: 


n\  =  ri  +  RT 


1 


1  +fo  +  yC 


■tv 


+  Lnl 


fo  +  yc 


N 


^+^  +  4*  +  yC'v/  1  +.4  +  4fc  +  yC‘ 


H  37> 


The  excess  chemical  potential  of  the  penetrant  with  respect  to  the  unperturbed  state  is  given  by: 


t>A  =  RT\ 


Lnl 


1+  fo  +  yC 


i  +/>  +  4*  +  yCs  /  i+fo  +  c>kk  +  yCN  i+To  +  yC^J 


(38) 
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6n\  is  the  excess  entropic  potential  of  the  diffusing  species  under  stress,  and  can  be  taken  as  our 


potential  field  function  WP. 

Rearranging  expression  (38)  and  noting  that  f0< sc  1  ,  tfkk  c  1  and  yCN  •«:  1  in  all 

polymer-penetrant  systems  of  interest  in  this  study  (small  strains,  minimal  swelling),  it  immediately 
follows  that: 


VWP~  -RT  \/kk  (39) 

and  the  stress-induced  flux  is: 

J°c*  +  DC  \/kk  (40) 

Note  that  4k  is  related  to  the  trace  of  the  stress  tensor  by  an  instantaneous  compressibility.  Thus, 
to  some  extent,  \t/kk  is  analogous  to  a  pressure  gradient.  However,  in  contrast  to  pressure,  which 
can  only  be  positive,  strain  or  stress  may  be  positive  or  negative,  leading  to  an  increase  or  decrease 
of  the  chemical  potential.  Therefore,  the  excess  chemical  potential  defined  earlier  is  an  algebraic 
quantity. 

Consider  a  polymer  constrained  between  two  rigid  adherends  held  at  a  fixed  distance.  If  the 
polymer  tends  to  swell  in  the  presence  of  moisture,  which  is  the  mechanical  component  of  free 
volume  dilatation,  becomes  negative  and  more  so  as  the  concentration  of  penetrant  increases  (See 
chapter  3  for  the  constitutive  relations  describing  the  hygro-mechanical  behavior).  Thus  VC  and 
\e/kk  have  opposite  signs  and  the  two  corresponding  fluxes  are  of  the  same  sign,  leading  to  acceler¬ 
ated  moisture  penetration.  This  prediction  is  consistent  with  a  number  of  results  reported  in  Ref¬ 
erence  7. 
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The  Stress  and  Temperature  Dependence  of  Solubility: 


The  stress  (strain)  dependence  of  solubility  must  be  addressed  if  one  wishes  to  obtain  correct 
concentration  predictions.  In  this  section,  Peteriin's  approach34  will  be  used  to  evaluate  the  solu¬ 
bility  S  in  relation  to  temperature  and  stress.  It  will  be  shown  that  the  predicted  sensitivity  of  S  to 
stress,  is  an  order  of  magnitude  lower  than  the  sensitivity  of  the  diffusion  coefficient  to  stress  and 
that  an  Arrhenius-type  of  relation  is  obeyed  for  temperature.  Once  again,  the  free  volume  concept 
will  prove  to  be  extremely  helpful  in  deriving  the  constitutive  equations. 

Peterlin34  used  the  Hidelbrand  treatment  of  the  thermodynamics  of  polymer-solvent  mixtures, 
in  which  the  entropy  of  mixing  appears  as  a  function  of  the  fractional  free  volumes  of  the  polymer 
(f2)  and  of  the  penetrant  (/|).  The  mixture  is  treated  as  a  two-component  Van  der  Waal  liquid  whose 
entropy  is  proportional  to  the  logarithm  of  the  volume  available  for  kinetic  movement,  which  is 
assumed  to  be  equal  to  the  free  volume  of  the  mixture33.  Assuming  that  the  resultant  free  volume 
fraction  of  the  two-component  system  can  be  obtained  by  a  linear  rule  of  mixture, 

/=  v,/,  +v2/2  (41) 


the  following  expression  for  the  chemical  potential  of  the  penetrant  was  obtained: 

HX=£  +  RTlxA  +  y\  +  Ln{  1  -  v2£)]  (42) 


where: 


v*= _ ki _ 

fl  +  (fl  ~fl)v2 


(43) 


Expression  (42)  is  another  form  of  the  Flory- Huggins  equation  shown  in  the  previous  section,  (the 
same  notation  is  used).  Note  that  /Jt  the  free  volume  fraction  of  the  polymer  is  equal  to  the^  pa- 
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rameter  introduced  earlier  and  that  the  free  volume  fraction  f  of  a  low  molecular  penetrant  is  gen¬ 
erally  larger  than  f2.  The  chemical  potential  of  the  sorbent  in  the  vapor  phase  is  given  by: 

VX=fi  +  RTLn-Z-  (44) 

F  T 

where  P/PT  is  the  activity  of  the  vapor  in  terms  of  partial  pressure.  At  equilibrium,  the  chemical 
potential  of  the  sorbent  in  the  vapor  phase  must  be  equal  to  the  cnemicai  potential  of  the  sorbent 
dispersed  in  the  polymer.  Expressing  this  equality  and  assuming  small  sorption  (v2  close  to  1),  gives: 

V'  =~A^r  exp[~(1  +  ;it1)]  (45) 


By  definition,  solubility  in  weight  fraction  can  be  expressed  as: 


P 1  vi 

P  P  '2  v2 


where: 


px  =  density  of  the  sorbent 
p2  =  density  of  the  polymer 


(46) 


It  follows  that  solubility  is  related  to  the  polymer  free  volume  fraction  by: 


g  _  fiP  1 
f\P%P  T 


exp[  —  (1  +  Xi)] 


(47) 


The  fractional  free  volume  appears  as  an  entropic  quantity  measuring  the  probability  for  the  cre¬ 
ation  of  a  sorption  site.  Thus,  by  changing  the  free  volume  of  the  polymer  component,  stress  acts 
predominantly  on  the  entropy  of  the  system  (see  expressions  (42)  and  (43)).  According  to  our  no¬ 
tation,  the  free  volume  fraction  of  the  polymer  under  strain  is  given  by: 

/2=/2°0+-4')  (48) 

h 
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where: 


k 


Substituting  expression  (48)  into  (47)  yields  our  final  result  for  the  strain  dependence  of  solubility: 


fo  Pi 
f\P\P  T 


exp[  -(1  +  xi)](l 


^kk 


fo 


(49) 


Now,  from  the  Flory-Huggins  theory,  we  know  that  the  interaction  parameter  is  related  to  the  heat 
of  mixing  A//,  by: 


X\  - 


A  Hx 
RTv\ 


(50) 


Substituting  (50)  in  (49)  and  noting  that  v2  is  close  to  1,  yields  an  Arrhenius  type  of  temperature 
dependence: 

A//.  _  Etfr  _ 

S  —  exp  (  )  [  1  +  -jr-  ]  (51) 

where:  SA  =  constant. 

A//,  =  activation  energy  for  solubility  (enthalpy  of  mixing). 

Note  that  in  the  case  of  a  vapor  in  contact  with  the  polymer,  the  activation  energy  for  solubility 
also  contains  the  heat  of  condensation  of  the  vapor.  This  is  due  to  the  fact  that  in  the  Hidelbrand 
theory,  the  reference  state  of  the  two  unmixed  components  is  the  liquid  state. 


Thus,  according  to  this  treatment,  the  relative  change  in  solubility  due  to  stress  is  equal  to  the 
mechanical  dilatation  of  the  free  volume  normalized  to  the  initial  free  volume  fraction: 


«5S  =_4 

So  f0 


(52) 


where:  5,  =  solubility  under  zero  strain  for  a  fixed  temperature. 
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An  alternate  and  more  direct  treatment  of  the  strain  dependence  of  solubility  is  to  imme¬ 
diately  assume  that  the  solubility  is  proportional  to  the  total  fractional  free  volume  of  the  polymer. 
This  model  is  based  on  the  assumptions  that  (1)  once  equilibrium  sorption  is  achieved  sorbent 
molecules  fill  all  the  available  free  volume,  including  the  strain-induced  free  volume,  and  that  (2) 
swelling  is  negligible: 

S  =  -y~  (6  +  4)  (53) 

Jo 

which  is  obviously  equivalent  to  expression  (52).  It  will  be  shown  in  Chapter  2  that  the  above  re¬ 
sult  agrees  quite  well  with  data  collected  on  the  carbon  dioxide-low  density  PE  system,  but  that  lack 
of  agreement  is  found  in  the  case  of  water  in  Ultem  1000. 

The  relative  change  of  the  diffusion  coefficient  due  to  a  small  strain  will  be  determined  in 
Chapter  2.  It  will  be  shown  that: 


SD  _  Bd  j 

D  si 


(54) 


Recalling  that  BD  is  close  to  1  and  that  f„  is  quite  small  (0.03),  we  conclude  that  diffusivity 
should  be  more  sensitive  to  strain  than  solubility,  by  one  order  of  magnitude  or  two.  Also,  note  that 
diffusivity  is  much  less  affected  by  strain  in  the  rubbery  range  than  in  the  glassy  range.  This  property 
may  be  attributed  to  the  larger  free  volume  and  to  the  virtual  incompressibility  of  rubbers. 


1.  A  MODEL  FOR  DIFFUSION  IN  POLYMERS.  GOVERNING  EQUATIONS: 


27 


Conclusion 


A  comprehensive  model  for  the  diffusion  of  gases  and  vapors  in  polymers  in  a  uniform  tem¬ 
perature  field  and  in  the  absence  of  a  non-mechanical  driving  force  has  been  discussed.  In  the 
isothermal  case,  all  our  results  can  be  summarized  in  a  system  of  three  governing  equations: 

-f  -  =  V.  {d  (4 C)  [VC  -  cv4]  }  (55) 

C  = - -  (56) 

Z)(4q=  D0  exp /kkjVC  A  (57) 

l*  /0  +  4  +  yC  J 

where:  C  =  normalized  concentration 

S„  =  solubility  in  the  reference  state 

Ps  =  experimentally  or  theoretically  determined  coefficient 

P  =  vapor  pressure 

D„  =  diffusion  coefficient  in  the  reference  state 
4k  -  dilatation  of  the  free  volume 
y  =  coefficient  of  swelling  expansion 
f0  =  reference  fractional  free  volume  of  the  polymer 
BD  =  parameter  inversely  related  to  the  critical  void  size 

It  is  now  firmly  established  that  the  diffusion  boundary  value  problem  of  interest  in  joint 
durability  studies  is  highly  non-linear  and  is  coupled  with  the  mechanical  response  of  the  polymer. 
The  fully  coupled  solution  for  a  number  of  boundary  value  problems  will  be  given  in  Chapter  3, 
using  finite  element  analysis.  Note  that  equations  (55)  to  (57)  are  implicitely  time  dependent.  It 
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follows  that  the  above  closed  form  representations  are  incorrect  when  not  properly  coupled,  in 
some  iterative  scheme,  with  the  equations  describing  the  viscoelastic  response.  (Convolution  no¬ 
tations,  as  in  expression  (23),  could  be  used,  but  would  only  add  confusion,  with  no  real  gain  in 
rigor) 

In  Chapter  2,  experimental  diffusion  data  on  a  number  of  polymer-penetrant  systems  will  be 
used  to  validate  the  diffusivity  and  solubility  models  developed  in  this  section.  Particular  emphasis 
will  be  placed  on  the  effect  of  stress  (strain)  on  the  transport  rate. 
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2.  A  MODEL  FOR  DIFFUSION  IN 
POLYMERS.  EXPERIMENTAL: 


Introduction: 

Governing  equations  for  the  diffusion  of  moisture  in  adhesives  and  more  generally,  for  the 
diffusion  of  small  molecules  in  polymers,  have  been  derived  in  Chapter  1.  A  theory  for  the  diffusion 
coefficient,  based  on  free  volume,  and  a  theory  for  solubility,  based  on  the  thermodynamics  of 
mixtures,  have  been  proposed.  In  Chapter  2,  the  diffusivity  and  solubility  models  are  validated, 
using  selected  experimental  data  from  the  literature  as  well  as  data  collected  by  the  author.  Since 
the  primary  goal  of  this  section  is  to  verify  transport  theories,  the  polymers  studied  need  not  be 
good  adhesives.  For  the  same  reason,  other  penetrants  than  moisture  can  be  used.  The  polymer- 
sorbent  systems  studied  are:  ethylbenzene-polystyrene,  toluene-polystyrene,  toluene- 
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polyvinylacetate,  PET-oxygen,  PET-carbon  dioxide,  PET-water,  Ultem  polyimide-water,  and  low 
density  PE-carbon  dioxide.  Only  the  polyimidc  was  selected  specifically  for  this  study. 

The  new  contributions  are  found  mainly  in  the  area  of  stress  effects.  Two  general  methods 
are  available  for  measuring  transport  properties  in  mechanically  loaded  polymers:  permeation  and 
sorption.  Data  obtained  with  both  techniques  will  be  presented.  Most  of  the  literature  data  fea¬ 
turing  stress  effects  seem  to  have  been  obtained  by  permeation  techniques.  Permeation  techniques 
however,  are  limited  to  tensile  stress  states  because  the  permeation  membranes  cannot  be  loaded 
in  compression  or  in  shear.  (The  small  membrane  thickness  results  in  buckling  under  those  condi¬ 
tions).  For  this  reason,  a  new  modified  sorption  technique  was  designed,  which  permits  the  study 
of  compressive  and  shear  stress  states. 


The  Temperature  and  Concentration  Dependence  of 
Diffusivity  in  Polystyrene  and  in  Polyvinylacetate  in  the 
Rubbery  Range : 


An  expression  for  the  temperature,  concentration  and  strain  dependence  of  diffusivity  above 
Tt  was  derived  in  Chapter  I .  In  the  absence  of  any  mechanical  volumetric  strain,  as  expected  in  the 
rubbery  state,  we  found  that  one  should  have: 


D, 


D  =  -=r-  T  exp< 


Bl 


N 


_  «(T—  T0)  4-  yC 

fo  f0  +  a{T-T0)  +  yCN 


(26’) 
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Expression  (26')  was  fitted  with  diffusion  coefficients  from  experimental  data  obtained  by  Vrentas 
et  a/.14'27  for  the  temperature  and  concentration  dependence  in  the  rubbery  range.  The  data  were 
extracted  from  step-change  sorption  experiments  following  a  complex  procedure  discussed  in  Ref¬ 
erence  27.  The  solute-polymer  sy'‘?ms  investigated  were-  ethylbenzene-polystyrene,  toluene- 
polystyrene,  and  toluene-polyvinylacetate.  Each  set  of  data  provides  diffusion  coefficients  for  a 
number  of  concentrations  and  temperatures. 

Expression  (26')  was  fitted  to  each  of  the  three  data  sets  using  a  nonlinear  least  squares  rou¬ 
tine.  D  was  expressed  in  m1/ sec  in  order  to  compare  the  fitting  parameters  with  the  known  WLF 
constants.  The  resulting  parametric  curves  with  the  corresponding  experimental  data  for  the 
ethylbenzene-polystyrene,  the  toluene-polystyrene,  and  the  toluene-polyvinylacetate  systems,  are 
shown  in  Figures  1,2  and  3  respectively.  By  itself,  the  excellent  correlation  obtained  only  shows 
that  the  mathematical  form  of  Expression  (26')  is  correct.  In  order  to  demonstrate  the  predictive 
capability  of  the  model,  the  fitting  parameters  were  compared  with  the  expected  constants  as  given 
by  the  free  volume  theory  of  Williams,  Landel  and  Ferry.  The  results  are  summarized  in  table  1  and 
2.  Within  a  small  error,  the  obtained  values  for  BDlft,  Dt  and  Tt  were  correct  and  were  independent 
of  the  nature  of  the  solvent.  In  addition,  the  obtained  Tt  values  for  the  two  polymers  were  quite 
close  to  the  actual  Tt' s. 

Since  the  amount  of  swelling  is  so  large  in  the  three  polymer-solvent  systems  investigated,  our 
model  was  in  fact  pushed  beyond  its  expected  solvent  mass  fraction  range.  In  principle,  the 
Vrentas-Duda  model  would  be  more  appropriate  when  such  high  solvent  concentrations  are  found 
in  the  polymer.  In  practice,  however,  the  number  of  fitting  parameters  is  so  large  in  the  Vrentas- 
Duda  model  (14  parameters,  versus  7  in  Expression  (26'),  that  from  a  purely  mathematical  point 
of  view,  its  curve  fitting  capability  is  not  improved  significantly  over  that  of  Expression  (26').  We 
conclude  that  Expression  (26')  is  an  excellent  phenomenological  representation  of  the  temperature 
and  concentration  dependence  of  the  diffusion  coefficient  in  the  rubbery  range.  Since  moisture 
concentrations  in  adhesives  are  relatively  low  compared  to  the  concentrations  encountered  in  the 
three  systems  above,  we  a-fortiori  expect  excellent  results  in  a  study  of  moisture  diffusion  in  adhe¬ 
sives,  in  a  comparable  temperature  range. 
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The  Stress  Dependence  of  Diffusivity  in  Amorphous  PET 
and  in  Semicrystalline  PET: 


The  predictive  capability  of  the  Free  Volume  Theory  for  stress  dependence  in  a  glassy 
polymer  was  verified  using  experimental  data  collected  by  the  authors  of  Reference  28.  The  dif¬ 
fusion  of  Oj,  CO  2  and  H2  0  through  strained  Poly  Ethylene  Terephtalate  films  and  through  pressurized 
PET  containers  was  studied.  The  measured  parameter  was  permeability  P,  which  is  related  to  the 
diffusion  coefficient  (£>)  and  solubility  (5)  by: 

P=DS  (56) 

Experimental  Procedure: 

A  standard  permeation  apparatus  manufactured  by  Modem  Controls  was  used.  A  schematic 
of  the  permeation  cell  is  given  in  Figure  4.  The  sorbent  and  carrier  gas  were  both  maintained  at 
one  atmosphere.  Thus,  the  driving  force  for  permeation  was  a  partial  pressure  differential;  not  an 
absolute  pressure  differential.  The  films  were  loaded  externally  and  sandwiched  between  the  two 
elements  of  the  permeation  cell.  Due  to  the  absence  of  internal  pressure,  a  barrier  of  silicon  grease 
at  the  film-cell  interface  was  sufficient  to  prevent  any  leakage.  The  PET  films  were  mounted  in  the 
permeation  cell  under  zero  stress  and  oxygen  was  allowed  to  reach  a  steady  state  flux  through  the 
membrane.  The  films  were  then  progressively  stressed  to  higher  levels;  steady  state  permeation  was 
re-established  at  each  level.  Two  types  of  PET  films  were  studied:  extruded  and  biaxially  oriented. 
The  stretched  films  were  produced  at  100  C  on  a  long  film  stretcher.  Orientation  was  4X  by  2.8X. 
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In  another  experiment,  two-liter  PET  bottles  (oriented  4X  by  2.8X  upon  blowing)  were 
pressurized  with  carbon  dioxide.  The  bottles  were  capped  with  a  gage/valve  assembly  and  pressur¬ 
ized.  The  sorbents  studied  were  carbon  dioxide  and  water,  ie:  not  necessarily  the  pressurizing  gas. 
A  permeability  value  was  calculated  from  the  sorbent  flux,  the  internal  pressure  and  the  average 
bottle  thickness.  The  bottles  were  assumed  to  be  cylindrical  and  stresses  were  calculated  using  thin 
wall  vessel  theory.  All  the  data  reported  in  this  section  were  collected  at  room  temperature,  by  the 
authors  of  Reference  28. 


Results  and  Discussion: 


Evaluation  of  Parameters  B,  ft  and  y  in  the  case  of  amorphous  PET: 

Since  our  constitutive  equations  have  been  derived  in  terms  of  the  volumetric  expansion  of 
the  free  volume,  we  need  a  stress-strain  relationship  in  order  to  use  the  available  permeability  data 
reported  as  a  function  of  the  stress  state.  For  a  fully  amorphous  material,  i/kk  is  related  to  the 
hydrostatic  stress  okk  by. 


Pf  ,c7 , 

=  —  °kk  (57) 

where:  Pf  =  compressibility  of  the  free  volume.  (As  explained  earlier,  the  occupied  volume  in  the 
sense  of  Turnbull  is  compressible).  From  References  13  and  30,  we  find  that:  0.033  and  we 

learn  that  for  amorphous  PET : 


kk 


fif 

—  ~  0.36  cm2lkg 

From  Reference  30:  af  ^  1.6510-4  °C~l.  The  above  value  for  af  was  given  for  a  semi-crystalline 
PET.  Assuming  af  is  unchanged  for  amorphous  PET  we  have: 
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pf=  0.6  10-5  cm *jkg  (4.17  10“7  inJ/lb) 


B  can  be  evaluated  from  af  and  the  WLF  parameters  C,  and  C2  of  the  shift  factor: 


B  =  2.303  ayC,  C2 


Using  the  standard  values  for  an  amorphous  polymer  (C,  =  17.44  and  C2  —  51.6)  wc  find: 


B  =  0.34 


Since  the  penetrant  molecules  of  interest  are  small,  we  will  further  make  the  assumption  BDaB. 
The  value  of  y  does  not  need  to  be  calculated  with  accuracy  in  this  analysis.  It  will  be  sufficient  to 
know  that  it  can  be  neglected  for  low  activity  penetrants.  (For  the  penetrants  considered  in  this 
section,  y  is  probably  of  the  order  of  l0-<). 


The  stress  dependence  of  D  and  P  in  the  glassy  state: 
Equation  (27)  can  be  rewritten  as  follows: 


r>_n  B^_  yC  B°  *kk 

=  0CXP  4  /,  +  /tt  +  rC"J  P_  4  /,  +  4  +  rC* 


Since  the  penetrants  of  interest  (02i  C03  and  H2 0)  have  a  low  activity  in  PET,  swelling  is  minimal 


and  the  following  approximation  is  valid: 


yC"  fg 


Furthermore,  for  low  strain  levels,  we  also  have: 


4*  «  fg 


Therefore  expression  (60)  may  be  rewritten  as: 
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At  low  strain  levels,  then: 


D  ^  D0 


(64) 


Combining  expressions  (57)  and  (64)  gives  the  relative  change  in  diffusivity  due  to  an  externally 
applied  load: 


AD  BD  Pf_ 

D  ~  A  3  °kk 

Jg 


(65) 


If  the  Hidelbrand-Peterlin  theory  is  obeyed,  S  is  a  much  weaker  function  of  akk  than  D  (See  the 
section  on  solubility  in  Chapter  1)  and  we  can  rewrite  equation  (65)  in  terms  of  a  relative  change 
in  permeability: 


AP  BD  £f_ 

PAT,  °kk 

h 


(66) 


Expressions  (65)  and  (66)  are  valid  for  small  stresses  only.  They  predict  the  initial  variation  of  P 
and  D  due  to  an  externally  applied  load. 


Correction  For  Crystallinity: 

As  a  first  approximation  the  following  correction  can  be  used  for  semi-crystalline 
polymers25: 


D  =  Da  (l-Xc) 


(67) 


or: 
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P  =  Pa  (1  —  2XC  +  X*) 


(68) 


where:  P„  and  Da  are  the  permeability  and  diffusivity  for  the  amorphous  phase  and  Xc  is  the  degree 
of  crystallinity.  The  above  correction  accounts  for  the  fact  that  the  amorphous  phase  dominates 
transport  properties.  Note  that  this  correction  does  not  appear  in  the  relative  changes 
or  A  more  important  correction  must  be  made  however:  In  a  semi-crystalline 

polymer,  the  stress  experienced  by  the  amorphous  phase  is  smaller  than  that  experienced  by  the 
material  as  a  whole.  This  is  due  to  the  fact  that  (1)  the  amorphous  regions  have  a  higher 
compressibility  than  the  crystalline  regions  and  (2)  displacement  continuity  at  the  interface  between 
the  amorphous  phase  and  crystalline  phase  must  be  satisfied23  .  Therefore  (66)  must  be  modified 
slightly  by  introducing  the  ratio  <j>  of  the  effective  stress  (experienced  by  the  amorphous  region)  to 
the  externally  applied  stress: 


A  P 
P 


h 

T  °kk 


(69) 


The  morphology  of  a  semi-crystalline  polymer  depends  on  the  circumstances  of  crystallization  and 
subsequent  thermal  and  mechanical  treatment'3.  As  a  result,  <j>  is  expected  to  be  a  strong  function 
of  the  specimen  thermal  and  mechanical  history,  making  predictions  of  <j>  nearly  impossible.  We 
will  now  check  the  agreement  of  expression  (69)  with  experimental  data  borrowed  from  Reference 
28.  These  data  show  the  dependence  on  stress  of  the  permeability  of  PET  to  02,  CO2  and  H20. 

The  Case  of  Extruded  PET: 

The  extruded  PET  studied  was  slightly  crystalline:  4%.  This  will  allow  us  to  verify  ex¬ 
pression  (69)  for  quasi-amorphous  PET  i.e.:  <f>  as  1. 

The  experimental  data  are  summarized  in  Figure  5.  Two  cases  were  studied:  a  balanced  biaxial 
stress  state  and  a  uniaxial  stress  state.  Permeability  was  found  to  be  independent  of  the  direction 
of  the  applied  stress,  indicating  that  the  extruded  films  were  isotropic. 
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akk  =  2 a  and  akk  =  a  for  a  biaxial  and  uniaxial  stress  state  respectively.  For  a  balanced  biaxial 
stress  state,  the  initial  permeability  changes  for  amorphous  PET  can  then  be  predicted  by: 


A P  -  B°  h 

—~2  -j  ~  " 

h 


(70) 


and  similarly,  for  a  uniaxial  stress  state: 


AP  _  B°  h_ 

PAT, 

h 


(71) 


Expression  (70)  gives: 
AP 


A  P 
Pa 


—  12.6  Gpa~x  (8.7  1 0_5psi_I)  and  expression  (71)  gives: 


Pa 


=  6.3  Gpcrx  (4.4  10  5psi-').  Figure  5  shows  that  the  agreement  with  experimental  data  is 


remarkable.  The  doubling  of  the  initial  slope  from  uniaxial  to  biaxial  stress  states  is  a  good  check 
of  the  proposed  free  volume-based  model.  Examination  of  (60)  suggests  that  the  permeability 
should  increase  more  slowly  at  higher  stresses.  This  is  not  borne  out  of  Figure  5  and  is  perhaps 
indicative  of  damage  formation  or  morphology  changes. 


The  Case  of  Oriented  PET: 

As  indicated  earlier,  permeability  tests  were  also  conducted  on  biaxially  oriented  PET.  The 
main  differences  in  transport  properties  with  extruded  PET  are  related  to  the  higher  degree  of 
crystallinity,  resulting  in  a  lower  value  of  <)>  (A  degree  of  ciystallinity  of  28%  was  reported  for  the 
oriented  PET  studied  in  this  section).  <j>  was  evaluated  on  semicrystalline  PET  by  the  authors  of 
Reference  30.  A  value  of  .694  was  reported  but  unfortunately  cannot  be  reused  for  specimens 
having  undergone  a  different  thermal  and  mechanical  history.  Since  the  available  data  points  for 
oriented  PET  were  obtained  from  pressurized  cylinders,  we  have:  akk  =  1.5  ah  ,  where:  ah  is  the 
hoop  stress. 

It  is  found  that,  in  order  to  obtain  a  good  fit  with  experiment,  <f>  must  be  close  to  0.46.  Re¬ 
sults  pertaining  to  oriented  PET  are  summarized  in  Figure  6.  Our  value  of  <f>  suggests  a  different 
crystalline  morphology  than  in  the  sample  studied  by  the  authors  of  Reference  30. 
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Summary: 


Our  findings  on  the  stress  dependence  of  D  and  P  can  be  summarized  as  follows: 

1 .  The  free  volume  approach  can  predict  the  initial  stress-induced  variation  of  the  permeability 
of  an  amorphous,  glassy  polymer.  The  slope  prediction  for  amorphous  PET  is  excellent. 

2.  The  doubling  of  the  slope  from  a  uniaxial  stress  state  to  a  biaxial  stress  state  constitutes  a  good 
check  of  the  free  volume  concept. 

3.  A  correction  can  be  made  to  account  for  the  degree  of  crystallinity.  This  is  due  to  the  fact  that 
the  volumetric  stress  level  experienced  by  the  amorphous  phase  (dominating  transport  prop¬ 
erties)  is  smaller  than  that  experienced  by  the  material  as  a  whole.  Unfortunately  <f>  is  difficult 
to  predict,  due  to  a  complex  dependence  on  the  sample  thermal  and  mechanical  history. 

4.  Expression  (65)  reveals  an  interesting  property  of  polymers: 

Materials  with  a  low  diffusivity  (i.e.  a  low  ft)  will  be  more  sensitive  to  stress  because  is 
inversely  proportional  to  ft 2 
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The  Strain  Dependence  of  Diffusivity  in  a  Thermoplastic 
Polyimide: 


Experimental  Procedure: 


The  polymer  selected  for  this  study  was  an  amorphous  thermoplastic  Polyetherimide  with  the 
commercial  name  of  Ultem31  1000  {Tg  =  217  C  ).  The  measurement  method  used  was  sorption. 
The  classical  sorption  experiment  features  a  plane  sorbate  sheet  of  known  thickness,  suspended  in 
an  atmosphere  maintained  at  constant  temperature  and  at  constant  pressure.  Mass  uptake  is  plotted 
as  a  function  of  square  root  of  time,  and  D  is  extracted  from  the  initial  slope.  In  our  study,  sample 
geometries  were  designed,  to  permit  deformation  of  the  polymer  in  tension,  compression,  or  shear; 
while  monitoring  small  weight  changes  simultaneously.  This  was  accomplished  by  using  a  plas¬ 
tically  deformable  substrate  as  a  light-weight  loading  fixture.  The  arrangement  for  the  sorption  cell 
was  an  adaptation  of  the  Me- Bain  apparatus6  ,  In  order  to  avoid  periodic  removal  of  the  sample 
from  the  sorption  cell,  the  sample  was  suspended  directly  from  a  digital  balance  by  a  Nylon  wire. 
Relative  humidity  was  controlled  by  circulating  moist  air  originating  from  a  vapor  exchanger  con¬ 
taining  a  saturated  salt  solution.  For  optimal  temperature  control,  the  sorption  cell  and  vapor 
exchanger  were  both  kept  immmersed  in  a  thermostated  water  bath.  The  digital  balance  was  a 
Sartorius  1 800  with  a  precision  of  0. 1  mg  and  with  a  serial  I/O  port  which  was  used  for  real  time 
data  acquisition.  A  schematic  of  the  set-up  is  given  in  Figure  7.  The  tests  were  conducted  at  60  C 
and  at  a  relative  humidity  of  87%,  generated  by  a  sodium  chloride  solution.  (Early  tests  of  the 
apparatus  showed  that  a  100%  RH  produces  condensation  on  the  specimen  during  small  temper¬ 
ature  fluctuations) 
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Specimen  Preparation: 


Ultem  1000  was  selected  for  four  main  reasons:  (1)  its  excellent  adhesion  to  aluminum  made 
it  attractive  for  the  kind  of  specimen  geometry  used,  (2)  its  thermoplastic  nature  facilitates  the  in¬ 
terpretation  of  the  results  in  light  of  the  Free  Volume  Theory,  (3)  the  relatively  high  diffusion  co¬ 
efficient  of  water  in  Ultem  permits  rapid  data  collection  (Z)  in  Ultem  is  roughly  10  times  higher  than 
in  an  epoxy),  and  (4)  stress  relaxation  is  not  a  concern  within  the  time  frame  of  the  tests,  due  to  the 
large  temperature  gap  between  the  test  temperature  (60  C)  and  Ts  (217  C).  In  order  to  obtain 
tensile  or  compressive  strains,  Ultem  films  were  bonded  to  the  surface  of  a  rectangular  aluminum 
beams  (A1  6061),  and  the  beams  were  then  plastically  deformed  to  the  desired  residual  curvature 
with  a  four  point  bending  fixture  (The  beam  substrate  concept  was  originally  proposed  by 
Putter32).  Beam  curvature  was  uniform  between  the  two  central  loading  points  and  could  yield 
residual  longitudinal  strain  values  on  the  outer  skin,  as  high  as  3.5%,  in  either  a  tensile  or 
compressive  mode.  At  this  point,  it  is  important  to  note  that  although  the  largest  residual  strains 
obtained  in  this  study  were  well  below  the  7-8%  tensile  elongation  yield  limit  of  Ultem31,  our  cal¬ 
culations  indicate  that  the  yield  strain  of  Ultem  1000  could  have  been  exceeded  during  the  defor¬ 
mation  operation  in  some  cases.  This  is  simply  due  to  the  fact  that  the  recoverable  elastic 
component  of  the  beam  deformation  typically  is  four  times  larger  that  the  plastic  component.  Thus 
the  lowest  residual  strain  which  may  have  been  associated  with  prior  plastic  yield  was  estimated  to 
be  around  1.5%.  Following  the  bending  procedure,  the  extremities  of  the  beam  with  non-uniform 
curvatures  were  cut  off.  The  specimen  dimensions  were  the  product  of  a  compromise  between  the 
following  requirements:  (1)  an  upper  limit  on  specimen  weight  was  imposed  by  the  capacity  of  the 
balance  used  (lOOg  in  our  case),  (2)  the  largest  possible  beam  thickness  was  desirable  in  order  to 
minimize  beam  deflections  and  lower  the  experimental  error  on  strain,  (3)  the  film  thickness  had 
to  be  small  compared  to  the  beam  thickness,  in  order  to  produce  a  uniform  through-the-thickness 
strain  distribution  and  in  order  to  limit  test  duration,  and  (4)  the  weight  of  the  film  had  to  be  large 
enough  so  that  the  mass  gain  would  not  too  small  compared  to  the  weight  range  of  the  balance. 
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The  beam  thickness  and  film  thickness  were  12.7  mm  and  0.38  mm  respectively.  The  specimen 
width  was  25.4  mm.  The  aluminum  substrate  was  vapor  degreased,  grit  blasted,  FPL  et.hed,  and 
primed  with  a  solution  of  L'ltem  in  methylene  chloride  (5g  of  L’ltem  per  100  ml).  The  him  was 
bonded  by  compression  molding  at  250  C  and  14  Mpa  (2000  psi)  in  a  hot  press.  This  bonding 
procedure  produced  a  polvmer-to-metal  bond  of  sufficient  quality  to  withstand  large  deformations 
as  well  as  environmental  attack  for  the  duration  of  the  test.  Note  that  for  small  strains,  an  elastic 
deformation  was  produced  in  the  film.  This  can  be  explained  by  the  fact  that  stress  relaxation  was 
minimal  in  the  time  frame  of  the  sorption  tests  (one  day  for  a  thickness  of  .38mm).  Finally,  all  the 
specimens  were  conditioned  in  a  vacuum  oven  at  60  C,  for  three  days,  prior  to  plastic  deformation 
and  testing. 

In  order  to  obtain  a  shear  strain  state,  samples  were  made  consisting  of  a  solid  aluminum 
cylinder  coated  with  a  thin  polyimide  film  of  uniform  thickness.  This  was  accomplished  by  molding 
a  square  block  of  of  l'ltem  around  the  cylinder,  which  was  then  machined  down  with  a  lathe  w'th 
a  tolerance  of  +  /-2 5  microns.  The  cylinders  were  then  plastically  deformed  in  a  torsional  jig, 
thereby  producing  a  quasi-uniform  shear  strain  in  the  film  bonded  to  their  outer  skin.  The  design 
requirements,  surface  treatments  and  bonding  technique  were  the  same  as  with  the  beam  specimens. 
The  cylinder  diameter  and  film  thickness  were  12.6  mm  and  0.5  mm  respectively.  Shear  strains  as 
large  as  4%  could  be  obtained.  All  the  torsional  specimens  were  conditioned  in  a  vacuum  oven  at 
60  C  for  five  days,  prior  to  plastic  deformation  and  testing. 


C 


Results  and  Discussion: 


A  typical  sorprion  curve  obtained  from  our  tests  is  shown  in  figure  8.  The  initial  mass  gain 
was  not  linear  when  plotted  against  time,  but  became  linear  when  plotted  against  sqm 'e  root  of 
time.  Figure  9  illustrates  this  result  by  showing  a  plot  of  the  fractional  weight  increase  (multiplied 
by  a  factor  of  two),  against  the  square  root  of  time  normalized  to  the  film  thickness.  The  straight 
line  up  to  70%  of  the  saturation  level  is  indicative  of  Fickian  behavior  with  a  constant  ddlusion 
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coefficient.  The  factor  of  two  on  the  ordinate  axis  originates  from  the  fact  that  only  one  side  of  the 
film  is  exposed  to  water  vapor.  Since  initial  moisture  penetrations  proceeds  as  in  a  semi-infinite 
medium,  the  initial  rate  of  mass  uptake  was  half  of  what  it  would  have  been  in  a  free  film.  The 
correction  factor  simply  enabled  us  to  extract  the  diffusion  coefficient  by  measuring  the  initial  slope 
according  to  the  classical  procedure. 

Further  examination  of  Figure  8  reveals  that,  as  the  maximum  mass  is  approached,  a  depar¬ 
ture  from  ideal  Fickian  behavior  occurs,  in  the  form  of  a  sudden  excess  sorption.  The  excess  mass 
was  consistently  observed  on  most  specimens.  Since  this  effect  is  also  observed  on  neat  films,  it 
must  be  related  to  some  kind  of  bulk  viscoelastic  response  of  Ultem.  The  exact  nature  of  the 
mechanism  involved  is  still  unknown  at  this  stage.  It  is  fortunate  that  this  relaxation  phenomenon 
did  not  come  into  play  at  an  earlier  stage  of  sorption,  when  it  would  have  most  likely  given  rise  to 
a  non- Fickian  kinetics.  The  dashed  line  in  Figure  8  represents  the  extrapolated  Fickian  sorption. 
Since  the  anomalous  behavior  is  transient  and  short,  the  extrapolated  asymptote,  and  not  the  actual 
maximum,  was  used  for  the  determination  of  the  solubility  of  water  vapor  in  Ultem. 

Figure  10  shows  the  diffusion  coefficient  of  water  in  Ultem,  normalized  to  its  value  at  zero 
strain,  versus  the  longitudinal  component  of  the  strain  (on  the  outer  skin  of  a  plastically  bent 
beam).  Positive,  as  well  as  negative  strain  values  are  plotted.  For  small  strains,  a  linear  response, 
consistent  with  our  theoretical  predictions  for  an  amorphous  polymer  in  the  elastic  range  can  be 
seen.  Note  that  the  compressive  strain  state  lowers  diffusivity,  whereas  the  tensile  stress  state  in¬ 
creases  diffusivity,  as  expected.  Further,  whenever  the  longitudinal  strain  exceeds  1  %  in  absolute 
value,  a  slope  reversal  occurs.  A  very  large  decrease  in  D  at  high  tensile  elongations,  similar  to  what 
we  encountered  in  our  study,  has  been  reported  in  polyethylene  drawn  at  60  C  and  tested  at  25  C. 
The  sharp  drop  was  attributed  to  orientation-induced  crystallization25  .  Because  of  the  lack  of 
symmetry  in  Ultem's  molecular  chains  however,  this  mechanism  is  unlikely  to  have  occured  here. 
As  for  the  increase  of  D  at  large  compressive  strains,  damage  growth  appears  to  be  a  plausible  ex¬ 
planation,  since  defects  are  known  to  enhance  the  transport  rate  in  polymers7.  However,  no  damage 
could  be  detected  visually  in  any  of  the  specimens.  At  this  point,  it  is  instructive  to  recall  that 
plastic  deformation  is  likely  to  have  developed  in  the  specimens  featuring  a  residual  surface  longi- 


2.  A  MODEL  FOR  DIFFUSION  IN  POLYMERS.  EXPERIMENTAL: 


43 


tudinal  strain  as  low  as  1.5%,  which  is  precisely  the  strain  level  at  which  we  see  our  theoretical 


predictions  fail.  The  precise  molecular  events  which  could  lead  to  such  an  abrupt  change  in  dif¬ 


fusion  behavior  under  plastic  deformation  are  not  known  at  this  time. 


Recalling  that  the  compressibility  of  the  free  volume  is  roughly  one  half  of  that  of  the  specific 


volume,  the  initial  slope  in  Figure  10  can  be  related  to  a  slope  with  respect  to  the  dilatational  strain 


of  the  free  volume,  by  the  following  approximation: 


AD  _  AD  2 

D„  Jbir  Do*xx  (1-vi-v.) 


V!  =  lateral  Poisson's  Ratio  (Plastic  deformation  of  the  beam) 


vj  =  Poisson's  ratio  of  the  polymer 


Poisson' s  ratio  v,  was  approximated  by  the  value  0.5  because  the  residual  deformation  in  the 


beam  (after  unloading  the  bending  fixture)  is  dominated  by  the  plastic  component  of  strain,  and 


plastic  flow  is  known  to  occur  with  no  volume  change.  Poisson's  ratio  v2  was  reported33  to  be  0.435. 


Once  multiplied  by  the  correction  factor  introduced  in  expression  (72),  the  experimentally  measured 


slope  for  small  strains  (in  Figure  10),  was  found  to  be  equal  to  1230  (dimensionless).  According  to 

BD 

the  Free  Volume  Theory,  this  value  must  be  equal  to  — -pp.  If  BD  is  assumed  to  be  equal  to  0.4,  a 

h 

equal  to  0.018  is  found.  If  BD  is  assumed  to  be  equal  to  1.0,  a  ft  equal  to  0.029  is  found.  (.4  to 


1  is  the  range  for  B  suggested  in  the  literature13)  Noting  that  these  two  values  are  close  to  the  lower 


and  upper  limits  quoted  for  ft  in  the  literature13,  we  conclude  that  expression  (64)  is  obeyed  at  low 
strain  levels  and  that  the  assumption  that  the  stress  dependence  of  D  below  Tt  is  free  volume  gov¬ 


erned,  is  verified. 


Figure  1 1  shows  the  diffusion  coefficient  of  water  vapor  in  Ultem,  normalized  to  its  value  at 


zero  strain,  versus  shear  strain  (on  the  skin  of  a  plastically  twisted  cylinder).  The  Lagrangian  defi¬ 


nition  of  strain  is  used.  Diffusivity  is  found  to  be  independent  of  the  level  of  shear  strain.  Since  a 


pure  shear  strain  state  does  not  alter  the  specific  volume  of  the  material,  the  above  result  constitutes 


a  crucial  verification  of  the  Free  Volume  Theory. 
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Sumisiary: 


Our  findings  on  the  strain  dependence  of  D  can  be  summarized  as  follows: 

1.  Expression  (64)  can  predict  the  initial  strain-induced  variation  of  diffusivity  in  the  water-Ultem 
system,  both  under  a  tensile  stress  state  and  a  compressive  stress  state. 

2.  The  initial  reduction  of  D  under  a  compressive  stress  state,  constitutes  an  excellent  check  of 
the  free  volume  concept. 

3.  Shear  strain  does  not  affect  D. 


The  Strain  Dependence  of  Solubility  in  Low  Density 
Polyethylene  and  in  a  Thermoplastic  Polyimide: 


The  Case  of  Carbon  Dioxide  in  Low  Density  Polyethylene: 


A  behavior  consistent  with  the  theories  proposed  in  Chapter  1  has  been  reported  in  Reference 
25,  where  the  authors  studied  the  permeation  of  C02  in  low  density  polyethylene  films  under  strain 
(see  Figure  12).  For  small  strains,  expression  (51)  is  able  to  describe  the  behavior  correctly  within 
the  experimental  error.  The  main  result  is  the  confirmation  that  solubility  is  fairly  insensitive  to  low 
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strains.  Practically,  this  means  that  if  a  polymer- penetrant  system  exhibiting  similar  behavior  was 
encountered  in  a  durability  study,  the  solubility  could  easily  be  assumed  to  be  stress  independent. 
For  large  strains,  expression  (51)  also  leads  to  correct  predictions25  (see  Figure  12),  but  it  is  difficult 
to  lend  any  significance  to  this  result  in  view  of  the  fact  that  crystallinity  appears,  leading  to  prop¬ 
erties  which  cannot  be  adequately  described  by  an  unmodified  free  volume  theory.  It  was  also 
shown  in  Chapter  1  that  diffusivity  is  much  more  sensitive  to  strain  than  solubility,  roughly  by  an 
order  of  magnitude;  this  is  also  quite  evident  in  the  experimental  data  presented  in  Figure  12.  We 
conclude  that  the  solubility  models  based  on  the  thermodynamics  of  mixtures  approach  and  the  free 
volume  occupation  approach  seem  to  work  reasonably  well  at  low  strains  for  the  Polyethylene 
carbon  dioxide  system.  The  fact  that  one  finds  agreement  below  Tg,  where  thermodynamic  equi¬ 
librium  is  not  achieved,  is  not  too  surprising  when  one  considers  that  the  derivation  were  based 
primarily  on  geometric  considerations  (either  lattice  theory  or  free  volume  occupancy). 


The  Case  of  Water  Vapor  in  Ultem  Polyimide: 


The  solubility  data  discussed  in  this  section  were  collected  along  with  the  diffusivity  data 
presented  in  Figures  10  and  11.  The  same  specimens  and  the  same  conditions  were  used.  Figure 
13  shows  the  relative  change  of  the  solubility  of  water  in  Ultem  films  subjected  to  varying  levels 
of  longitudinal  strains.  Re-using  expression  (72)  and  replacing  D  by  S  yields  a  slope  with  respect 
to  the  dilatation  of  the  free  volume,  equal  to  1077  (dimensionless),  which  is  roughly  the  slope  we 
found  for  the  diffusion  coefficient.  According  to  expression  (51),  this  slope  should  be  close  to  40. 
The  reason  behind  this  sharp  departure  from  the  theory  could  be  related  to  the  fact  that  the  hy¬ 
drogen  bond  network  of  the  polyimide  is  gradually  broken  up  under  strain,  thereby  releasing  an 
increasing  number  of  hydrogen  bond  sites  for  water.  Unfortunately,  this  explanation  is  not  con¬ 
sistent  with  our  data  for  D  under  a  shear  strain.  Figure  14  shows  the  relative  change  of  the  solu¬ 
bility  of  water  in  U'tem  film*  subjected  to  varying  levels  of  shear  strains.  As  in  the  case  of 
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diffusivity,  solubility  is  found  to  be  independent  of  the  level  of  shear  strain.  This  result  was  ex¬ 
pected,  since  a  pure  shear  strain  state  does  not  change  the  specific  volume  of  the  material;  but  it  is 
not  consistent  with  the  hydrogen  bond  argument  developed  above  because  hydrogen  bonds  would 
also  be  broken  up  under  shear. 

Summary: 

1.  Expression  (51)  can  predict  the  initial  variation  of  solubility  under  a  uniaxial  strain,  in  the 
carbon  dioxide-low  density  PE  system,  but  not  in  the  water-Ultem  system. 

2.  A  shear  strain  state  does  not  affect  solubility  in  the  water-Ultem  system,  as  expected. 


Conclusion: 

With  the  exception  of  the  strain  dependence  of  solubility  in  Ultem,  good  agreement  has  been 
found  between  the  theory  and  experiment.  This  leads  us  to  conclude  that  the  models  presented  in 
Chapter  1  hold  very  promising  predictive  capabilities  for  joint  durability  studies.  The  lack  of 
agreement  with  the  solubility  model  for  certain  materials  can  be  easily  remedied,  simply  by  replacing 
expression  (51)  by  an  equation  determined  empirically.  Such  a  procedure  is  perfectly  legitimate  in 
the  framework  of  an  engineering  solution  of  the  durability  problem. 

Thus  far,  the  data  were  mainly  a  description  of  the  short  term  response  of  the  polymer.  Our 
level  of  confidence  is  now  sufficient  to  embark  in  more  complex  simulations,  including  the 
viscoelastic  behavior  of  the  polymer  over  longer  times.  Chapter  3  will  include  a  validation  problem 
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featuring  a  time  dependent  diffusion  coefficient,  as  well  as  numerical  simulations  of  moisture  pen¬ 
etration  in  a  butt  joint. 
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Figure  I.  Validation  for  Temperature  and  Concentration  Dependence:  Ethylbenzene-Polystyrene 
System.  Data  from  Vrcntas  &  Duda  |27|. 
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Figure  2.  Validation  for  Temperature  and  Concentration  Dependence:  Toluene-Polystyrene  System. 
Data  from  Vrcntas  &  Duda  [27|. 
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Figure  5.  Stress  Dependence  of  Permeability:  Oxygen  in  Extruded  PET  [28|. 
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Figure  6.  Effect  of  a  Biaxial  Stress  on  Permeability:  Oxygen,  Carbon  Dioxide  and  Water  in  Oriented 
PET  128]. 
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Figure  8,  Sorption  Curve  for  a  Torsional  Specimen:  The  Ultem  Film  i 
Level  of  .83%.  (60  C). 
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Figure  12.  EfTect  of  Strain  in  the  Carbon  Dioxide-Low  Density  PE  system:  Normalized  D,  S,  and  P 
vs.  Unidirectional  Strain.  (25] 
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3.  A  MODEL  FOR  THE  DIFFUSION  OF 


MOISTURE  IN  ADHESIVE  JOINTS. 
NUMERICAL  SIMULATIONS: 


Introduction: 


Governing  equations  for  the  diffusion  of  moisture  in  adhesive  joints  have  been  derived  in 
Chapter  1  and  validated  experimentally  in  Chapter  2.  In  these  preceding  sections,  the  coupling  of 
diffusion  kinetics  with  the  mechanical  behavior  was  established.  It  was  shown  that  the  highly 
nonlinear  nature  of  the  governing  equations,  as  well  as  their  implicit  time  dependence  made  the  use 
of  an  iterative  numerical  solution  necessary. 

In  Chapter  3,  a  two-dimensional  solution  of  the  fully  coupled  diffusion  problem  is  obtained 
using  the  finite  element  code  NOVA.  NOVA  has  been  under  continuous  development  at  Virginia 
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Tech  by  Dr.  Reddy  and  S.  Roy  for  the  past  two  years.  The  objective  is  to  provide  a  more  accurate 
analysis  of  adhesively  bonded  joints3*-36 .  In  NOVA,  the  mechanical  response  of  the  adhesive  layer 
can  be  modelled  using  Schapery's  or  Knauss'  nonlinear  single  integral  constitutive  laws  for  multi- 
axial  states  of  stress.  Penetrant  permeation  is  modelled  using  the  diffusion  equations  derived  in 
Chapter  1  (An  extensive  validation  study  of  the  various  original  features  of  the  code  is  presented  in 
References  34  to  36).  In  this  section,  one  of  the  most  interesting  capabilities  of  NOVA  will  be 
demonstrated  by  a  numerical  simulation  reproducing  experimental  results  featuring  a 
time-dependent  diffusion  coefficient  in  a  Polystyrene  film  strained  uniaxially.  Further,  the  effect  of 
the  various  forms  of  strain  coupling  on  diffusion  kinetics  will  be  studied  by  simulating  a  butt  joint 
undergoing  moisture  intake  from  the  edges.  The  accompanying  evolution  of  the  stress  and  strain 
fields  in  the  adhesive  layer  will  also  be  presented.  Chapter  3  is  the  fruit  of  a  team  effort  with  S.  Roy. 


Knauss '  Nonlinear  Viscoelastic  Theory: 


Since  the  diffusion-governing  equations  are  coupled  with  the  mechanical  response  via 
volumetric  strain,  constitutive  equations  for  the  viscoelastic  behavior  of  the  adhesive  are  also  needed 
in  order  to  solve  the  fully  coupled  diffusion  problem.  Knauss'  nonlinear  viscoelasticity  theory  is 
the  most  natural  choice  for  the  present  study  because  it  employs  the  same  phenomenological  de¬ 
scription  as  the  diffusion  constitutive  behavior  proposed  in  Chapter  1.  Free  volume  in  the  sense 
of  Turnbull,  is  used  as  a  unifying  parameter  to  describe  changes  in  the  time  scale  of  the  viscoelastic 
response.  Specifically,  the  theory  states  that  mechanical  dilatational  strain,  temperature  and 
sorbent,  simultaneously  act  as  a  time  accelerating  parameter  by  dilating  the  free  volume. 
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In  the  Knauss  approach,  the  basic  single  integral  formulation  of  linear  viscoelasticity  is  used. 
The  stress  and  strain  tensors  are  related  by  the  Stieltjes  convolution  integrals  and  the  material  is 
assumed  to  undergo  small  deformations16.  Consider  an  isotropic  polymer.  Let  J(t)  and  B(i)  be  the 
shear  compliance  and  bulk  compliance  respectively.  The  constitutive  equations  in  the  framework 
of  linear  viscoelasticity  are  given  by: 

1  r‘  dSy 

ey  =  7  J_^  J(‘  ~  f )  -gf  dt  (73) 

1  fr 

•tt-y  J  if-  ^  (74> 

where:  e u  =  dcviatoric  strain 

tkk  -  volumetric  strain  (Mechanical  component) 

S,j  =  deviatoric  stress 
okk  =  volumetric  stress 

and  the  indicial  summation  convention  is  assumed.  The  equation  numbers  are  continued  from 
Chapters  1  and  2. 

Knauss  introduced  the  nonlinearity  by  expressing  that  the  time  scale  of  viscoelastic  response 
as  described  by  7(0  and  B(t)  is  a  strong  function  of  free  volume  dilatation  A/.  The  total  free  volume 
dilatation  A/ is  expressed  as  the  sum  of  a  mechanical  strain  contribution,  a  thermal  expansion 
contribution  and  a  sorbent  expansion  contribution: 

A/=  tfkk  +  a  A T+yC  (75) 

where:  a  =  coefficient  of  thermal  expansion  of  the  free  volume 

y  =  coefficient  of  swelling 

=  volume  dilatation  of  the  free  volume  due  to  external  loads 
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This  approach  has  been  used  in  our  derivation  of  an  expression  for  the  diffusion  coefficient 
in  Chapter  1.  Furthermore,  Knauss  assumed  that:  (1)  free  volume  changes  do  not  change  the  dis¬ 
tribution  function  of  retardation  times  in  J(l)  and  B(t),  implying  that  the  material  is 
thermorhcologically  simple,  and  (2)  the  same  shift  factor  can  be  used  for  both  shear  and  bulk 
properties.  Mathematically,  the  time  is  simply  shifted  by  using  the  differential  expression: 


dt'  = 


dt 

a(4/) 


(76) 


where  t‘  denotes  the  reduced  time  and  the  shift  factor  a(AJ)  is  given  by: 


Ln  {u(4/)}  =  --f 

JO 


A f[T)  +  /kk  +  yC 
fo  +  A f[T)  4-  4*  +  yC 


(77) 


where  A/[T)  is  a  function  describing  thermally-induced  changes  in  the  free  volume  in  the  temper¬ 
ature  domain  of  interest.  For  example,  A/[T)  =  a(T  —  Tt)  above  Tt. 

Thus,  Knauss'  shift  factor  is  an  immediate  generalization  of  the  WLF  shift  factor  (expression 
20).  Note  that  for  =  0  and  yC  =  0,  the  standard  WLF  equation  is  recovered.  Knauss'  nonlinear 
constitutive  behavior  can  now  be  summarized  as  follows: 


1  V  dSn 

M~r)  if dr 


(78) 


f  Bw-r) 


d<Tkk 

dr 


dr 


(79) 


with: 


3.  A  MODEL  FOR  THE  DIFFUSION  OF  MOISTURE  IN  ADHESIVE  JOINTS.  NUMERICAL 
SIMULATIONS: 


67 


c 


O 


««-  L  <so) 

L  (81) 

An  obvious  advantage  in  using  Knauss'  constitutive  model  in  our  study  is  that  the  shift  factor 
subroutine  can  be  shared  by  both  the  diffusion  code  and  the  viscoelasticity  code.  Only  void  size 
factors  B  and  BD  differ  in  principle  in  the  two  boundary-value  problems. 

It  was  noted  earlier  that  NOVA  has  provision  for  either  the  Schapery  or  Knauss  nonlinear 
viscoelastic  theories.  The  basic  form  of  the  constitutive  equations  is  similar  in  both  theories;  they 
only  differ  by  the  form  and  number  of  nonlinearizing  functions.  Fortunately,  it  is  possible  to  reduce 
Knauss'  model  as  a  particular  case  of  Schapery's  model  by  changing  the  expression  for  the  shift 
factor  and  setting  the  remaining  nonlinearizing  functions  to  1.  This  scheme  was  used  in  the  finite 
element  formulation  of  the  viscoelastic  behavior  in  NOVA. 


Numerical  Scheme: 


The  coupling  between  diffusion  and  viscoelasticity  can  be  easily  seen  by  cross-examining  the 
governing  equations  for  the  Diffusion  Boundary-  Value  Problem  (DBVP,  Expressions  (55)  to  (57)) 
and  the  Constitutive  equations  for  the  Viscoelasticity  Boundary-Value  Problem  (VBVP,  Ex¬ 
pressions  (78)  to  (81)).  The  diffusion  equations  are  strain  dependent,  while  the  viscoelastic  behavior 
is  affected  by  a  sorbent  concentration  term  in  the  shift  factor.  The  mechanical  response  is  also  af- 
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fected  by  sorbent  concentration,  in  that  mechanical  strains  can  be  generated  due  to  swelling.  This 
effect  is  incorporated  in  the  governing  equations  by  stating  that  the  total  strain  contains  a  sorbent 
expansion  component. 


The  interplay  between  the  DBVP  and  the  VBVP  is  illustrated  in  Figure  15.  Numerically,  this 
coupling  can  be  implemented  by  solving  the  diffusion  problem  and  viscoelasticity  problem  con¬ 
currently  in  order  to  continually  update  the  coupling  variables,  namely,  the  volumetric  strain  and 
sorbent  concentration.  This  scheme  has  been  encoded  in  NOVA  and  is  summarized  by  the  flow 
chart  shown  in  Figure  16.  Note  that  the  strain  dependence  of  the  solubility  is  not  included  in  this 
study.  The  experimental  data  presented  in  Chapter  2  seem  to  justify  this  approximation  as  long 
as  the  penetrant  does  not  show  much  chemical  affinity  with  the  polymer. 

The  numerical  treatment  presented  herein  introduces  an  implicit  time  dependence  in  the 
diffusivity  and  in  the  solubility,  thereby  allowing  the  simulation  of  more  complex  sorption  behav¬ 
iors  known  as  "anomalous*  or  "case  II".  Case  II  sorption  refers  to  a  case  where  the  rate  of  transport 
is  entirely  governed  by  the  rate  of  viscoelastic  relaxation7.  Anomalous  diffusion  refers  to  an  inter¬ 
mediate  case  between  Fickian  and  Case  II  diffusion. 


Finite  Element  Formulation  of  the  Problem: 


Viscoelasticity  Formulation: 


lire  total  strain  and  stress  components  in  a  material  can  be  written  as  the  sum  of  deviatoric 
and  volumetric  components: 
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cu  =  eu  +  y  6 


—  S,,  +  —  o 


°u =  +  7 


For  a  viscoelastic  material,  the  constitutive  equations  discussed  earlier  can  be  rewritten  in  the  fol 
lowing  form,  if  the  effects  due  to  loading  history  prior  to  t  =  0  are  negligible: 


*kk  =  y  5(°)  VkkU)  +  y  J  M  W  ~  V)  y:  ***(•*■)  • 


?(/  =  y  m  Sip)  +  -}  J'  A/  (*  -  **)  £  Sip)  dr 


where,  <p  and  \{/'  are  given  by  expressions  (80)  and  (81)  respectively. 

Using  results  from  equation  (84)  in  equations  (82)  and  83),  one  obtains  for  two  dimensional  anal 


{e}  =  [C]  •  {c/a} 


where: 


[C]  = 


{e}  -  {su,  e22,  y  12>  e33) 


(  —  +  —  )(  —  -  — )0  (  —  -  —  ) 
1  9  3  A  9  6  ;  1  9  6  ' 

(  —  +  — )  0  (  — -  — ) 
v  9  3  '  v  9  6  ’ 


symmetric 


(*■+') 
v  9  3  ’ 
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.N-  -."V.CV  -.V-/  V  V  •  '  ,-.V  •  -  "  - -  -  -  -  •  v-  O  '  O  O.OO.O  M 


riTWiWiwc 


{da}  =  {dan,  do22,  drl2  da33} 


and  where  the  symbol  (♦)  denotes  the  convolution  operator. 

If  the  transient  compliances  are  now  written  in  the  form  of  a  Prony  series,  then: 


r  I 

A BW)  =  ^  Br  (1  -  exp(  -  iA/rr)) 


Ai(iA)  =  Jr{  1  -  exp(  -  i Hlr)) 


Substitutmg  equations  (86a)  and  (86b)  in  equation  (85a)  results  in  a  matrix  equation  given  by: 


{«}  =  [(V]  {«*}  +  {//} 


where  the  matrix  [N]  contains  the  instantaneous  compliances  at  time  t,  and  the  vector  {//}  contains 
the  components  of  the  hereditary  strains.  Pre-multiplying  equation  (87)  by  [AG-1  and  rearranging, 


{a}  =  [L]  ({c)  -  {//}) 


where: 


[L]  =  [^]- 


The  finite -element  equilibrium  equations  may  be  established  by  invoking  the  principle  of 


virtual  work: 


{<5u}r  [*]7{<t}  dV~{F^  =  {0) 


Using  results  from  equation  (88)  in  (90)  yields, 
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(91) 


(c5u}r^|  IBY  \_L]  [z]  dV  -  J  [5]r[L]  {H)dV-{F)j  = 


Writing  the  strain-displacement  relations  as: 


to  =  [5]  {m} 


and  substituting  in  (91): 


{Su}T  [5]r[L]  [5]  dV{u)  -  J  [5]r[L]  {H)dV- {/}^  = 


or  simply, 


(<5m)7  {/?}  =  {0} 


Noting  that  equation  (93)  contains  a  source  of  nonlinearity  imbedded  in  the  definition  of  the 
shift  factor,  the  Newton-Raphson  iteration  technique  is  employed  to  solve  for  the  displacements. 
For  the  i1”  iteration,  the  incremental  displacements  {Am,}  are  obtained  from: 


{AmJ  =  -  IktT1  {/?,} 


to)  =  {q_i)  +  {Am,-} 


where: 


Ur]=  [  [S]r[L]  [5]  dV 

J  V 
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Moisture  Diffusion  Analysis: 


The  finite  element  formulation  of  Fick's  Law  in  two  dimensions  is  developed  using  the  weak 
variational  form54'36;  the  "weak*  form  of  the  variational  procedure  weakens  the  continuity  require¬ 
ments  on  the  displacements,  by  allowing  discontinuity  in  the  displacement  gradients: 


I  V- 


(d--  +  kdc 

ot  ox  \  ox  ox 


)-i(D 


~  +  KDC^~\  dxdy  =  0(98) 
dy  dy  I 


For  the  plane  strain  case,  the  volumetric  strains  can  be  expressed  as 


dll  .  dV  >— t  A  rp 

EWt= U+ a7~vC_aAr 


Assuming  that  the  moisture  concentration  may  be  approximated  by: 

n 

c(x,y,t)  =  Yj 

j= i 

and  the  test  function  in  expression  (98)  may  be  set  to  be  equal  to  the  interpolation  function: 

Substitution  of  Eq.  (100)  into  the  weak  form  of  Eq.  (98)  (See  Reddy34)  gives 

CMe]  {Q  +  lKe]  {Q  =  {i^} 

where: 


[M*]  =  J  4>(\l/j  dxdy 
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UH-D( 


d\l/t  d\ l>j  dtj/j  d\j/j 
dx  dx  +  dy  '  dy 


d\  dC  d\T  \  #< 
dxdy  y  dx  *  dx  )  dx 


(  d\ 

,  02v 

dC 

dAT  ^ 

dh  11 

l  dxdy 

'  *5  y 

dy  “ 

dy  j 

dy  Jj 

(105) 

dxdy 


4(f" ^  +*c  i^)"' 


The  time  derivative  {C}  is  approximated  by  using  the  0  -  family  of  approximation  for  the  nth  time 
step, 

0{Qn+1  +  (1  -  0){C„}  =  ({Qn+1  -  {C]„)/At„+1  for  0^  1  (107) 

where  0  is  a  weighting  parameter.  From  equations  (99)  and  (106),  we  obtain  for  each  element, 

Ue]  (Qn+1  -  lBel  {C}n  -  V Pe}„  =  {0}  (108) 


where: 


[/4e]  =  [A/e]  4-  0Ar„+]  lKel 
[5s]  =  [A/s]  —  (1  —  0)  At„+1  OT 
[5s]  =  At„+1(0  +(1-0)  {/«}„) 
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Recognizing  that  a  source  of  nonlinearity  in  the  form  of  the  diffusion  coefficient  D  is  imbedded  in 
the  matrix  [AT*3»  the  Newton-Raphson  technique  is  employed  to  solve  for  the  moisture  concen¬ 
trations  {Qn+1  at  each  time  step. 


Validation  Problem:  Viscoelastic  Diffusion  Through  a 
Polystyrene  Film: 


Smith  et  a/.37-38  conducted  a  permeation  experiment  to  study  gas  transport  in  polystyrene  and 
found  that  the  diffusion  coefficient  for  C02,  Ar  and  Xe  decreased  with  time  when  the  polystyrene 
film  was  subjected  to  a  constant  uniaxial  strain.  Figures  17  and  18  reproduce  experimental  results 
from  Reference  37.  Figure  17  illustrates  the  time  dependence  of  the  diffusion  coefficient  for  C02  in 
a  Trycite  film  (biaxially  oriented  polystyrene  film)  at  different  strains  at  50"C.  Figure  18  illustrates 
the  time  dependence  of  the  diffusion  coefficient  for  C03  and  Xe  in  a  Trycite  film  at  the  same  strain 
level  at  50°C. 

The  reason  these  experimental  data  were  found  to  be  convenient  for  a  test  of  NOVA  is  that 
polystyrene  is  adequately  characterized  mechanically.  Reference  13  gives  the  viscoelastic  shear 
compliance  and  the  viscoelastic  bulk  compliance  of  polystyrene  around  the  glass  transition  tem¬ 
perature  (100°C).  Bulk  properties  are  rarely  documented  in  the  literature.  This  is  due  to  a  lack  of 
recognition  for  their  importance,  compounded  by  complexities  in  the  measurement  techniques. 
Yet,  our  study  demonstrates  that  the  diffusion  behavior  as  well  as  the  nonlinear  viscoelastic  be¬ 
havior  are  intimately  related  to  bulk  properties.  In  fact,  the  methodology  developed  here-in  estab¬ 
lishes  that  bulk  characterization  is  an  absolute  necessity  for  improved  diffusion  predictions. 
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Although  the  shear  compliance  and  bulk  compliance  are  given  for  two  different  molecular 
weights  (500,000  and  600,000  respectively),  the  properties  were  assumed  to  be  usable  for  any  high 
molecular  weight  polystyrene.  This  assumption  is  valid  as  long  as  the  molecular  weight  is  large 
compared  to  some  critical  value  (38,000  for  Polystyrene),  corresponding  to  the  onset  of 
entanglement  coupling13. 

The  compliances  around  100°C  given  in  Reference  13  were  curve-fitted  with  a  Prony  Series 
using  a  nonlinear  least-square  fitting  routine.  Since  predictions  must  be  compared  to  data  at  50°C, 
a  time-temperature  shift  had  to  be  performed.  It  is  well  known  that  the  shift  factor  below  the  glass 
transition  temperature  is  governed  by  an  activation  energy  instead  of  an  activation  volume,  as  in 
the  case  of  the  rubbery  state.  The  shift  factor  is  given  by: 


a  =  exp 


AH  (  1 

R  \  Tr 


where:  AH  -  activation  energy 

R  —  gas  constant 
Tr  =  reference  temperature 
T  =  temperature 


(112) 


Values  of  AH  between  30  kcal/mole  and  40  kcal/mole  have  been  suggested  by  Ferry13  and 
Matsuoka17.  A  value  of  35  kcal/mole  was  taken  in  our  study,  which  corresponds  to  a  shift  fretor 
of  1500.  This  means  that  the  retardation  times  at  50°C  arc  1500  larger  than  a  100°C,  reflecting  the 
slower  mechanical  response  of  the  material  at  a  lower  temperature  (time  in  sec).  The  same  shift 
factor  was  used  for  the  bulk  compliance  and  the  shear  compliance17  . 

The  reference  temperature  T„  for  the  diffusion  problem  was  fixed  at  50°C  and  the  reference 
diffusion  coefficient  D0  was  taken  from  the  measured  values  at  50°C  in  Reference  37.  Since  all  the 
data  points  in  reference  37  are  given  at  50°C,  the  temperature  term  in  the  shift  factor  was  set  to  zero 
in  all  the  computations.  The  swelling  coefficient  of  expansion  y  was  set  to  zero  because  the 
sorbents  used  in  the  permeation  studies  have  little  chemical  affinity  with  polystyrene. 


3.  A  MODEL  FOR  THE  DIFFUSION  OF  MOISTURE  IN  ADHESIVE  JOINTS.  NUMERICAL 
SIMULATIONS: 


76 


At  this  point,  it  is  useful  to  recail  that  the  volumetric  strain  4k  contained  in  the  diffusion 
equation  corresponds  to  the  dilatation  of  the  free  volume  and  not  to  the  mechanical  volumetric 
strain  tkk ,  used  in  the  governing  equations  for  viscoelasticity.  Kovac's  model  presented  in  Reference 
13  postulates  that  the  free  volume  dilatational  strain  4k  is  equal  to  the  transient  component  of  the 
mechanical  strain  tkk.  Kovac's  approach  was  used,  and  the  free  volume  dilatation  strain  at  any  time 
was  given  by; 

^kk  ~  Ekk  ~  om  (1 13) 

At  infinite  times,  one  finds  that  4k  is  approximately  equal  to  one  half  of  ckk,  which  is  con¬ 
sistent  with  remarks  made  previously  in  Chapter  2. 

Figures  19,  20  and  21  show  numerical  simulations  produced  by  NOVA,  of  the  problem 
studied  by  Smith  et  al  in  Reference  37.  Figure  19  shows  the  variation  of  the  diffusion  coefficient 
with  time  for  the  three  strain  levels  indicated  in  Figure  17.  From  Figure  19,  it  is  evident  that,  in¬ 
dependently  of  the  strain  level,  the  diffusion  coefficient  reaches  a  peak  at  t  =  1  hour  and  then  slowly 
decays  back  to  the  reference  value  D„.  This  behavior  can  be  attributed  to  an  initial  increase  in  the 
free  volume  due  to  the  application  of  the  uniaxial  strain,  followed  by  a  continuous  recovery  in  free 
volume  at  constant  strain,  as  the  polystyrene  film  undergoes  shear  relaxation  which  gradually  re¬ 
duces  the  hydrostatic  stress  level.  A  larger  applied  strain  produces  larger  initial  dilatation,  resulting 
in  a  higher  peak  in  the  diffusion  coefficient.  Figure  19  also  reveals  that  the  rate  of  change  of  the 
diffusion  coefficient  (to  be  paralleled  with  the  rate  of  change  of  the  free  volume)  becomes  larger  as 
the  applied  strain  level  increases.  This  strain  level  dependence  in  the  response  rate  is  an  excellent 
illustration  of  the  nonlinear  nature  of  the  mechanical  response,  described  here  by  Knauss'  model. 
It  should  be  noted  that  the  initial  increase  in  the  diffusion  coefficient  is  largely  due  to  the  fact  that 
bulk  creep  occurs  on  a  much  shorter  time  scale  than  creep  in  shear  (Our  data  indicate  a  shift  of 
about  five  decades  between  the  bulk  retardation  spectrum  and  the  shear  retardation  spectrum). 
Upon  examining  the  time  scale  of  Figure  19,  it  becomes  clear  that  the  initial  increase  in  the  diffusion 
was  not  experimentally  accessible  to  Smith  et  al  .  For  this  reason,  the  initial  transient  response 
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was  seen  as  instantaneous  by  Smith  et  al.  and  only  the  relaxation  part  of  the  curve  is  shown  in 
Figure  17.  For  the  same  reason,  the  diffusion  coefficient  in  Figure  17  is  normalized  with  respect 
to  its  value  at  a  time  accessible  experimentally,  instead  of  being  normalized  with  respect  to  D0,  as 
in  Figure  18.  Bearing  these  details  in  mind,  it  becomes  apparent  that  the  right  hand  side  of  Figure 
19  correlates  very  nicely,  although  not  exactly,  with  the  experimental  results  of  Figure  17.  The  good 
agreement  between  the  simulation  and  the  experimental  data  is  a  strong  indication  that  NOVA 
holds  real  promise  for  joint  durability  prediction. 

The  influence  of  penetrant  molecule  size  on  the  diffusion  coefficient  for  gases  in  polystyrene 
was  qualitatively  studied  by  varying  the  magnitude  of  void  size  parameter  BD.  Levita  and  Smith37 
explained  that,  due  to  molecular  shape  effects,  the  effective  critical  void  size  of  a  C02  molecule  is 
larger  than  that  of  Xenon.  It  follows,  from  our  discussion  in  Chapter  1,  that  BD  for  C02  should  be 
smaller  than  that  of  Xenon  ( BD  and  the  critical  void  size  Vc  are  inversely  related).  Figure  20  shows 
predictions  obtained  from  NOVA  for  two  values  of  BD  and  a  strain  of  1.8%  in  a  polystyrene  film. 
By  comparing  Figure  20  with  Figure  18,  it  can  be  seen  that  NOVA  correctly  predicts  that  a  smaller 
penetrant  molecular  size  leads  to  a  faster  rate  of  decrease  in  the  diffusion  coefficient.  Since  the  ac¬ 
tual  Bn  values  for  C02  and  Xe  are  unknown,  the  above  comparison  is  only  qualitative  at  this  point. 

The  effect  of  physical  aging  on  the  diffusion  coefficient  was  studied  by  implementing  equation 
(31)  of  Chapter  1  in  NOVA.  The  penetrant  was  C02  and  the  values  of  temperature  strain  anu  t, 
were  set  at  50°C,  1.8%  and  24  hours  respectively.  Figure  21  shows  that  a  faster  rate  of  physical 
aging  denoted  by  a  higher  value  of  parameter  n  ,  causes  the  diffusion  coefficient  to  decay  more 
slowly.  This  behavior  is  expected,  since  a  smaller  free  volume  causes  molecular  relaxation  processes 
to  occur  over  a  larger  period  of  time.  Unfortunately,  no  experimental  data  are  available  for  com¬ 
parison  here. 
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Diffusion  Kinetics  in  a  Butt  Joint 


In  this  section,  a  number  of  numerical  simulations  of  moisture  diffusion  in  the  adhesive  layer 
of  a  butt  joint  are  presented.  The  main  emphasis  is  on  evaluating  the  relative  importance  of  the 
various  coupling  effects  between  diffusion  and  mechanical  strain  (stress).  In  order  to  achieve  this 
goal,  a  parametric  study  was  carried  out  using  realistic  values  for  the  adhesive  properties.  The  ad¬ 
hesive  mechanical  properties  were  those  of  polystyrene  at  50°C.  The  fact  that  polystyrene  is  a  poor 
adhesive  is  of  no  concern  for  this  kind  of  study.  In  order  to  study  the  effect  of  swelling,  a  large 
volume  for  the  coefficient  of  swelling  expansion  (y)  was  needed.  Since  y  for  polystyrene  is  negli¬ 
gible,  a  fictitious  material  combining  the  mechanical  properties  of  polystyrene  and  the  coefficient 
of  swelling  expansion  of  a  moderately  hydrophilic  polymer  (y  =  0.01)  was  used. 

Table  3  summarizes  the  four  cases  investigated  in  the  parametric  study.  (The  numbering  of 
the  cases  run,  should  not  be  confused  with  the  fundamental  modes  of  sorption  known  as  "Case  I' 
and  "Case  II")  The  joint  geometry  and  finite  element  discretization  are  given  in  Figure  22.  Note 
that  the  mechanical  boundary  conditions  used  were  a  uniform  axial  displacement  resulting  in  a 
uniform  strain  of  2%  in  the  adhesive  layer.  The  normalized  moisture  concentration  at  the  free  edge 
of  the  adhesive  layer  is  unity  and  the  initial  concentration  throughout  the  adhesive  layer  is  zero. 

Figure  23  shows  the  moisture  concentration  profiles  within  the  adhesive  layer  at  three  differ¬ 
ent  times  when  there  is  no  coupling  (case  1).  In  this  case,  the  diffusion  coefficient  remains  constant 
with  time,  that  is,  D  =  De.  Case  I  is  the  classical  Fickian  diffusion  in  a  plane  sheet  for  which  a 
closed  form  solution  is  available.  For  short  times,  as  in  our  simulations,  the  sheet  may  be  consid¬ 
ered  as  a  semi-infinite  medium  and  the  following  approximation  is  convenient: 

C(t,y)=  erfc  t)  (114) 

where  erfc  is  the  error  function  complement. 
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Figure  24  shows  the  moisture  concentration  profiles  for  the  case  where  there  is  viscoelastic 
coupling  only,  that  is,  when  both  the  diffusion  coefficient  and  the  shift  factor  arc  dependent  on  the 
transient  component  of  the  dilatational  strain  (Case  2). 

Figure  25  depicts  the  case  where  there  is  full  coupling  in  the  diffusion  coefficient  and  in  the 
shift  factor.  This  means  that  the  diffusion  coefficient  and  the  shift  factor  are  both  a  function  of  the 
dilatational  strain  and  of  the  moisture  concentration  at  any  given  point  in  the  adhesive.  Note  the 
dramatic  increase  in  penetration  rate  relative  to  the  two  previous  cases,  as  well  as  the  convex  shape 
of  the  concentration  profiles.  The  moisture  profiles  for  the  fully  coupled  diffusion  problem  (Case 
4)  are  not  shown  because  they  are  indistinguishable  from  those  plotted  in  Figure  25  (Case  3).  In 
addition  to  all  the  features  of  case  3,  case  4  includes  the  contribution  of  the  strain-induced  driving 
force  VWp.  This  last  case  required  the  use  of  135  time  steps,  compared  to  about  50  time  steps  for 
case  1. 

lastly,  Figure  26  presents  the  results  for  each  of  the  four  previous  cases,  for  comparison  at 
time  t  =  60  minutes.  Note  that  case  3  and  case  4  are  still  indistinguishable  after  60  minutes.  Only 
strain  and  concentration  coupling  in  the  diffusion  coefficient  seem  to  have  a  significant  accelerating 
effect.  Of  course,  the  effect  of  differential  swelling  would  have  appeared  if  a  larger  value  of  the  co¬ 
efficient  of  moisture  expansion  had  been  used.  High  magnitudes  for  y,  promoting  case  II  kinetics, 
are  encountered  when  the  penetrant  is  a  good  solvent  of  the  polymer.  For  structural  adhesives  in 
the  presence  of  moisture  however,  y  is  expected  to  be  close  to  0.01  (An  estimate  of  y  from  data  in 
Reference  39  yields  a  value  of  0.014  for  a  3501-6  epoxy).  The  various  types  of  diffusion  behaviors 
in  polymers  have  been  graphically  summarized  by  Hopfenberg  and  Frish6  in  terms  of  penetrant 
activity  and  temperature  (see  Figure  30).  Examination  of  this  chart  confirms  that  non-Fickian 
driving  forces  are  not  operative  below  Tf  as  long  as  penetrant  activity  is  small,  and  that  a  concen¬ 
tration  dependent  diffusion  coefficient  is  then  sufficient  to  describe  the  transport  behavior. 

The  main  findings  from  this  study  may  be  summarized  as  follows: 

1 .  The  mechanical  and  hygroelastic  coupling  terms  in  the  diffusion  coefficient  must  be  included 
in  a  durability  analysis,  due  to  their  substantial  accelerating  effect  on  penetration  kinetics. 
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2.  The  effect  of  differential  swelling  can  be  counted  as  negligible  for  a  typical  structural  adhesive. 
Polymers  for  which  this  extra  driving  force  becomes  important  are  not  likely  to  be  selected  for 
adhesive  applications,  due  to  their  poor  performance  in  moist  environments. 


Evaluation  of  the  Stress  and  Strain  Fields  in  a  Butt  Joint 
Upon  Moisture  Penetration: 


The  results  presented  in  this  section  were  generated  by  running  case  3  (see  Table  3)  with 
NOVA.  The  bond  line  discretization,  mechanical  boundary  conditions  and  material  properties 
were  the  same  as  earlier  (same  y  and  same  mechanical  properties).  The  stress  and  strain  fields  were 
plotted  at  three  different  elapsed  times:  0  min.,  30  min.  and  60  min.  On  each  plot  presented 
here-in,  free  edge  OB  of  the  butt  joint  is  located  at  the  right-hand  side  and  distance  along  the  bond 
line  is  normalized  with  respect  to  the  semi-width  of  the  joint. 

Figures  27  and  28  illustrate  the  effect  of  moisture-induced  swelling  on  the  mechanical  com¬ 
ponents  of  normal  strain  e„  and  transverse  strain  tyy,  respectively.  In  Figure  27,  the  absorbed 
moisture  induces  swelling  strains  within  the  adhesive,  resulting  in  a  lower  value  of  e„  .  Note  that 
e„  refers  to  the  mechanical  component  of  strain.  In  the  BVP  considered  here,  three  types  of  strains 
are  used:  the  mechanical  strain,  the  hygroscopic  strain  and  the  kinematic  strain.  The  kinematic 
strain  (2  %  here)  always  remains  equal  to  the  sum  of  the  mechanical  and  hygroscopic  strains.  Since 
the  boundary  conditions  used  impose  a  uniform  positive  kinematic  normal  strain  of  two  percent 
initially,  the  lower  values  of  mechanical  strain  e„  in  the  presence  of  moisture,  remain  positive.  Note 
that  in  an  unloaded  butt  joint,  e„  would  have  become  negative,  reflecting  the  constraining  effect 
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of  the  stiff  adherends  surrounding  the  adhesive.  By  comparing  Figure  27  with  Figure  25,  it  appears 
quite  dearly  that  the  region  with  a  lowered  normal  strain  corresponds  to  the  domain  occupied  by 
moisture. 

Figure  28  shows  a  negative  initial  transverse  strain  resulting  from  the  transverse  Poisson's 
contraction  upon  loading.  In  the  absence  of  moisture,  the  algebraic  value  of  tyy  would  tend  to  in¬ 
crease  initially  and  then  decrease,  due  to  the  transient  bulk  response  of  the  polymer.  The  initial  bulk 
creep  can  be  seen  here  in  the  unperturbed  region  at  time  =  30  min.,  causing  the  Poisson's  ratio  to 
decrease.  The  subsequent  bulk  relaxation  can  be  seen  near  the  free  edge  as  the  Poisson's  ratio  in¬ 
creases  again  at  longer  times.  As  in  the  case  of  the  normal  strain,  moisture-induced  swelling  causes 
an  abrupt  change  in  the  distribution  of  s>y,  and  the  transverse  strain  field  in  the  dry  region  seems 
to  be  unchanged. 

Figure  29  depicts  the  effect  of  moisture-induced  swelling  on  the  distribution  of  normal  stress 
cr„.  The  comments  made  for  e„  in  Figure  27  can  be  repeated  here.  The  only  difference  in  behavior 
is  found  in  the  dry  region  ahead  of  the  moisture  front,  where  a  stress  reduction  due  to  viscoelastic 
relaxation  can  be  seen. 


Conclusion : 


Governing  equations  for  the  diffusion  of  small  molecules  in  polymers  have  been  derived  and 
incorporated  into  a  two-dimensional  Finite  Element  code.  In  the  case  of  an  adhesive  joint,  simu¬ 
lations  have  shown  that  the  main  accelerating  effects  in  the  diffusion  rate  of  moisture  through  the 
bond  line  can  be  traced  to  the  strain  and  concentration  dependence  of  the  diffusion  coefficient,  the 
latter  playing  a  predominant  role. 
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Invariably,  the  development  of  powerful  predictive  capabilities  raises  the  question  of  the  ad¬ 
equacy  of  material  characterization.  A  main  finding  of  this  study  is  that  the  bulk  properties  (also 
called  volumetric  properties)  of  the  adhesive  play  a  central  role  in  the  coupling  mechanisms  between 
the  diffusion  behavior  and  the  viscoelastic  response.  Unfortunately,  the  bulk  characterization  of 
polymers  has  not  received  proper  attention  in  experimental  mechanics.  This  trend  must  be  reversed 
if  rigorous  predictions  of  the  hygromechanical  behavior  of  adhesives  are  sought. 
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:  Moisture  concentration 
:  Diffusion  coefficient 
:  A  material  constant 
:  Interpolation  function  for  /"  node 
:  Moisture  concentration  at  f*  node 
:  Coefficient  of  expansion  due  to  moisture  sorption 
:  Coefficient  of  thermal  expansion 
:  Temperature  change 
:  Step-size  for  the  nth  time  step 
:  x  —direction  cosine  for  unit  vector  normal  to 
the  boundary 

:  y  -direction  cosine  for  unit  vector  normal  to  the 
boundary 
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Table  3.  Summary  of  butt  joint  simulations. 


Case  # 

Strain 
dependence 
in  D  &  a 

Concentration 
dependence 
in  D  &  a 

Strain-induced 
driving  force 

Mechanical  BCs 

(longitudinal 

strain) 

1 

No 

No 

No 

2% 

2 

Yes 

No 

No 

2% 

3 

Yes 

Yes* 

No 

2% 

4 

Yes 

Yes* 

Yes 

2% 

(*y  =  0.01) 
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Figure  18.  Effect  of  Penetrant  Size  on  P  and  D:  Trycite  Film  at  50  C,  Subjected  to  a  Strain  of  0.018 
|37j.  (P  and  D  were  measured  one  hour  after  stretching) 
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Figure  19.  Time  Dependence  of  D  in  a  Polystyrene  Film  at  50  C  and  Influence  of  the  Strain 
Level:  Numerical  Simulation  with  BD  =  .25. 
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Figure  21.  Time  Dependence  of  D  in  a  Polystyrene  Film  at  SO  C  and  Influence  of  Physical  Aging; 
Numerical  Simulation  with  B°  »  0.25  and  «  =  .018. 
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Figure  22.  Butt  Joint  Geometry  and  Bond  Line  Discretization. 
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Figure  23.  Simulation  of  diffusion  in  a  Butt  Joint  (Case  1):  Moisture  Profiles  Within  the  Adhesive  in 
the  Absence  of  Coupling. 
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Figure  24.  Simulation  of  Diffusion  in  a  Butt  Joint  (Case  2):  Moisture  Profiles  Within  the  Adhesive 
With  a  Strain  Dependence  in  D. 
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Figure  26.  Comparison  Between  Cases  l,  2,  3  and  4:  Moisture  Profile  Within  the  Adhesive  at 
t  **  60  min. 


3.  A  MODEL  FOR  THE  DIFFUSION  OF  MOISTURE  IN  ADHESIVE  JOINTS.  NUMERICAL 
SIMULATIONS: 


« 


jsesi 


a 


$ 

I 

‘l* 

*•> 


$ 

:: 


3 


CONCENTRATION 

DEPENDENT 

DIFFUSION 


EFFECTIVE  Tg 


CASED 

TRANSPORT 


^  diffusion 


\  CRAZING 

\ 

_  V 


CONCENTRATION 

INDEPENDENDENT 

DIFFUSION 


PENETRANT  ACTIVITY 


Figure  30.  Hopfenberg-Frish  Chart  of  Anomalous  Diffusion:  Effect  of  Penetrant  Activity  on  Transport 
Kinetics  |6J. 
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4.  BOUNDARY  EFFECTS  ON  MOISTURE 


TRANSPORT  IN  ADHESIVES  JOINTS: 


Introduction: 


r 


The  purpose  of  this  section  is  to  examine  interfacial  transport,  namely  diffusion  through  the 
boundary  region  bridging  the  adherend  and  the  adhesive.  It  was  postulated  at  the  beginning  of  our 
discussion  on  diffusion,  that  a  hydrophilic  oxide  could  in  principle  alter  the  free  energy  of  the  water 
molecules  locally.  The  resulting  perturbation  in  the  flux  was  initially  included  in  component  J*  of 
the  flux.  Our  goal  is  now  to  show  that: 

(1)  Although  a  hydrophilic  substrate  gives  rise  to  an  interfacial  pressure,  component  J*  of  the  flux 
can  be  neglected. 

(2)  Recent  experimental  findings  suggest  that  other  types  of  boundary  effects,  which  involve  mo¬ 
bility  changes,  rather  than  a  thermodynamic  driving  force,  are  more  likely. 
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The  Disjoining  Pressure: 


It  will  be  assumed  that  no  cavity  or  crack  is  present  in  the  interphase  region;  that  diffusion, 
not  wicking,  is  the  primary  transport  mode  through  the  interphase.  Further,  it  will  be  assumed  that 
no  contaminant  is  present  in  the  interphase,  excluding  osmotic  effects.  Let  us  consider  a  transition 
region  between  a  high  energy  oxide  layer,  which  is  impermeable  to  water,  and  a  polymer  containing 
diffusing  water  molecules.  The  high  energy  surface  causes  the  water  molecules  to  feel  a  stronger 
intermolecular  field  in  the  vicinity  of  the  substrate  than  in  the  bulk  polymer.  The  resulting  drop 
in  free  energy  is  expected  to  provide  a  driving  force  for  local  perturbations  in  the  diffusional  tlux. 
Penetrant  molecules  such  as  H20  will  tend  to  drift  along  the  direction  which  reduces  their  po¬ 
tential  energy,  that  is  ,  perpendicular  to  the  substrate  surface.  Once  equilibrium  is  reached,  this 
thermodynamic  force  normal  to  the  surface  gives  rise  to  an  intcrfacial  pressure  known  in  the  liter¬ 
ature  as  disjoining  pressure40^12  .  Water  molecules  are  still  attracted  to  the  surface  but  their  mean 
free  path  is  now  much  smaller.  As  a  result,  an  excess  concentration  is  present  locally,  causing  an 
interfacial  pressure.  Thus  the  disjoining  pressure  is  a  measure  of  the  hydrophylic  character  of  the 
substrate.  From  a  thermodynamical  point  of  view,  a  parallel  can  be  made  between  disjoining 
pressure  and  surface  tension: 

•  Surface  tension  arises  because  molecules  feel  a  weaker  intermolecular  field  at  the  boundary  of 
their  diffusion  domain.  It  is  a  measure  of  the  work  required  to  transport  a  molecule  from  the 
bulk  region  to  the  surface.  Surface  tension  tends  to  reduce  the  dimension  of  the  interface  re¬ 
gion. 

•  The  disjoining  pressure,  then  is  just  the  opposite  of  surface  tension:  it  arises  because  molecules 
feel  a  stronger  intermolecular  field  at  the  boundary  of  their  diffusive  domain  then  they  do  in 
the  bulk.  It  is  a  measure  of  the  work  required  to  transport  a  molecule  from  the  boundary  to 
the  bulk  region.  Unlike  surface  tension,  the  disjoining  pressure  tends  to  expand  the  size  of  the 
boundary  region. 
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In  the  following  development,  the  disjoining  pressure  will  be  defined  in  the  framework  of 
equilibrium  thermodynamics;  the  contribution  of  the  surface  potential  to  the  diffusion  flux  will  be 
assessed  and  a  governing  equation  for  interfacial  diffusion  will  be  derived.  The  interphase  will  be 
defined  as  the  boundary  region  where  the  intermolecular  field  of  water  is  different  from  that  inside 
the  bulk  polymer.  Since  the  oxide  is  assumed  to  act  as  a  diffusion  barrier,  it  follows  that  the 
interphase  is  confined  to  the  polymer  side  of  the  interface.  Note  that  this  particular  definition  of 
the  interphase  is  purely  arbitrary.  It  has  value,  only  in  a  context  limited  to  the  next  few  pages. 

Let  superscript  a  denote  a  variable  associated  with  the  interphase.  The  interphase  will  be 
treated  as  a  separate  phase  in  the  thermodynamics  sense.  For  a  small  reversible  change  dU°  in  the 
internal  energy  of  the  interphase,  one  has: 

dU°  =  TdS°  +  EM.rfrtf  ~PdV°  (115) 

where:  T  =  Temperature,  S  =  Entropy,  n  =  Chemical  Potential,  n  =  Number  of  moles,  P  = 
Pressure,  and  V  -  Volume.  Since  S,  n,  V  are  extensive  variables,  one  can  write: 

Ua  =  TSa  +  -  PV°  (1 16) 

Differentiating  (116)  and  subtracting  (115)  yields: 

S°dT  -  V’dP  +  'Ln°dtil  =  0  (117) 

Since  the  water  is  found  in  mixture  with  the  polymer,  the  interphase  can  be  treated  as  a  two- 
component  system.  Thus,  the  free  energy  of  this  system  may  be  written  as: 

Sn ?dnt  =  riffj0  dnH20  + dnR  (118) 

where  subscript  R  stands  for  polymer.  Expression  (117)  is  analogous  to  the  expression  derived  by 
Gibbs  for  a  surface40  ,  except  the  terms  involving  surface  tension  and  surface  area  are  replaced  by 
terms  involving  interphase  pressure  and  interphase  volume  respectively. 

At  constant  temperature,  combining  (117)  and  (118)  will  yield  an  isotherm: 
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Or,  if  n%  is  dose  to  unity,  which  is  valid  when  the  penetrant  is  at  low  concentrations: 


d^R 

^PHiO 


V* 


(120) 


Expression  (120)  indicates  that  a  change  in  the  chemical  potential  of  water  may  cause  a  change  in 


interphase  pressure.  Assuming  that  changes  in  fiH  0  are  mainly  associated  with  the  nature  of  the 


surface  not  with  the  polymer,  it  can  be  stated  that  : 


=  0.  Expression  (120)  then  becomes: 


<a>. 


y<J  d^H20 


(121) 


Note  that  V’ln°Hj0  is  the  partial  molar  volume  of  water,  also  denoted  9Hj0.  In  the  problem  of  in¬ 
terest,  the  absolute  value  of  dP  is  the  disjoining  pressure.  Consequently,  the  disjoining  pressure  is 
given  by: 

n  =  -  AP°  (122) 

The  significance  of  expression  (122)  is  that  the  interphase  behaves  as  though  its  boundaries  were 
repelled  by  a  disjoining  pressure  n.  The  phenomenon  mainly  results  from  three  different  types  of 
interactions:  electrical,  dispersion  and  polar.  The  disjoining  pressure  may  be  written  as  the  sum  of 
three  components41: 

**  =  ^electric  ^dispersion  "b  polar  (123) 


n  is  the  pressure  differential  between  the  bulk  and  the  interphase  due  to  the  presence  of  a  stronger 
intermolecular  field  in  the  interphase: 
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n  =  p 


surface 


-Pi 


bulk 


(124) 


At  this  point,  another  analogy  can  be  made  with  osmotic  pressure:  In  an  osmotic  device,  the 

solvent  is  transferred  to  a  solution  through  a  semi-permeable  membrane,  due  to  the  presence  of  a 

gradient  in  free  energy  between  the  two  contiguous  chambers;  the  net  flow  stops  when  equilibrium 

pressure  A P  =  — -  is  established  in  the  solution  chamber.  In  a  joint,  the  interphase  region 

^h2o 

can  be  seen  as  the  analog  of  the  solution,  and  the  bulk  polymer  region  can  be  seen  as  the  analog 
of  the  semi-permeable  membrane. 

The  chemical  potential  /rHj0  of  water,  is  not  constant  across  the  thickness  of  the  interphase, 
rathe1"  it  is  a  function  of  the  distance  z  to  the  surface.  Noting  that:  %j0  n°„20  =  G°,  we  find  that40: 


n(jr)“"  )A'T'y 


(125) 


where:  A  =  surface  area 

It  becomes  clearer  then,  that  the  disjoining  pressure  is  a  measure  of  the  work  expanded  to  move  a 
molecule  of  water  from  the  surface  to  the  bulk  region: 


J'OO 

A  n(z)  dz  =  Gf  -  Gf  (126) 

o 

where  Gj  and  G%  are  the  free  energies  at  the  surface  and  in  the  bulk  respectively.  G%  —  G’s 
then,  is  a  boundary  free  energy,  very  much  like  a  surface  free  energy  at  the  surface  of  a  liquid. 


Potential  Assisted  Diffusion: 


Consider  a  planar  interface  defined  by  three  perpendicular  axes  Ox,  Oy  and  Ot.  Let  O,  be  the 
out-of-plane  axis  with  the  origin  on  the  interface  plane  and  the  positive  direction  pointing  toward 
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the  polymer.  According  to  the  previous  reasoning,  a  potential  field  will  exist  in  the  vicinity  of  any 
hydrophilic  surface,  in  which  water  molecules  will  “fall”.  The  potential  energy  Wf  of  a  water 
molecule  will  be  a  function  of  z  only.  The  average  drift  velocity  V  of  the  water  molecules  in  the 
presence  of  such  an  attraction  can  be  written  as: 

V  =  —  B\Wp  (127) 


where:  B  =  local  molecular  mobility 


The  molecular  flux  J$  resulting  from  the  potential  field  is  given  by: 

J$  =  -CBVWf 


(128) 


In  the  most  general  case,  the  diffusion  coefficient  is  related  to  mobility  by43: 

D  =  BkTfl  (129) 


The  thermodynamic  factor  n  becomes  important  only  when  the  energy  of  mixing  of  the  system  is 
high.  Again,  this  is  unlikely  in  the  context  of  the  problem  at  hand,  ft  will  be  treated  simply  as  an 
unknown  constant. 

For  the  sake  of  clarity,  the  stress  induced  flux  J*  will  be  omitted  in  the  flux  equation.  This 

simplification,  will  in  no  way  affect  the  final  conclusion.  Using  Cartesian  coordinates  and  noting 
dWt  dW$ 

that  —5 —  =  — - —  =  0,  the  total  flux  becomes: 
dx  dy 


J  = 


i-D 


dC 

dy 


dC  ftC 
dz  kT 


3Wf 

dz 


k 


(130) 


Expressing  the  conservation  of  mass  and  combining  with  the  flux  equation,  yields  the  governing 
equation  for  interfacial  diffusion: 
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A  Functional  Form  for  the  Potential  Energies  Associated 
with  the  Polar ,  Dispersion  and  Electric  Interactions: 


The  following  discussion  will  deal  with  the  force  field  experienced  by  an  isolated  gas  molecule 
placed  in  the  immediate  vicinity  of  a  surface.  This  means  that  the  gas  behaves  ideally,  or  stated 
another  way,  that  an  individual  molecule  does  not  feel  the  intermolecular  field  created  by  the  species 
belonging  to  the  same  phase.  Although  this  assumption  does  not  adequately  represent  the  condi¬ 
tions  experienced  by  penetrant  molecules  mixed  with  a  polymer  matrix,  useful  information  can  be 
obtained  on  the  range  of  the  interfacial  force  field.  The  potential  energy  Wf(z)  represents  the  work 
required  to  move  an  isolated  molecule  from  the  surface  of  the  oxide  to  a  distance  z  from  the  sur¬ 
face: 


dWf  * 

-jf-dz=wf(z)-Wf(o)  (132) 

4> 

(A  random  orientation  is  assumed  for  the  water  molecules) 

It  has  been  mentioned  earlier  that  the  disjoining  pressure  is  the  net  of  several  components: 
dispersion,  polar  and  electrical.  As  a  first  step,  only  the  first  two  interactions,  collectively  known 
as  Van  der  Waal  forces,  will  be  considered.  A  list  of  three  intermolecular  forces  contributing  to  the 
Van  der  Waal  component  of  the  disjoining  pressure  is  given  in  Table  4.  Note  that  the  power  de¬ 
pendency  is  -6  in  all  three  cases.  When  passing  from  a  molecule-molecule  interaction  to  a 
molecule-surface  interaction,  the  total  interaction  energy  becomes  the  sum  of  the  individual  inter¬ 
actions.  Since  oxide  dipoles  (eg:  hydroxyl  groups)  are  scattered  on  the  surface  of  the  oxide, 
dipole-dipole  interactions  are  to  be  integrated  over  the  entire  surface  S: 
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(133) 


where:  C,  =  constant 

The  dispersion  force  exists  between  any  pair  of  molecules  or  atoms  regardless  of  polarity.  Conse 
quently  dispersion  interactions  are  to  be  integrated  over  the  entire  volume  V  of  the  substrate: 


Iff  Tdy 


(134) 


where:  C2  =  constant 

In  performing  the  above  integrations,  the  exponent  on  z  drops  by  two  units  in  the  case  of  dipole- 
dipole  interactions  and  by  three  units  in  the  case  of  dispersion  interactions: 

<  =  -~  (135) 

where  subscript  1  corresponds  to  dipole-  dipole  interactions;  subscript  2  corresponds  to  dispersion 
interactions;  a>,  is  a  constant,  n,  =  4,  and  n,=  3. 

To  avoid  having  an  inf'  .iite  potential  at  the  surface,  which  would  be  meaningless  physically,  it  is 
convenient  to  write  the  potential  Wf  as: 


wl 

( z  +  af‘ 


where  a  may  be  physically  interpreted  as  the  molecular  radius  of  water. 


(136) 
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Let  us  now  evaluate  o,  from  the  disjoining  pressure  FI°  at  the  surface,  that  is,  at  z  =  0.  Ac 
cording  to  expression  (125)  we  have: 


n°— t  "l11*-0  <'37> 

Since  a  variation  in  the  potential  energy  of  the  molecules  corresponds  to  an  equal  variation  in  the 
free  energy  of  the  system15,  we  have: 

dG°  =  n^o  Na  dWp  (138) 

with:  N„  =  Avogadro's  Number. 


And  combining  expressions  (136),  (137)  and  (138)  yields: 


n;V'+! 
Wl~  4  nva 


where  T  is  the  surface  concentration  of  moisture  at  equilibrium. 


(139) 


Finally,  the  potential  energy  can  be  written  as  follows: 


wil  =  -  n?  4 


n,+1 


4  TNa  (a  +  zf1 


+  WP 


where:  Wp°  =  An  arbitrary  integration  constant 

n?  =  polar  component  of  the  disjoining  pressure 
n»  =  dispersion  component  of  the  disjoining  pressure 


The  z  component  of  the  interphase  flux  may  then  be  written  as: 


r4>  _  DOC 

2  ~  kT 


1 

(a  +  z)n'+' 


(140) 


(141) 
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From  equation  (141),  it  becomes  apparent  that  J*  drops  by  97  %,  two  molecular  radii  away  from 
the  surface,  in  the  case  of  polar  forces;  and  drops  by  94  %,  two  molecular  radii  away,  in  the  case 
of  the  dispersion  force.  Note  that  the  distance  z  at  which  J*  vanishes  is  very  close  to  the  atomic 
scale.  In  fact,  within  the  small  boundary  region  where  J*  vanishes,  an  integration  performed  over  a 
continuous  domain  as  in  equations  (133)  and  (134)  yields  incorrect  results.  In  a  region  so  close 
to  the  surface,  a  discrete  summation  of  all  interactions  would  be  necessary.  In  spite  of  this  limita¬ 
tion  of  the  continuum  treatment,  an  important  conclusion  can  be  drawn:  Jf  due  to  Van  der  Waals 
forces  may  be  neglected  for  all  practical  purposes.  Note  that  this  does  not  mean  that  these  inter¬ 
actions  or  the  disjoining  pressure  are  negligible  in  absolute  terms  within  the  domain  size  considered. 
What  is  inferred  from  our  reasoning  is  simply  that  the  gradients  in  the  corresponding  potential  en¬ 
ergies  vanish  at  an  extremely  short  distance,  thereby  causing  negligible  flux  perturbations,  at  least 
as  long  as  the  penetrant  is  not  found  in  a  condensed  phase. 

Let  us  now  examine  the  effect  of  the  electrical  component  of  the  disjoining  pressure.  By 
performing  an  integration  similar  to  that  in  expressions  (133)  and  (134),  it  can  be  shown  that  the 
electric  field  <1>  in  the  vicinity  of  a  uniformly  charged  plate  (same  as  a  plate  of  infinite  dimensions), 
is  independent  of  distance  z  and  is  given  by44  : 


(142) 


where:  a  -  Surface  charge  density 

e0  =  Dielectric  constant  of  the  medium 

Since  water  molecules  behave  as  small  dipoles,  let  us  consider  their  potential  energy  in  an  electric 
field.  The  potential  energy  Wf  of  randomly  oriented  dipoles  placed  in  electric  field  <I>  is  given 
by40 

Wf3  =  -  <n><t>  (143) 


where  <  n  >  is  the  average  dipole  moment  of  the  molecule 
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It  is  reasonable  to  assume  that  a  few  monolayers  away  from  the  surface  plane,  the  dipoles  are 
undergoing  random  thermal  agitation,  or  stated  an  other  way,  that  <  n  >  is  independent  of  z.  If 
the  variation  in  the  dielectric  constant  is  not  too  steep,  this  means  that  Wf  is  almost  constant  and 
its  gradient  in  the  z  direction  very  small.  Again,  the  conclusion  is  that  the  local  flux  perturbation 
is  likely  to  be  minimal. 

A  more  rigorous  treament  of  the  physics  of  fluids  near  solid  surfaces  has  been  given  by 
Teletzke  and  co workers45.  The  equation  of  state  of  a  non-polar  fluid  near  a  solid  was  obtained  by 
summing  all  the  potential  energies  of  the  fluid-solid  interactions  and  fluid-fluid  interactions,  and 
by  allowing  fluid  density  to  vary  as  a  function  of  distance  to  the  solid.  For  closed  isothermal  sys¬ 
tems,  the  equilibrium  density  profiles  and  disjoining  potential  profiles  were  obtained  by  solving  for 
the  minima  of  the  total  free  energy  of  the  system  and  requiring  conservation  of  mass.  The  effect  of 
bulk  density  and  fluid  temperature  on  the  boundary  density  profiles  was  discussed  in  detail  on  the 
basis  of  a  parametric  numerical  analysis.  The  finding  which  is  of  interest  to  us  is  that,  in  all  the 
cases  considered,  any  excess  density  or  any  disjoining  potential  fluctuation  was  found  to  vanish  at 
a  distance  from  the  solid  corresponding  to  10-20  (fluid)  molecular  radii.  Thus,  we  have  some  ad¬ 
ditional  indication  that  the  surface  energy  of  the  substrate  alone  cannot  account  for  flux  perturba¬ 
tions  beyond  a  distance  corresponding  to  a  few  monolayers  of  penetrant. 


Effect  of  Interphase  Polymer  Rheology: 


It  has  been  widely  suggested  that  polymer  conformation  is  likely  to  be  altered  in  the  boundary 
region  near  the  interface.  Recent  results  obtained  from  the  study  of  thin  polysulfone  films  on 
anodized  A1  2090  substrates44  seem  to  indicate  that  segmental  mobility  is  reduced  in  a  region  ex¬ 
tending  as  far  as  a  few  microns  away  from  the  substrate  /  adhesive  interface.  For  example,  a  higher 
than  normal  TJt  as  well  as  variations  in  dielectric  loss  factor,  mechanical  loss  factor,  and  rheological 
activation  energy  have  been  reported  in  the  boundary  region44,47.  In  addition,  these  property 
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changes  have  been  found  to  be  sensitive  to  surface  pretreatment.  Note  that  a  local  drop  in  polymer 
viscosity,  if  it  is  related  to  lower  free  volume,  is  likely  to  decrease  the  mobility  (or  diffusion  coeffi¬ 
cient)  of  penetrant  molecules  near  the  interface  plane.  In  another  study,  Knollman4*  has  measured 
shear  modulus  variations  in  the  interfacial  region  of  a  AF  453  (Epoxy)  /  A1  6061  bond,  using  an 
ultrasonic  wave  technique.  He  showed  that  the  shear  modulus  of  the  epoxy  starts  to  drop  from  its 
bulk  value,  as  far  as  0.3  mm  away  from  the  substrate.  These  various  property  changes  strongly 
suggest  that  the  diffusion  coefficient  in  the  interphase  may  deviate  from  its  bulk  value. 

Changes  in  molecular  mobility  inside  the  interphase  can  be  translated  mathematically  by 
using  a  diffusion  enhancement  factor  D".  With  such  an  approach,  the  interfacial  diffusion  coeffi¬ 
cient  D ’  will  be  related  to  the  bulk  diffusion  coefficient  D  by: 

Da(T,  <r,  C)  =  Da  D[T,  v,C)  ( 144) 

Note  that  the  functional  form  of  the  diffusion  coefficient  is  unchanged,  only  its  magnitude  differs. 
This  simple  approach  has  been  used  successfully  in  metallic  diffusion  studies,  in  order  to  incorpo¬ 
rate  the  effect  of  high  diffusivity  paths  such  as  dislocation  pipes  and  grain  boundaries4’.  It  can  be 
easily  implemented  in  a  finite  element  code  by  assigning  different  values  of  D°  in  various  boundary 
regions.  The  intent  of  this  parameter  is  not  to  give  an  arbitrary  fitting  parameter,  but  rather  to  be 
a  method  to  account  for  molecular  mobility  differences  in  the  polymer  as  a  function  of  proximity 
to  the  boundary.  Variations  in  diffusivity  across  the  transition  region  can  be  represented,  simply 
by  refining  the  Finite  Element  mesh  near  the  boundary  and  by  varying  the  diffusion  enhancement 
factor  in  a  stepwise  manner.  The  main  challenge  obviously,  lies  in  developing  techniques  for 
measuring  D  on  very  thin  films  in  contact  to  a  substrate. 
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Conclusion : 


ft 


Although  the  surface  energy  of  the  substrate  can  give  rise  to  an  interfacial  pressure,  it  cannot 
alone  account  for  flux  perturbations,  mainly  because  the  range  of  the  interfacial  force  field  is  ex¬ 
tremely  short.  Whence,  if  interfacial  diffusion  effects  exist,  they  must  be  related  to  changes  in  the 
molecular  mobility.  A  number  of  recent  findings  indirectly  suggest  such  an  effect,  and  also  indicate 
that  the  thickness  of  the  polymer  layer  affected  can  be  quite  large.  It  is  not  clearly  known  at  this 
time  whether  D ’  is  smaller  or  greater  than  one.  However,  if  future  studies  show  that  D°  is  greater 
than  one  (ie:  lower  diffusivity  locally),  this  will  mean  that  bulk  diffusion  is  rate  controlling,  and  that 
interface  effects  are  not  very  important  in  the  context  of  a  durability  study.  It  is  instructive  to  re¬ 
member  that  we  have  assumed  we  were  dealing  with  a  neatly  fabricated  joint,  namely  that  no 
damage  or  contaminants  were  present  in  the  boundary  region.  These  restrictions  exclude  wicking 
or  osmotic  effects  in  interfacial  moisture  transport. 
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Table  4.  Inter  molecular  Forces  Contributing  to  The  Disjoining  Pressure  |40]. 
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5.  THE  CRITICAL  CONCENTRATION 
MODEL  FOR  JOINT  DURABILITY 
PREDICTION 


Introduction: 

In  Chapters  1,  2  and  3,  a  general  treatment  for  the  diffusion  of  small  molecules  through 
polymers  has  been  proposed,  with  the  objective  of  providing  better  predictions  of  environment 
penetration  into  adhesive  joints.  In  our  treatment  of  diffusion,  no  particular  restriction,  except 
molecular  size,  has  been  placed  on  the  nature  of  the  penetrant  species.  Thus,  in  principle,  we  are 
now  in  possession  of  a  general  predictive  tool  for  calculating  the  rate  of  environment  digression. 

The  second  stage,  in  the  study  of  durability,  involves  the  loss  of  integrity  of  the  interfacial 
regions.  Contrary  to  diffusion,  the  precise  mechanisms  of  joint  degradation  are  extremely  dependent 
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upon  a  particular  adhesive  /  substrate  /  environment  combination.  Environmental  attack  is  known 
to  cause  loss  of  strength  by  four  types  of  mechanisms,  operating  in  an  isolated  manner  or  in 
conjunction5041: 


(1)  By  physically  and  /  or  chemically  attacking  the  adhesive;  e.g:  swelling-induced  loss  of  mechan¬ 
ical  strength,  polymer  hydrolysis. 

(2)  By  causing  a  sign  change  in  the  work  of  adhesion  (bond  line  displacement  by  water). 

(3)  By  breaking  interfacial  covalent  bonds;  e.g:  hydrolysis. 

(4)  By  physically  and  /  or  chemically  degrading  the  substrate;  e.g:  corrosion  or  hydration  of  the 
oxide  layer. 

Therefore,  a  specific  failure  criterion  for  environmental  failure,  as  well  as  specific  test  procedures, 
must  be  developed  for  each  adhesive  /  substrate  /  environment  combination  of  interest.  This  un¬ 
fortunate  state  of  affairs  forces  us  to  limit  ourselves  to  particular  cases  and  to  favor  an  empirical 
approach  over  a  general  theoretical  approach.  Nonetheless,  a  simple  failure  criterion  has  been  dis¬ 
covered  by  Gledhill  et  alp,  which  makes  quantitative  durability  predictions  possible  for  at  least  one 
important  category  of  adhesive  joints:  epoxy-to-steel  joints  in  humid  environments.  The  purpose 
of  this  chapter  is  mainly  to  review  this  Gledhill's  proposal  and  to  show  how  it  can  be  used  in 
conjunction  with  diffusion  data,  to  calculate  residual  strength  as  a  function  of  exposure  time. 
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The  Failure  of  Epoxy-to-steel  Joints  in  High  Humidity 
Environments: 


The  failure  criterion  proposed  by  Gledhill  et  alP  is  based  on  data  collected  earlier  on  the 
behavior  of  mild  steel-epoxy  joints  in  humid  environments53-54.  The  main  results  from  these  previ¬ 
ous  studies  may  be  summarized  as  follows: 

*  AES  and  XPS  shows  that  dry  failure  is  primarily  cohesive53  (close  to  the  oxide  but  not  interfacial). 

*  AES  and  XPS  shows  that,  after  exposure  to  water,  the  locus  of  failure  is  exactly  between  the  metal 
oxide  and  the  adhesive53  (100  %  adhesive). 

*  Substrate  corrosion  is  not  the  cause  of  environmental  failure  but,  rather  a  post  failure 
phenomenon54. 

*  As  long  as  chemisorption  and  interdiffusion  across  the  interface  are  absent,  the  mechanism  for 
environmental  failure  is  bond  line  displacement.  Further,  thermodynamic  considerations  based  on 
surface  free  energy  measurements  can  predict  interface  stability  for  a  number  of  liquids54. 

*  Below  a  certain  relative  humidity  level,  no  degradation  by  moisture  is  observed,  even  after  long 
times. 

These  results  led  Gledhill  et  al.  to  formulate  a  number  of  working  hypotheses  listed  below: 

(1)  There  exists  a  critical  level  of  water  concentration  Cc  in  the  adhesive,  below  which  joint  strength 
remains  intact. 

(2)  At  the  point  where  concentration  C  exceeds  Cc,  there  exists  a  sharp  boundary  between  an  outer 
weakened  region  where  failure  is  interfacial,  and  an  inner  region  where  failure  is  cohesive  in  the 
adhesive. 

(3)  Below  Q.  fracture  strength  is  equal  to  dry  fracture  strength. 
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(4)  The  outer  region,  or  region  of  incipient  interfacial  debonding,  has  negligible  residual  fracture 
strength  and  may  be  assimilated  to  an  environmental  crack. 

According  to  this  type  of  reasoning,  durability  appears  to  be  entirely  governed  by  transport 
kinetics.  Three  parameters  play  a  crucial  role:  (1)  the  rate  diffusion,  (2)  adhesive  bulk  fracture 
strength  and  (3)  the  critical  moisture  concentration.  A  simple  characterization  along  these  lines 
was  undertaken  by  Gledhill  et  al.  .  The  diffusion  coefficient  of  water  in  the  adhesive  (unstressed) 
was  measured  at  several  temperatures  by  sorption.  Bulk  adhesive  toughness  was  determined,  using 
a  compact  tension  specimen.  The  critical  concentration  was  determined  indirectly;  namely,  after 
measuring  fracture  strength  on  an  exposed  butt  joint  specimen,  an  environmental  crack  dimension 
was  calculated  from  simple  fracture  mechanics.  Then,  by  comparing  crack  size  with  the  calculated 
moisture  distribution  at  the  time  conditioning  was  interrupted,  the  critical  concentration  was  de¬ 
duced.  The  value  of  Cc  for  the  steel/epoxy  system  studied  was  found  to  be  1.35  %  in  weight.  At 
this  point,  it  is  useful  to  note  that  Cc  was  obtained  from  combined  fracture  and  diffusion  data  at 
one  temperature  level  only. 

In  order  to  verify  the  model  in  a  conclusive  manner,  residual  butt  joint  strength  versus  ex¬ 
posure  time  was  calculated  at  a  number  of  other  temperature  levels,  and  compared  with  experiment. 
As  demonstrated  by  Figure  31,  fairly  good  agreement  was  found.  Note  that  when  the  saturation 
level  exceeds  the  critical  concentration,  temperature  affects  the  failure  rate  by  changing  the  diffusion 
coefficient  of  water  in  the  adhesive.  When  the  saturation  level  is  below  the  critical  concentration 
(i.e.,  20  C,  55%  RH),  no  debond  is  seen. 

Although  Gledhill  et  al.  ascertained  Cc  implicitely,  a  more  direct  method  can  be  used.  For 
example,  residual  strength  can  be  measured  after  conditioning  several  specimens  under  varying  rel¬ 
ative  humidities,  until  equilibrium  sorption  is  achieved.  After  measuring  solubility,  it  then  becomes 
possible  to  extract  Cc  by  locating  the  knee  in  the  oy  vs.  C  curve.  This  procedure  was  used  by 
Comyn  and  coworkers5S  in  a  study  of  the  durability  of  pretreated  aluminum  /  FM  1000  bonds.  In 
the  next  section,  the  direct  method  will  be  used  to  determine  the  residual  fracture  strength  as  a 
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function  of  moisture  content  in  the  case  of  a  Ultem  1000-to-aluminum  bond.  The  effect  of  surface 
pre-treatment  will  be  studied. 


The  Case  of  Ultem-to- Aluminum  Bonds: 

Specimen  Fabrication  and  Testing: 

Double  cantilever  beam  specimens  with  Aluminum  2024  adherends  and  with  an  Ultem  1000 
bond  line  were  used  in  this  study.  The  specimens  were  fabricated  by  compressing  Ultem  sheets 
between  pretreated  aluminum  beams,  at  a  temperature  of  240  C  (460  F)  and  with  a  pressure  of  7 
kPa  (1000  psi).  Bond  thickness  was  set  to  a  value  of  .254  mm  (.010')  by  using  shims.  The  di¬ 
mensions  of  the  beams  were  3.2  mm  x  15  mm  x  139  mm  (1/8'  x  .59'  x  5.5").  These  dimensions 
were  chosen  to  prevent  plastic  yield  within  the  expected  load  range.  Beam  width  was  calculated 
from  the  diffusion  coefficient  obtained  in  chapter  2  in  order  for  equilibrium  sorption  to  be  attained 
after  one  month  of  exposure  to  humidity.  All  the  beams  were  carefully  polished  with  3M  W2  AN4 
sand  paper  prior  to  chemical  treatment.  Three  types  of  surface  preparations  were  used;  (1)  solvent 
degreasing  only,  (2)  solvent  degreasing  and  an  FPL  etch,  and  (3)  solvent  degreasing  and  a  sodium 
hydroxide  etch  recommended  in  Reference  56  (3  to  10  min.  in  568g/l  DI  water,  rinse  in  DI  water 
and  air  dry). 

All  the  DCB  specimens  were  conditioned  at  50C  and  under  a  controlled  relative  humidity. 
The  principle  of  the  conditioning  fixture  was  exactly  identical  to  that  described  in  Figure  7  of 
Chapter  2,  with  the  obvious  difference  that  no  mass  measurement  was  needed.  A  total  of  8  con¬ 
ditioning  compartments  were  used  simultaneously,  each  with  a  different  saturated  salt  solution  and 
a  different  vapor  activity.  The  list  of  the  various  salts  used  and  the  corresponding  relative  humidities 
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at  50  C  is  given  in  Table  5.  Each  compartment  contained  six  specimens  (two  of  each  surface 
treatment  were  represented).  In  addition,  six  more  specimens  were  placed  in  a  vacuum  oven  at  50 
C  for  the  same  length  of  time,  so  as  to  obtain  fully  dried  control  specimens.  After  40  days  of 
conditioning  residual  Mode  I  fracture  toughness  for  a  given  relative  humidity  was  determined  by 
averaging  ten  measurements  from  two  specimens  (five  measurements  per  specimen).  These  tests 
were  performed  on  an  MTS  hydraulic  machine,  under  a  constant  crosshead  rate  of  .001  inch  per 
second.  For  every  measurement,  crack  length  was  obtained  from  the  beam  compliance  prior  to 
the  onset  of  debonding  and  the  fracture  energy  value  was  computed  from  the  crack  length  and  the 
critical  debond  load  by  using  elementary  beam  theory  (see  chapter  6). 


Results  and  Discussion: 

Upon  visual  examination,  all  the  failure  loci,  including  those  observed  on  dried  specimens 
seemed  to  be  interfacial  between  Ultem  and  the  Aluminum  oxide.  All  the  data  points  from  the 
fracture  study  are  shown  in  Figure  32.  Each  data  point  is  the  average  of  ten  measurements.  The 
averaging  scheme  was  quite  useful  in  that  it  reduced  the  large  scatter  found  in  the  toughness 
measurements.  Against  our  expectations,  Figure  32  does  not  show  any  evidence  of  a  critical  relative 
humidity  for  any  of  the  surface  treatments.  (Note  that  the  critical  relative  humidity  is  related  to  the 
critical  concentration  by  solubility)  According  to  Gledhill,  the  critical  concentration  corresponds  to 
the  transition  between  cohesive  and  interfacial  failure.  It  follows  that  the  absence  of  critical  con¬ 
centration  in  our  results  can  probably  be  attributed  to  the  fact  that  the  bonds  were  not  of  sufficient 
quality  to  exhibit  cohesive  strength  at  low  relative  humidity.  It  also  appears  that  the  Sodium 
Hydroxide  treatment  proposed  in  Reference  56  did  not  lead  to  any  improvement  in  bond  strength 
as  compared  to  a  solvent  degreased  substrate. 
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Conclusion 


Gledhill  and  his  coworkers  have  developed  a  simple  procedure  to  predict  the  failure  of 
epoxy-to-steel  joints  exposed  to  humidity.  The  method  is  based  on  the  critical  concentration  con¬ 
cept.  The  good  agreement  between  theoretical  predictions  and  experiment  indicates  that  the  dif¬ 
fusion  rate  of  moisture  is  the  rate  controlling  parameter  for  durability  in  this  class  of  adhesive  bonds. 
The  characterization  required  for  the  proposed  analysis  includes  the  diffusion  coefficient  of  water 
in  the  adhesive,  the  solubility  of  water  in  the  adhesive,  the  bulk  fracture  strength  of  the  adhesive 
and  the  concentration  of  incipient  failure  Cc.  Unfortunately,  our  experimental  results  on  the  Ultem 
to  Aluminum  system  remain  inconclusive. 
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Table  5.  List  of  the  Salt  Solutions  Used  for  Relative  Humidity  Control. 


SOLUTION 

(SATURATED) 

RELATIVE  HUMIDITY 
AT  50°C  (%) 

MAGNESIUM 

CHLORIDE 

mm 

35 

POTASSIUM 

CARBONATE 

(Ks.C03)2H20 

47 

CALCIUM 

CARBONATE 

Ca(N0a)2 

60 

SODIUM 

NITRITE 

Na(NOz) 

66 

SODIUM 

CHLORIDE 

NaCI 

85 

AMMONIUM 

SULFATE 

<NH4)2S04 

89 

SODIUM 

SULFATE 

(No2SO4)IOH20 

98 

PURE 

WATER 

h2o 

100 
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Figure  31.  Butt  Joint  Fracture  Stress  Versus  Time  in  Environment:  Points  are  Experimental  Meas¬ 
urements,  the  Solid  Line  is  the  Curve  Fit  from  Which  the  Critical  Concentration  Was  De¬ 
termined,  the  Dashed  Lines  Are  Theoretical  Predictions  |52|. 
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6.  DEBOND  PREDICTIONS  IN  THE 


PRESENCE  OF  AN  AGRESSIVE 
ENVIRONMENT:  THE  CASE  OF 
RUBBER-TO-STEEL  JOINTS  IN  SEAWATER. 


Introduction: 


In  Chapters  1  to  5,  a  general  treatment  of  environment-assisted  failure  has  been  proposed.  It 
was  shown  that  there  exists  strong  evidence  that  such  an  approach  can  yield  good  results  for  a  large 
class  of  structural  joints,  when  exposed  to  water  as  the  primary  degradation  agent.  Unfortunately, 
as  we  will  see  with  the  case  of  rubber-to-steel  joints  in  marine  environment,  other  penetrants  than 
water  can  play  an  active  part  in  bond  degradation.  These  penetrants  can  greatly  complicate  the 
analysis  of  the  problem  by  inducing  chemical  reactions  with  the  adhesive,  which  in  turn  can  dras- 
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tically  alter  the  diffusive  process.  In  such  a  case,  novel  approaches  specifically  tailored  to  the 
problem  at  hand  must  be  devised. 

This  chapter  will  illustrate  this  point  by  addressing  the  special  case  of  rubber-to-metal  bonds 
subjected  to  a  cathodic  potential,  in  seawater.  The  failure  mode  will  be  explained,  a  fracture  test 
will  be  described,  and  a  scheme  to  predict  debond  rate  will  be  suggested.  The  emphasis  will  be 
placed  on  the  mechanical  analysis  of  the  fracture  process  in  rubber-to-steel  joints. 


The  Cathodic  Debonding  of  Rubber-to-Steel  Joints  in 
Marine  Environment: 


Elastomer  to  metal  bonding  is  an  important  adhesion  problem  for  a  variety  of  modem 
structures  and  components.  Major  applications  include  the  automotive,  tire,  shipbuilding,  and  off 
shore  drilling  industries.  In  these  applications,  the  rubber  may  serve  to  inhibit  corrosion,  to  seal 
the  component  from  intrusion  of  unwanted  substances,  or  as  a  load  bearing  structure.  Although 
acceptable  bonds  are  routinely  produced  for  many  of  these  functions,  bond  durability  can  be  a 
problem  under  certain  environmental  exposure.  Of  current  concern  are  the  rubber  to  metal  bonds 
for  marine  applications.  Although  these  bonds  have  been  shown  to  be  quite  durable  in  seawater 
environments,  the  presence  of  small  amounts  of  cathodic  potential  can  rapidly  deteriorate  the  ad¬ 
hesive  bond57.  Cathodic  potential  may  be  provided  by  the  galvanic  action  of  dissimilar  metals. 
This  situation  is  intentionally  created  on  many  marine  structures  by  the  use  of  sacrificial  zinc 
anodes.  Ironically,  these  anodes  serve  to  protect  the  steel  hull  but  degrade  the  performance  of 
rubber  bonds,  paint,  and  other  organic  coatings.  Even  if  the  metal  substrates  are  initially  electrically 
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insulated  from  the  anodes,  corrosion  products  and  sediment  can  build  up  on  the  components  to 
permit  the  flow  of  electrons. 

Cathodic  debonding  processes  have  been  widely  studied  by  researchers  in  the  area  of  thin 
organic  coatings58-59.  When  metallic  substrates  are  exposed  to  aqueous  salt  solutions  and  supplied 
with  free  electrons,  reduction  of  dissolved  oxygen  or  water  at  the  adhesive  to  metal  oxide  interface 
occurs  and  locally  increases  the  pH.  When  these  reaction  products  are  sequestered  inside  a  small 
debonded  region  or  adhesive  layer,  an  extremely  high  pH  can  result.  Several  mechanisms  have  been 
proposed  to  explain  the  relatively  rapid  debonding  that  results  from  the  alkaline  environment58-6'. 
Although  the  mechanisms  responsible  for  the  bond  deterioration  are  complex  and  not  entirely  un¬ 
derstood,  our  observations  have  indicated  that  the  high  pH  hydrolyzes  the  polymeric  primer  leaving 
a  weakened  bond.  The  weakened  zone  may  remain  intact  unless  a  sufficient  force  is  applied. 

A  variety  of  tests  have  been  advocated  for  evaluating  the  adhesion  of  a  flexible  layer  bonded 
to  a  rigid  substrate.  The  cone  test  (ASTM  D-429)  has  been  widely  used  for  measuring  the  adhesion 
of  rubber  to  metal  substrates.  The  design  of  this  test  causes  debonding  to  initiate  at  the  interior 
of  the  specimen,  however,  making  it  inappropriate  for  tests  subjected  to  environmental  exposure. 
A  number  of  different  specimens  have  been  based  around  the  peel  concept.  These  include  the 
climbing  drum  peel  test  (ASTM  D-1781),  the  floating  roller  peel  test  (ASTM  D-3167),  and  the 
simple  peel  test,  using  180°  peel  (ASTM  D-903)  or  other  angles.  Unfortunately,  the  fixtures  asso¬ 
ciated  with  these  tests  are  quite  cumbersome  to  use  in  a  harsh  liquid  environment. 

The  various  peel  tests  have  been  utilized  for  measuring  fracture  energies  of  rubber  to  metal 
bonds,  but  they  have  a  number  of  limitations.  All  of  these  tests  tend  to  induce  very  large  defor¬ 
mations  at  the  debond  tips.  These  extreme  distortions  produce  significant  material  and  geometric 
nonlinearities  which  make  the  specimens  very  difficult  to  analyze.  If  one  wishes  to  measure  time 
dependent  fracture  energies,  the  viscoelastic  dissipation  in  the  relatively  large  rubber  bulk  can  alter 
the  results.  Although  fabric  can  be  embedded  in  the  rubber  to  minimize  the  strain  energy 
stored62,  localized  large  deformations  can  still  induce  substantial  energy  dissipation  in  a  nonlinear 
viscoelastic  material.  To  minimize  these  difficulties,  a  specimen  is  needed  in  which  the  defor¬ 
mations  are  smaller  and  in  which  the  viscoelastic  dissipation  is  reduced. 
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In  an  effort  to  avoid  some  of  the  problems  associated  with  these  specimens,  several  improved 
techniques  have  been  developed  and  one  of  these  will  be  reported  herein.  The  objective  of  the 
current  work  has  been  to  develop  tests  to  measure  the  critical  strain  energy  release  rates  for  the 
rubber  to  metal  bonds  in  various  stages  of  deterioration.  A  test  for  determining  Ge  for  elastomer  to 
rigid  adherend  bonds  subjected  to  aqueous  environments  will  be  described.  Although  the  technique 
is  equivalent  to  the  widely  used  Double  Cantilever  Beam  Specimen,  special  problems  and  attributes 
are  important  for  the  extension  of  this  technique  to  evaluate  elastomer  to  meted  bonds. 


The  Double  Cantilever  Sandwich  Beam: 


The  double  cantilever  sandwich  beam  (DCSB)  illustrated  in  Figure  33  is  a  modification  of  the 
double  cantilever  beam  (DCB)  originally  proposed  by  Ripling,  et  a/.“.  The  specimen  is  shown  with 
an  extensometer  mounted  on  the  unit  for  measuring  the  displacements  during  the  calibration  study. 

The  DCSB  specimen  consists  of  a  thin  layer  of  elastomer  bonded  between  two  metal 
adherends.  Although  the  elastomer  appears  much  like  a  thick  adhesive,  the  intent  is  to  study  the 
adhesive  between  the  steel  and  the  rubber.  Thus,  the  elastomeric  layer  is  to  be  seen  as  a  soft 
adherend  in  a  constrained  state.  This  type  of  geometry  greatly  reduces  viscoelastic  dissipation  in 
the  rubber  for  three  reasons:  (1)  the  volume  of  the  elastomer  is  minimal,  (2)  the  strain  levels  are 
smaller  than  in  other  specimens,  thereby  minimizing  the  nonlinear  effects,  and  (3)  the  elastomer  is 
highly  constrained,  forcing  it  to  deform  in  a  bulk  rather  than  shear  mode.  (The  bulk  behavior  of 
polymers  tends  to  be  substantially  less  time  dependent  than  the  shear  behavior.)  In  the  following 
section,  it  will  be  shown  that  the  specimen  can  be  analyzed  with  a  closed  form  beam-on-elastic- 
foundation  solution. 


6.  DEBOND  PREDICTIONS  IN  THE  PRESENCE  OF  AN  AGRESSIVE  ENVIRONMENT:  TIIE 
CASE  OF  RUBBER-TO-STEEL  JOINTS  IN  SEAWATER. 


130 


Design  and  Analysis  of  the  Double  Cantilever  Sandwich 
Beam: 


The  experimental  phase  of  the  work  was  conducted  on  two  different  specimen  geometries 
with  dimensions  given  in  Table  6.  Type  2  is  a  refined  version  of  type  1,  having  thicker  adherends 
and  a  thinner  layer  of  rubber.  The  steel  adherends  were  vapor  degreased,  grit  blasted  with  steel  grit, 
and  vapor  degreased  again  in  tri-chloroethane  1,1,1  in  preparation  for  bonding.  Two  coats  of 
Chemlok  205  primer  (blend  of  chlorinated  rubber  and  phenol  formaldehyde)  and  two  coats  of 
Chemlock  220  topcoat  (proprietary)  were  brushed  on  the  adherends,  allowing  each  to  dry  prior  to 
applying  the  next.  After  preheating  the  mold  in  a  platen  press,  specimens  along  with  uncured 
5109-S  neoprene  rubber  were  inserted.  The  adhesive  and  rubber  were  vulcanized  at  154°C  (310°F) 
for  1.5  hr  at  a  nominal  pressure  of  3.45  MPa  (500  psi). 

In  order  to  characterize  the  toughness  of  adhesive  bonds,  the  relationship  between  the  strain 
energy  release  rate,  load,  and  crack  length  must  be  known.  This  relationship  has  been  predicted 
analytically  and  determined  experimentally.  Using  an  MTS  servo-hydraulic  load  frame,  load- 
deflection  behavior  was  obtained  for  the  two  specimen  configurations  as  a  function  of  crack  length. 
Deflections  were  measured  with  an  extensometer  and  increments  in  crack  length  were  made  by  us¬ 
ing  a  fine-toothed  saw  to  cut  the  rubber.  Although  tension  specimens  of  the  rubber  had  exhibited 
considerable  creep,  the  DCSB  specimens  did  not  show  measurable  time  dependence  presumably 
because  of  the  minimal  amount  of  rubber  loaded  in  a  triaxial  fashion. 

Figure  34  illustrates  the  measured  compliance  of  the  type  1  specimen.  The  dashed  line  re¬ 
presents  the  predicted  compliance  based  on  neglecting  the  deformations  in  the  rubber.  It  is  the 
simple  cantilever  beam  solution  commonly  used  for  the  DCB  specimen  and  assumes  that  the  beam 
rests  on  a  rigid  foundation.  The  compliance  is  based  on  the  relationship 

C  =  8a3  /  (Ebh3)  (145) 
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in  which  a  is  the  crack  length,  £  is  the  Young's  Modulus  of  the  adherend,  b  is  the  specimen  width, 
and  h  is  the  thickness  of  each  metallic  adherend.  Since  the  measured  compliance  is  consistently 
higher  than  the  simple  prediction,  an  improved  approximation  is  required.  To  achieve  this,  the  thin 
rubber  layer  was  assumed  to  act  as  a  classical  elastic  foundation  and  a  closed  form  solution  based 
on  the  finite  length  beam  on  an  elastic  foundation  solution  given  by  Hetenyi64  was  developed. 
K.<tnmnen“  has  employed  a  similar  beam  on  elastic  foundation  solution  for  the  monolithic  double 
cantilever  beam  specimen  to  account  for  transverse  deformations. 

Referring  to  Figure  35,  one  can  calculate  the  total  beam  deflection  at  point  B  as: 

yB  =  yA  +  0A-  a  +  £a3/3£/  (146) 


The  slope  and  deflection  at  point  A  can  be  calculated  by  utilizing  Hetenyi's  solution  for  a  finite 
beam  on  an  elastic  foundation.  For  the  case  of  equal  applied  forces  on  the  two  beams,  symmetry 
requires  that  there  be  no  horizontal  (frictional)  forces  and  the  fracture  mode  is  pure  mode  I,  since 
the  nonsymmetry  due  to  the  rubber  remaining  on  one  adherend  produces  negligible  error. 

The  applied  force  P  at  point  B  is  translated  to  point  A  and  an  equivalent  moment  is  also 
applied  at  point  A.  The  characteristic  root  of  the  governing  differential  equation  for  a  beam  on  an 
elastic  foundation  is: 


(147) 


For  beams  of  finite  length,  the  solution  depends  on  the  nondimensional  quantity  [A(L  —  a)].  For 
the  current  problem,  this  quantity  exceeds  it  and  according  to  Hetenyi,  the  beam  may  be  classified 
as  a  long  finite  length  beam.  For  an  applied  force  P  at  point  A,  the  slope  and  deflection  at  A  are 
given  by: 


*a 


2PX2 

k 


•  1  0 

sinh  Xp  +  sin  Xp 

2  2 

sinh  Xp  —  sin  Xp 


(148) 
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= 


IPX 

k 


sinh  Xp  cosh  Xp  —  sin  Xp  cos  Xp 

2  .2 
sinh  Ap  —  sinMp 


with:  p  =  L  —  a 


(149) 


and  for  an  applied  moment  of  (Pa)  at  point  A,  the  slope  and  deflection  are  given  by: 


eA 


4  PaX3 
k 


sinh  Xp  cosh  Xp  +  sin  Xp  cos  Xp 

•  2  2 
sinh  Xp  —  sin  Xp 


(150) 


y a  = 


2  PaX1 
k 


sinh22p  +  sin2^p 

*  2  o 

sinh  Xp  —  sin  Xp 


(151) 


By  superimposing  these  results  and  substituting  them  into  equation  (146),  one  obtains  a  closed 
form  solution  for  the  deflection  of  the  DCSB  specimen. 

The  remaining  task  for  completing  the  analysis  is  the  determination  of  the  foundation 
stiffness,  k.  This  is  a  difficult  parameter  to  determine  analytically  because  of  the  complex  triaxial 
stress  state  in  the  rubber.  As  a  first  step,  k  was  determined  from  experimental  compliance  data. 
Although  the  calculated  values  of  foundation  stiffness  show  a  slight  decrease  as  the  debond  increases 
in  length  the  values  are  relatively  stable  and  good  agreement  with  experimental  data  is  achieved  as 
illustrated  in  Figure  34.  Again  one  should  note,  however,  that  the  value  of  k  is  chosen  to  give  the 
best  agreement  with  the  compliance  data  rather  than  being  determined  independently.  At  this  stage 
the  foundation  stiffness  may  be  viewed  as  a  curve  fitting  parameter.  If  one  takes  the  observed 
foundation  stiffness  and  multiplies  by  the  thickness  of  the  rubber  layer  and  divides  by  the  width  of 
the  beam,  and  effective  elastic  modulus  of  the  rubber  layer  is  determined.  Values  obtained  are 
presented  in  Table  7  and  as  would  be  expected,  due  to  the  triaxial  constraint,  the  effective  modulus 
of  the  rubber  is  strongly  affected  by  the  thickness  of  the  rubber  layer.  The  measured  values  of  ef¬ 
fective  modulus  lie  between  the  tensile  modulus  of  the  rubber,  4. 14  MPa  (600  psi),  and  the  bulk 
modulus,  876  MPa  (127  ksi),  and  increase  rapidly  as  the  rubber  thickness  decreases. 
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As  a  second  step  in  ascertaining  the  foundation  characteristics,  it  was  shown  that  the  foun¬ 
dation  stiffness  can  be  predicted  analytically  using  a  procedure  developed  by  A.  N.  Gent  et  a/.M>6\ 
The  overall  deformation  of  the  block  of  rubber  is  assumed  to  result  from  the  superposition  of  the 
deformation  of  the  unconstrained  block  and  a  shear  deformation  necessary  to  restore  displacement 
continuity  along  the  bonded  interfaces.  The  normal  force  may  be  written  in  the  form: 


F  —  fcAE'i  (152) 

in  which  A  is  the  cross-sectional  area  of  the  block,  E‘  is  Young's  modulus  of  the  rubber,  e  is  the 

normal  strain  ir.  the  rubber  and  fc  is  a  corrective  factor  accounting  for  the  constraining  effect  of  the 

adherends.  In  the  case  of  a  rectangular  block  of  infinite  length,  A  =  bt,  e  =  and 
4  1 

fe  —  -y  +  -j-  b*jt2,  where  b  is  block  width,  t  is  block  thickness  and  y  is  deflection  in  each  beam. 
Whence,  for  a  foundation  of  unit  width  the  normal  stress  would  be: 

<r  =  ffi'zt  (153) 

or 

®  =  2 fcE'y  (154) 

Now  in  the  theory  of  Beams  on  Elastic  Foundations,  the  foundation  stiffness  k  is  defined  by: 

<t  =  kylb  (155) 

Taking  6=  1  and  comparing  (154)  and  (155)  yields: 

k  =  2/c£'  (156) 

Table  7  shows  that  the  foundation  stiffness  obtained  from  the  above  relationship  and  ob¬ 
tained  previously  by  curve  fitting  the  experimental  data  are  in  reasonable  agreement.  It  should  be 
noted  that  the  solution  to  the  deflection  equation  is  quite  insensitive  to  the  value  used  for  k  .  As 
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a  result,  excellent  agreement  with  the  experimental  compliance  data  is  obtained,  even  with  vari¬ 
ations  of  k  as  large  as  30%. 

It  must  be  noted  that  the  solution  given  by  Gent  for  the  infinitely  long  block  of  finite  width 
is  not  completely  appropriate  for  the  solution  to  the  beam  on  elastic  foundation  because  it:  (1)  does 
not  account  for  the  end  occurring  at  the  debond  tip  and  (2)  does  not  account  for  the  damped 
sinusoidal  displacements.  Both  of  these  factors  should  result  in  an  overprediction  of  the  foundation 
stiffness.  Nonetheless,  the  errors  introduced  do  not  seem  to  strongly  affect  the  calculations  of  the 
strain  energy  release  rate  and  these  approximations  seem  adequate  to  estimate  the  foundation 
stiffness. 

To  further  verify  the  model,  a  strain  gage  was  mounted  on  the  outer  surface  of  one  of  the  steel 
adherends  and  the  strain  was  recorded  for  various  debond  lengths.  The  results  from  this  are  pre¬ 
sented  in  Figure  36.  The  gage  was  positioned  127  mm  (5  in)  from  the  loading  end.  When  the  de¬ 
bond  tip  is  near  end  B,  the  strains  measured  are  very  small.  As  the  debond  approaches  the  gage, 
the  strains  become  positive  where  the  beam  oscillations  induce  a  tensile  mode  in  the  rubber.  As 
the  oscillation  passes  on  ahead  of  the  gage,  the  compressive  strains  induced  by  bending  in  the 
adherends  increase  up  to  a  maximum  value  equal  to  that  predicted  by  simple  beam  theory  for 
M  —  P*  127  mm.  The  predicted  values  of  strain  from  the  beam-on-elastic-foundation  analysis  are 
seen  to  agree  very  well  with  the  experimental  data. 

To  better  understand  how  various  parameters  affect  the  performance  of  the  DCSB, 
parametric  deflection  studies  are  presented  in  Figures  37  to  39.  Figure  37  illustrates  the  effect  of  the 
debond  length  on  the  normalized  deflections  for  the  type  1  specimen.  Figure  38  shows  how 
changing  the  rubber  thickness  would  change  the  deflections  for  a  given  adherend  thickness.  De¬ 
creasing  the  thickness  reduces  the  amount  of  bending  oscillation  and  shortens  the  characteristic 
oscillation  distance,  in  addition  to  reducing  the  time  dependence  as  mentioned  before.  Clearly  the 
rubber  sandwich  should  be  made  as  thin  as  practical.  Figure  39  shows  that  increasing  the  adherend 
thickness  reduces  dramatically  the  oscillations,  in  addition  to  increasing  the  maximum  available 
strain  energy  for  debond  propagation. 
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In  order  to  test  loaded  specimens  in  a  harsh  environment,  it  is  highly  desirable  to  use  a  simple, 
self-loading  device.  The  Boeing  wedge  specimen  is  a  form  of  the  DCB  which  is  loaded  by  driving 
a  wedge  between  the  adherends  to  produce  a  constant  displacement  test.  Unfortunately,  the  avail¬ 
able  G  for  this  type  of  load  is  highly  dependent  on  the  debond  distance,  obeying  an  inverse  fourth 
power  relationship  with  crack  length.  Although  this  type  of  loading  can  be  used,  the  technique  is 
very  sensitive  to  small  errors  in  debond  measurement,  and  the  available  strain  energy  decreases  so 
rapidly  that  the  measurements  must  be  taken  at  very  small  increments  in  debonding  to  obtain  ac¬ 
curate  information.  In  an  effort  to  avoid  these  difficulties,  a  loading  device  which  would  impose  a 
relatively  constant  value  of  G  would  be  beneficial. 

To  investigate  a  constant  G  test,  one  can  write  a  closed  form  solution  for  the  required  force 
necessary  to  produce  a  constant  G  rate: 


pi  _  2Gb 
dC/da 


(157) 


The  expression  for  dC/da  is  given  in  Appendix.  Figure  40  represents  parametric  iso-G  curves  for 
specimen  type  2.  The  indicated  fracture  energies  are  typical  of  Gt  for  the  weakened  bonds.  The 
vertical  portion  of  the  limit  of  validity  domain  represents  a  detached  length  of  216  mm  (8.5  in)  after 
which  the  remaining  support  is  too  short  for  the  analysis  to  be  valid.  The  curved  portion  of  the 
limiting  domain  represents  rotations  of  the  loaded  end  of  the  beam  which  exceed  8  degrees  and 
may  result  in  beam  foreshortening  errors. 

Iso-G  loading  can  be  accomplished  quite  simply  by  using  a  spring  to  load  the  specimen. 
Figure  41  illustrates  the  loaded  fixture  in  place  in  the  testing  environment.  Simple  helical  com¬ 
pression  springs  were  used  to  provide  the  energy  to  the  system.  Figure  42  illustrates  the  force 
available  from  a  given  spring  for  several  different  preload  rates  superimposed  on  a  desired  Iso-G 
curve.  By  properly  selecting  the  spring  stiffness  and  the  preload,  one  can  obtain  an  Iso-G  test 
window  with  a  width  of  50  mm  (2  in).  Thus  one  is  able  to  set  up  the  loading  device  to  provide  a 
relatively  constant  debond  driving  force  over  a  relatively  long  debond  distance.  With  a  small  as¬ 
sortment  of  springs  and  with  judicious  choices  of  the  preload,  one  can  follow  a  required  G  curve 
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over  most  of  the  length  of  the  specimen.  A  compact  load  cell  and  an  extensometer  can  be  used  to 
monitor  the  actual  applied  G  to  correct  for  small  deviations.  By  holding  G  relatively  constant, 
however,  the  collection  of  debond  growth  data  is  greatly  simplified. 


Specimen  Design  Guidelines: 


In  designing  DCSB  specimens,  the  following  guidelines  may  prove  helpful: 

•  Make  the  soft  layer  as  thin  as  practical  to  minimize  the  viscoelastic  effects. 

•  Select  a  specimen  width  which  is  large  compared  to  the  elastomer  thickness  in  order  to  increase 
the  constraining  action  on  the  elastomer. 

•  Determine  the  appropriate  G  range  for  the  adhesive  system  and  select  a  beam  thickness  which 
can  provide  the  desired  G  values  while  remaining  in  the  small  deflection  envelope. 

•  The  remaining  bonded  length  should  be  large  enough  that  the  long  finite  length  beam  as¬ 
sumption  is  valid. 

•  For  the  Iso-G  loading,  superimpose  the  spring  load  decay  curves  on  the  Iso-G  parametric 
curves  and  select  an  appropriate  initial  load  and  spring  constant  to  give  the  desired  result. 
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Difficulties  with  Interaction  Between  Diffusion  and 
Fracture: 


Because  the  specimen  was  designed  to  measure  the  fracture  energy  of  debonding  under  con¬ 
ditions  of  cathodic  delamination,  the  adherends  were  not  thick  enough  to  permit  testing  of  the  dry 
bond.  Preliminary  tests  were  conducted  in  seawater  or  NaOH  solutions  with  cathodic  potential 
applied.  Slow  debond  rates  were  measured  but  difficulties  were  encountered  because  the  results 
were  not  consistent.  When  specimens  were  broken  open  after  testing,  a  significant  chevron  or  re¬ 
verse  tunneling  effect  was  noted.  The  reasons  for  this  phenomena,  which  is  opposite  of  that  expe¬ 
rienced  in  standard  DCB  specimens,  is  that  the  peel  stresses  are  higher  at  the  outside  of  the 
specimen  because  of  the  "poker  chip'  effect,  and  because  diffusion  occurs  from  these  edges  into  the 
specimen.  Because  the  chevron  effect  was  so  severe,  meaningful  debond  distances  could  not  be 
obtained. 

As  other  tests  on  the  neoprene  to  steel  bond  were  conducted,  it  became  apparent  that 
moisture  diffused  into  the  specimen  and  left  the  bond  intact  but  severely  weakened.  When  choosing 
between  adhesive  systems  for  applications  exposed  to  cathodic  potential,  the  two  most  important 
parameters  to  measure  appear  to  be  the  rate  at  which  diffusion  occurs  and  the  retained  strength  of 
the  weakened  bond.  Although  the  DCSB  appeared  to  be  limited  for  measuring  fracture  parameters 
while  diffusion  into  the  specimen  was  occurring,  it  was  decided  that  the  important  fracture  param¬ 
eter  to  measure  was  the  critical  strain  energy  release  rate  of  the  weakened  bond.  Several  specimens 
were  conditioned  in  environment  to  produce  a  weakened,  but  intact,  cathodically  degraded  bond. 
By  imposing  potential  to  only  one  adherend,  this  bond  was  preferentially  weakened.  Once  the 
specimens  were  conditioned  they  were  tested  in  a  tension  test  machine.  This  weakened  bond 
strength  was  found  to  be  approximately  .054  J/m2.  This  fracture  energy  is  four  orders  of  magnitude 
smaller  than  the  60  J/m2  value  corresponding  to  the  dry  strength  of  the  joint  (cohesive  in  rubber). 
Estimates  for  the  weakened  bond  strength  were  also  obtained  by  an  alternate  approach6*  and  were 
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found  to  be  within  20%.  This  is  considered  to  be  quite  good  agreement  because  of  the  difficulty 
in  uniform  fabrication  and  conditioning  of  specimens,  and  the  extremely  small  values  of  Gc. 


A  Technique  for  Predicting  Debond  in  Mechanically 
Loaded  joints: 


The  residual  fracture  strength  data  discussed  in  the  previous  section  were  obtained  after  con¬ 
ditioning  specimens  in  the  absence  of  an  external  mechanical  load.  Thus,  this  type  of  information 
is  inherently  insufficient  for  predicting  debond  rate  in  actual  service  conditions  where  stress  is 
present.  The  reason  is  that  depending  on  the  applied  stress  level,  the  debond  tip  may  coincide  with 
the  advancing  weakened  bond,  or  may  lag  behind.  If  the  applied  stress  is  high  enough,  debond  can 
proceed  faster  than  the  diffusion  process  (with  the  Chemlok  205/220  system,  failure  is  cohesive 
when  the  bond  has  not  been  weakened).  In  the  absence  of  applied  or  residual  stress,  diffusion  of 
the  environment  still  occurs,  although  adherend  separation  may  not  occur.  At  intermediate  stress 
levels,  the  debond  tip  coincides  with  the  advancing  weakened  bond;  stress  opens  the  weakened  re¬ 
gion,  allowing  convection  to  transport  species,  rather  than  the  slower  diffusion  process. 

Whence  it  appears  that  the  driving  force  for  the  opening  mode  of  fracture  (mode  I),  because 
it  governs  the  relative  position  of  the  debond  tip  with  respect  to  the  diffusion  front,  must  be  the  rate 
determining  parameter.  We  conclude  that  a  crucial  piece  of  information  in  a  durability  study  is  the 
debond  rate  versus  G,  curve  (This  kind  of  dynamic  fracture  data  may  be  obtained  with  the  DSCB 
specimen,  or  with  a  variety  of  other  specimens  based  on  the  blister  concept).  A  complete  debond 
rate  versus  G,  curve  offers  the  advantage  of  scanning  the  entire  range  of  rate  controlling  situations 
in  a  purely  empirical  way.  The  rate  of  debond  of  a  real  life  structure  can  then  be  calculated  simply 
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by  iteratively  computing  the  current  mode  I  component  of  the  energy  release  rate  and  updating  the 
debond  tip  position,  using  an  experimental  data  base  obtained  by  the  techniques  described  herein. 


Conclusion: 


The  DCSB  offers  distinct  advantages  over  other  specimens  for  measuring  the  critical  strain 
energy  release  rate  for  elastomers  bonded  to  a  rigid  adherend  because  it  minimizes  the  viscoelastic 
dissipation.  It  is  easily  fabricated  and  can  be  analyzed  with  a  closed  form  beam  on  elastic  founda¬ 
tion  solution.  Foundation  stiffness  can  be  predicted  with  a  rubber  elasticity  analysis.  A  simple 
spring  loading  device  can  be  used  to  provide  near  constant  G  loading  over  a  relatively  wide  test 
window.  When  using  the  specimen  in  an  aqueous  environment,  diffusion  at  the  sides  can  cause 
spurious  results.  For  these  cases  however,  specimens  conditioned  to  equilibrium  can  be  tested  with 
good  results.  Finally,  a  simple  procedure  is  available,  for  using  experimental  debond  data  to  predict 
the  service  life  of  joint  structures,  even  when  the  precise  degradation  mechanisms  are  not  known. 
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Table  6.  DCSB  Specimen  Dimensions 


• 

Specimen 

Type 

Bond 

Thickness 

Beam 

Length 

Adherend 

Thickness 

Beam 

Width 

#1 

1.27  mm 
(0.05  in) 

273  mm 
(10.75  in) 

3.18  mm 
(0.125  in) 

25.4  mm 
(1.00  in) 

• 

#2 

0.76  mm 
(0.0a  in) 

267  mm 
(10.5  in) 

6.35  mm 
(0.250  in) 

25.4  mm 
(1.00  in) 

Table  7.  Foundation  Stiffness  and  Effective  Modulus  of  the  Rubber  Layer 


Specimen 

Type 

Bond 

Thickness 

Foundation 
Stiffness 
k  =  2f'E 
(Predicted) 

Foundation 

Stiffness 

k 

(Experimental) 

Effective 

Modulus 

khlb 

(Experimental) 

#1 

1.27  mm 
(0.05  in) 

1.11  MPa 
(162  ksi) 

0.68  MPa 
(98.9  ksi) 

34.4  MPa 
(5.00  ksi) 

#2 

0.76  mm 
(0.03  in) 

3.07  MPa 
(446  ksi) 

3.45  MPa 
(500  ksi) 

103.5  MPa 
( 1 5.0  ksi) 
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The  compliance  of  a  DCSB  specimen  is  given  by:  C  =  2 yBIP 


£  _  1  sinh  2p  cosh  Xp  —  sinh  Xp  cos  Xp  a  sinh Xp  +  sin  Xp 

ED?  (sinh2,lp  —  sin22p)  EIX1  ( sinh2ip  —  sin2Ap) 


2a2  sinh  Xp  cosh  Xp  4-  sin  Xp  cos  Xp  2a3 
(sinh2Ap  —  sin2)p)  ^£7 


where:  p  =  (L-  a) 


The  derivative  of  the  compliance  with  respect  to  debond  length  is  given  by: 


¥■'»'  with  D  —  i(L  —  a) 


“  _  (L  ~  DIX)1  H - y-  {( cosh2D  +  sinh2/)  —  cos 2D  +  sin 2D)(  sinh2/)  —  sin2) 


-2 (SinhD  CoshD  —  sinD  cosD)2}/(  sinh2/)  —  sin2/))2 


H - y  {[2(S7nAD  CoshD  +  sinD  cosD)(L  —  DjX)  —  1/A(  sinh2D  +  sin2/))](  sinh2D  —  sin 2/>) 

EIX 


—  2(  sinh2/)  4-  sm2D)(L  —  D)X.)(SinhD  CoshD  —  sinD  cosD )}  /  ( sinh2/)  —  sin2/)) 


2  rv  rv\2 


+  yyy  {[( coshz£>  +  sinh2D  4-  cos2D  -  sin2 D)(L  -  DjX) 
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-  2jX(L  -  Djl)(SinhD  CoshD  +  sinD  cosD )]  (  sinh2D  —  sin^D) 


-2 (L  -  DIX)\  sinh2/)  cosh2/)  -  sin2/)  cos2/))}/ ( sinh2/)  -  sin2/))2 
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Figure  34.  Experimental  and  Predicted  Compliance  vs.  Debond  Distance. 
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Figure  36.  Comparison  of  Experimental  and  Predicted  Surface  Strain:  Normalized  Surface  Strain  at 
Position  x  *•  127  mm  (5  in)  as  a  Function  of  Debond  Length  (Specimen  2). 
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Figure  37.  Predicted  Normalized  Deflection*  Tor  Different  Debond  Lengths. 
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NORMALIZED  DEFLECTION- y/P  (in/klp) 
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IOO  200 

POSITION -X  (mm) 


Figure  39.  Predicted  Normalized  Deflections  for  Various  Beam  Thicknesses. 
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In  the  case  of  polymer-to-metal  bonds,  water  diffusion  through  the  adhesive  is  a  rate  con¬ 
trolling  factor  for  durability.  The  diffusion  of  small  molecules  such  as  water  in  a  polymer  can  be 
understood  at  the  molecular  level  by  invoking  the  concept  of  free  volume  already  used  extensively 
in  the  field  of  polymer  rheology.  The  free  volume  treatment  of  diffusion  in  polymers  is  very  helpful 
in  understanding  the  effect  of  temperature,  mechanical  deformation  and  physical  aging  on 
diffusivity  and  solubility.  In  addition,  by  recognizing  that  the  free  volume  controls  both  the  rate 
of  segmental  rearrangement  and  the  transport  rate  of  small  foreign  molecules,  diffusion  and 
viscoelasticity  can  be  unified  into  one  single  theoretical  treatment.  This  unified  continuum  approach 
greatly  facilitates  the  modelling  of  coupling  effects  between  diffusion  and  stress  relaxation.  Further, 
we  found  that  the  observed  decrease  in  the  durability  of  metal-to-polymer  bonds  in  humid  envi¬ 
ronments  due  to  increased  temperature  and  stress  can  be  explained  by  the  accelerated  transport  rate 
of  water  under  these  conditions.  Boundary  diffusion  is  also  a  possible  path  for  water  transport  in 
adhesive  joints.  Theoretical  considerations  as  well  as  recent  experimental  results  show  that  flux 
perturbations  near  the  boundaries  of  the  adhesive  layer  are  likely  to  be  related  to  local  fluctuations 
in  polymer  properties,  not  surface  energetics. 

Another  key  requirement  for  improved  durability  prediction  is  the  development  of  failure 
criteria  for  interfacial  failure  in  the  presence  of  water  and  stress.  Existing  models  are  still  very  em- 
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pineal  and  limited  to  specific  adherend-interphase-polymer  systems.  In  the  case  of  mild  steel-to- 
epoxy  bonds,  the  critical  concentration  concept  seems  to  yield  very  promising  results. 

Finally,  a  convenient  and  well  understood  specimen  geometry  to  study  interfacial  debonding 
experimentally  is  the  Double  Sandwich  Cantilever  Beam  specimen.  The  mechanical  analysis  and 
test  procedure  for  the  DSCB  specimen  has  been  adapted  to  the  study  of  the  interfacial  debonding 
of  rubber-to-steel  joints  in  harsh  environments.  The  soft  polymeric  layer  increases  the  overall 
compliance  of  the  specimen  and  necessitates  a  beam-on-elastic-foundation-type  of  analysis. 

The  following  recommendations  for  joint  design  can  be  made  from  our  study: 

(1)  Select  an  adhesive  with  low  water  solubility  to  prevent  water  from  reaching  the  interface  in 
sufficient  concentration.  Another  beneficial  effect  of  low  solubility  is  to  minimize  swelling  stresses 
and  their  adverse  effect  on  diffusion  rate  and  interface  fracture. 

(2)  If  recommendation  (1)  cannot  be  met,  select  an  adhesive  with  a  low  water  diffusivity  so  as  to 
slow  down  the  sorption  rate. 

(3)  Minimize  peel  stresses  in  order  to  lower  the  diffusion  coefficient  and  slow  down  mode  I  crack 
growth,  which  is  extremely  detrimental  to  durability. 

(4)  Promote  stable,  receptive  oxides  capable  of  developing  strong  chemical  bonds  with  the  polymer. 
Such  bonds  are  not  easily  displaced  by  the  physical  intermolecular  forces  acting  around  water 
molecules. 
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